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given  to  those  methods  most  applicable  to  ion  collisions  with  silane  (SiH4)  and  to  those  molecules  where  estimates 
of  the  general  applicability  and  accuracy  of  the  methods  could  be  defined. 

During  Phase  II  of  our  research  program,  detailed  quantum  mechanical  calculations  of  the  potential  energy 
surfaces  for  prototype  ion-molecule  reactions  were  carried  out.  The  reactions  that  were  examined  included 
dissociative  charge  transfer  for  SiH4  interactions  with  H"*",  H3'*',  He'*'  and  Ar'*'.  The  minimum  energy  reaction 
pathways  that  were  obtained  from  these  ab  initio  studies  were  then  used  as  the  input  to  our  R-matrix  kinetics  code, 
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carried  out.  In  addition,  a  careful  thermodynamic  analysis  of  SiH4‘  and  SiFjj'  anions  was  performed  !r>  examine  the 
role  of  these  species,  which  may  have  long  residence  times  in  gaseous  discharges. 

This  investigation  was  carried  out  using,  mainly,  ab  initio  quantum  mechanical  methods.  Computational 
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This  report  was  prepared  by  the  United  Technologies  Research  Center,  East  Hartford, 
Connecticut,  under  Contract  F33615-87-C-2718.  The  research  was  funded  by  the  Air  Force 
Wright  Labortory,  Wright  Patterson  Air  Force  Base,  Ohio  45433-7650. 

Inclusive  dates  of  research  were  10  June  1987  through  28  March  1993.  Captain  Janies  R. 
Shoemaker  and  Dr.  Alan  Garscadden  have  served  as  Project  Managers  for  this  contract. 

Very  useful  discussions  with  Captain  Peter  Haaland  (Air  Force  Institute  of  Technology)  and 
Professor  Mark  Gordon  (Iowa  State  University)  are  also  acknowledged. 

All  aspects  of  the  research  work  reported  herein  were  aided  by  the  skilled  help  of  Judith  B. 
Addison  (UTRC),  who  assisted  in  the  analysis  of  the  calculated  data  and  in  the  preparation  of  the 
final  report. 

The  objectives  of  this  research  program  were:  1)  to  examine  and  expand  existing 
computational  techniques  that  are  applicable  to  studies  of  the  kinetics  of  ion-molecule  reactions,  2) 
to  compute  the  potential  energy  surfaces  for  several  prototype  ion-molecule  reactions  which  include 
electron,  cation  and  anion  collisions  with  silane,  3)  to  predict  the  reaction  rates,  their  energy 
dependence  and  product  distribution  for  these  prototype  systems,  4)  to  examine  the  thermodynamic 
database  for  SiHn-,  SiFn-  and  Si2Hn-  anions,  and  5)  to  study  possible  sequential  clustering 
reactions  of  silane  anions.  This  program  was  directed  toward  improving  the  technical  information 
base  in  support  of  the  modeling  of  silane  glow  discharges. 
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SECTION  1 


INTRODUCTION 

The  deposition  and  characterization  of  thin  rilms  have  received  enormous  attention  owing  to 
their  importance  in  the  fabrication  of  microelectronic  devices  and  solar  cells.  Considerable 
attention  has  been  given  to  silicon  films  deposited  on  quartz  or  other  substrates  and,  more  reccntiy, 
to  in  situ  silicon  doping  of  compound  semiconductors  such  as  GaAs.  One  process  for  film 
deposition  is  to  induce  a  gas  phase  reaction  by  either  pyrolytic  or  photolytic  means,  the  film  growth 
being  determined  mainly  by  gas  phase  decomposition  kinetics.  Such  chemical  vapor  deposition 
(CVD)  processes  are  often  slow  and  inefficient,  particularly  if  laser  irradiation  is  required  to  induce 
the  decomposition  reactions. 

Another  route  toward  thin  film  production  involves  the  plasma  enhanced  deposition  of  silicon 
(Reference  1)  from  a  gas  mixture  containing  a  silicon-bearing  molecule  such  as  silane  (SiH4), 
disilane  (Si2H^.  or  a  halogen-substituted  silicon  compound  such  as  silicon  tetrafluoride  (SiF4).  In 
particular,  glow  discharges  of  silane  or  silane-noble  gas  mixtures  have  attracted  considerable 
interest  as  a  means  of  producing  device-quality  hydrogenated  antorphous  silicon  (a-Si:H)  films 
(Reference  2).  The  growth  of  such  films  depends  on  the  complicated  chemistry  that  is  occurring  in 
the  cathode  sheath  region  of  the  typical  capacitively  coupled  RF  discharge.  A  good  deal  of 
progress  has  been  made  over  the  past  few  years  in  understanding  the  electron  energy  distribution 
and  the  primary  collision  processes  that  occur  in  the  sheath  region  (References  3  -  5).  However,  a 
quantitative  description  of  all  the  major  reactions  that  occur  in  the  discharge  region,  including 
electron  impact  dissociation  and  ionization,  and  the  chemical  role  of  both  cations  and  anions 
(formed  by  dissociative  attachment)  is  required  to  permit  reliable  modeling  of  RF  discharge 
sheaths.  Goals  of  this  research  program  were  to  identify  the  ion  chemistry  occurring  in  silane 
discharges,  to  establish  a  reliable  thermochemical  database  for  silicon-containing  anions  and  to 
predict  kinetic  reaction  rates  and  product  distributions.  Emphasis  was  placed  on  the  development 
of  sound  theoretical  methods  to  predict  the  collisional  energy  dependence  of  the  reaction  rates  and 
to  examine  the  kinetics  of  ion-molecule  reactions  and  electron  attachment  processes. 

Several  qualitative  pictures  of  both  ion  and  neutral  fragment  production  in  silane  discharges 
have  been  developed.  Haller  (Reference  6)  has  analyzed  the  current  status  of  ion  species  of  the 
type,  Si/H^'*’,  using  mass  spectrometric  techniques.  Secondary  ion-silane  reactions  to  form 
disilane  ions  were  proposed  to  explain  the  observed  ion  currents.  A  separate  model  was  proposed 
by  Turban,  et  al.  (Reference  7),  in  which  primary  silane  ions  are  formed  by  low-energy  electron 
impact  dissociation,  followed  by  secondary  ion-molecule  reactions.  Chapman  and  Gallagher 
(Reference  5),  however,  have  criticized  the  applicability  of  this  model  to  cathode  sheath  conditions 
where  it  is  believed  that  high  energy  electron  processes  (10-100  eV)  govern  the  overall  reaction 
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kinetics.  Their  interpretation  of  the  discharge  ion  chemistry  is  that  electron  impact  dissociation  of 
silane  occurs  as  a  primary  process  forming  SiHn'*’  (n=0,3),  with  n=2  or  3  heavily  favored  over 
more  highly  stripped  silane  Neutral  SiH^  radical  production  as  a  primary  process  must  also  be 
considered.  Based  on  'he  known  rate  of  H-atom  abstraction  from  silane  and  the  direct  electron 
impact  dissociatir  icaction,  SiH3  radicals  are  believed  to  predominate  in  the  discharge  region. 
SiH2''‘  and  more  highly  stripped  silane  cations  are  probably  driven  to  SiH3  by  fast  reactions  of  the 
type  +  SiH4  — >  SiHn+i  +  SiH3'*'. 

The  largest  uncertainties  at  present,  in  the  chemistry  of  the  discharge  sheath  region,  are  with 
the  secondary  ion-molecule  reactions  occurring  in  pure  silane-noble  gas  mixtures.  The  possible 
silane  ion-molecule  reactions  are  well  known  (Reference  8),  but  in  most  cases,  the  rate  coefficients 
are  known  experimentally  only  in  the  low-energy  region  (<  1  eV)  or  have  been  estimated  using  the 
Langevin-Gioumousis-Stevenson  formula  for  the  collision  cross  sections.  This  simple  model 
predicts  that  the  ion-molecule  reaction  rate  coefficients  are  independent  of  collision  energy. 
Dissociative  charge  transfer  rates  for  collisions  between  SiMn"*"  ions  and  silane  are  mainly 
unknown  at  present  since  knowledge  of  the  detailed  mechanism  of  the  ion-molecule  collision  is 
required  to  determine  branching  into  the  dissociation  channels.  The  somewhat  simpler  reactions 
between  silane  and  H+,  H3'‘‘  or  the  noble  gas  ions,  which  are  more  amenable  to  theoretical 
analysis,  are  also  mainly  unknown  at  present  in  a  mechanistic  sense.  In  addition,  a  preliminary 
analysis  of  the  role  of  silicon  hydride  anions,  SiHn",  has  recently  been  reported  (Reference  4). 

This  study  suggested  that  electron  dissociative  attachment  to  silane  may  exhibit  a  large  cross 
section,  particularly  if  the  silane  molecule  is  vibrationally  or  rotationally  excited.  The  role  of  anion 
chemistry  in  describing  the  kinetic  processes  in  the  sheath  region  has  received  little  attention  to 
date.  As  an  example,  Garscadden  (Reference  9)  and  Haaland  (Reference  10)  have  discussed  the 
possible  role  of  silane  anions  in  plasma  assisted  CVD  devices.  A  significant  production  rate  for 
SiH„"  would  greatly  influence  our  understanding  of  the  overall  chemistry  since  fast  mutual 
neutralization  reactions  between  cations  and  atrions  could  occur. 

The  intent  of  this  technical  program  was  to  address  these  uncertainties  in  the  chemistry  of 
silane  ion-molecule  reactions  through  a  theoretical  analysis  of  the  governing  reaction  rates.  Since 
dissociative  reactions,  with  or  without  accompanying  charge  transfer,  may  occur,  our  approach 
was  based  on  detailed  quantum  mechanical  studies  of  the  reaction  surfaces  to  yield  all  possible 
reaction  product  channels.  Subsequent  to  this  analysis  of  the  reaction  energetics,  cross  sections  for 
the  various  ion-molecule  processes  were  calculated  within  a  quantum  framework.  During  the 
course  of  this  research  investigation,  several  new  areas  for  research  became  apparent.  These  areas 
mainly  involved  the  formation  and  reactions  of  negatively  charged  species  whose  chemistry  and 
relative  importance  in  silane  plasmas  is  uncertain.  In  contrast  to  positive  ions,  which  are  swept  to 
the  walls  by  the  ambipolar  fields  in  the  discharge,  negative  ions  tend  to  become  trapped  with 
resultant  long  residence  times.  Garscadden  (Reference  1)  has  suggested  that  negative  ions  may 
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form  clusters  through  a  series  of  ion-molecule  reactions  similar  to  that  discussed  by  Mandich, 
Reents  and  Jarrold  (Reference  1 1)  and  Raghavachari  (Reference  12),  for  positive  ion  cluster 
growth.  It  was  not  clear  whether  sequential  clustering  reactions  of  negative  silane  ions  would 
encounter  a  bottleneck  at  a  small  cluster  size  or  whether  the  anion  chemistry  is  sufficiently  different 
that  growth  continues  to  the  critical  sizes  needed  for  spontaneous  formation  of  large  silicon 
particulates.  In  addition,  the  thermochemistry  of  silicon-containing  negative  ions  was  poorly 
defined  and  was  identified  as  an  area  for  additional  research  under  the  contract.  These  areas  and 
further  studies  of  the  formation  and  decomposition  pathways  for  negative  ion-molecule  reactions 
such  as  H~  -t-  SiH4  and  F"  -i-  SiF4  were  considered  as  part  of  this  research  program  under  a 
contract  extension. 

The  general  composition  of  this  report  is  as  follows.  In  Section  2,  we  present  a  description 
of  the  mathematical  methods  which  were  employed  in  this  research.  Included  in  Section  2  are  sub¬ 
sections  which  deal  with  the  construction  of  electronic  wavefunctions,  the  calculations  of 
expectation  properties,  the  calculation  of  electronic  wavefunctions  using  ab  initio  methods,  a 
description  of  the  development  of  R-matrix  codes,  and  a  description  of  several  electronic  structure 
codes  that  were  used  during  the  course  of  this  research  program.  The  calculated  results  and 
pertinent  discussion  are  presented  in  Section  3.  The  first  part  of  this  research  effort  focussed  on 
the  energy  dependence  and  reaction  product  distributions  for  dissociation  charge  of  Silij  with  H'*’, 
H3'*‘,  He"*"  and  Ar"^.  The  extensive  detailed  calculations  and  discussion  of  this  research  are 
contained  in  Report  WL-TR-91-2022  and  are  not  duplicated  herein.  A  concluding  summary  of  the 
current  program  is  given  in  Section  4.  A  list  of  all  technical  presentations  and  publications  that 
have  resulted  from  this  research  program  arc  given  in  Appendix  A.  Abstracts  of  the  pertinent 
published  papers,  reports  and  technical  talks  are  collected  in  Appendix  B. 
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SECTION  2 


MATHEMATICAL  BACKGROUND 

Central  to  these  theoretical  studies  are  the  actual  quantum-mechanical  calculations  which 
must  be  carried  out  for  the  atomic  and  molecular  species.  In  particular,  the  methodology  for 
computing  the  minimum  energy  reaction  pathways  and  testing  for  vibrational  stability  must  be 
carefully  analyzed.  For  added  clarity,  various  aspects  of  these  calculations  are  discussed  below. 

Much  evidence  on  diatomic  and  polyatomic  systems  indicates  the  inadequacy  of  minimum 
basis  sets  for  constructing  quantitatively  conect  molecular  wavefiinctions  (References  13  and  14). 
This  means  inner-shell  and  valence-shell  orbitals  of  quantum  numbers  appropriate  to  the  atoms 
(Is,  2s,  2p,  for  C,  N,  O;  etc.).  The  main  deficiency  of  the  minimum  basis  set  is  its  inability  to 
properly  describe  polarization  and  the  change  of  orbital  shape  for  systems  which  exhibit  large 
charge  transfer  effects.  Values  of  the  screening  parameters  for  each  orbital  can  either  be  set  from 
atomic  studies  or  optimized  in  the  molecule;  the  latter  approach  is  indicated  for  studies  of  higher 
precision.  When  high  chemical  accuracy  is  required,  as  for  detailed  studies  of  the  ground  or  a 
particular  excited  state  of  a  system,  a  more  extended  basis  must  be  used.  Double-zeta  plus 
polarization  functions  or  optimized  MOs  are  usually  required  as  a  minimum  representation  for 
reliable  calculated  results  of  chemical  accuracy. 

The  chosen  basis  sets  give  good  results  only  when  used  in  a  maximally  flexible  manner.  This 
implies  the  construction  of  perturbation  expansions  or  the  use  of  Cl  wavefunctions  with  all  kinds 
of  possible  orbital  occupancies,  so  that  the  correlation  of  electrons  into  overall  states  can  adjust  to 
an  optimum  form  at  each  geometrical  conformation  and  for  each  state.  Except  when  well-defined 
pairings  exist,  as  for  closed  shell  and  exchange  dominated  systems,  a  single-configuration  study 
(even  of  Hartree-Fock  quality)  will  be  inadequate. 

Proper  electronic  states  for  systems  composed  of  light  atoms  should  possess  definite 
eigenvalues  of  the  spin  operator  as  well  as  an  appropriate  geometrical  symmetry.  The  geometri¬ 
cal  symmetry  can  be  controlled  by  the  assignment  of  orbitals  to  each  configuration,  but  the  spin 
state  must  be  obtained  by  a  constructive  or  projective  technique.  Formulas  have  been  developed 
(Reference  15)  for  projected  construction  of  spin  states  from  orthogonal  orbitals,  and  programs 
implementing  these  formulas  have  been  in  routine  use  at  UTRC  for  several  years.  One  of  the  least 
widely  appreciated  aspects  of  the  spin-projection  problem  is  that  the  same  set  of  occupied  spatial 
orbitals  can  sometimes  be  coupled  to  give  more  than  one  overall  state  of  a  given  spin  quantum 
number.  It  is  necessary  to  include  in  calculations  all  such  spin  couplings,  as  the  optimum  coupling 
will  continuously  change  with  changes  in  the  molecular  conformation.  This  is  especially  important 
in  describing  degenerate  or  near-degenerate  excited  electronic  states. 
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In  the  sections  below,  we  describe  the  several  mathematical  approaches  that  are  applicable 
to  calculation  of  the  electronic  structure  of  molecules,  to  the  calculation  of  potential  energy 
surfaces  for  chemically  reacting  systems,  and  to  the  subsequent  calculation  of  radiation  and 
collisional  processes.  Since  several  different  approaches  are  indicated,  we  describe  their  expected 
areas  of  applicability. 

2.1  Method  of  Ab /n/f/o  Calculation 

2.1.1  Born-Oppenheimer  Separation 

For  a  system  of  n  electrons  and  N  nuclei,  and  considering  only  electrostatic  interactions 
between  the  particles,  we  have  for  the  total  Hamiltonian : 


V  ^ 


N  N  N  n  n  n 

X  +  Xv.v,+  X  j  Vi-Vj 

•»*j 


(1) 


where ; 


Md  =  -^iv?+V''(r„R„)  (2) 

"Wifi 

and  where  Wf.,  aia ,  Mj,  are  the  masses  of  the  electron,  atom,  and  combined  system  mass,  respec¬ 
tively.  Now  since  the  ratios  me/nv, ,  and  nie/Mx  are  both  small,  (2  ■  10^  -  5  •  10"^)  we  can  effect 
a  separation  of  the  electronic  and  nuclear  coordinates  treating  the  total  wavefunction  as  u  product 
of  a  nuclear  and  an  electronic  part.  We  have : 


V<rn,RN)  =  Z^k(RNM(?n>RN) 

k 


(3) 


where  Rn)  is  an  electronic  wavefunction  parametric  in  the  nuclear  coordinates  as  given  in 

Equation  (3)  and  Xk(RN)  are  nuclear  motion  wavefunctions  which  satisfy  (neglecting  terms  of  the 
order  of  me/ma ): 
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The  cross  term  in  Vq  •  can  be  eliminated  by  a  proper  change  of  variables  and  Equation 
(4)  then  reduces  to  a  3N-3  dimensional  Schrodinger  equation. 

For  most  systems,  where  the  velocity  of  motion  of  the  nuclei  is  slow  relative  to  the  electron 
velocity,  this  decoupling  of  electron  and  nuclear  motion  is  valid  and  is  referred  to  as  the  adiabatic 

approximation.  Equation  (3)  thus  defines  an  electronic  eigenstate  V’kOn.  Rn).  parametric  in  the 

nuclear  coordinates,  and  a  corresponding  eigenvalue  EkCRs)  which  is  taken  to  represent  the 
potential  energy  curve  or  surface  corresponding  to  state  k. 

2.1.2  Variational  Methods 

By  an  ab  initio  method  is  meant  one  that  starts  from  a  zero-order  Hamiltonian  which  is 
exact  except  for  relativistic  and  magnetic  effects,  and  which  involves  the  evaluation  of  electronic 
energies  and  other  relevant  quantities  for  wavefunctions  which  are  properly  antisymmetrized  in 
the  coordinates  of  all  the  electrons.  For  a  system  containing  n  electrons  and  M  nuclei,  the 
zero-order  Hamiltonian  depends  parametrically  on  the  nuclear  positions  and  is  of  the  form : 


•f  N  0  M  M  y  ^ 

2itl  itljtlri-Rj  i^^jRi-Rj 


(5) 


where  Z\  and  Ri  are  the  charge  and  position  of  nucleus  i,  Fj  is  the  position  of  electron  j,  and  Vj^ 

is  the  Laplacian  operator  for  electron  j.  All  quantities  are  in  atomic  units,  i.e., lengths  in  bohrs, 
energies  in  hartrees  (1  hartree  =  2  Rydbergs).  The  many-electron  wavefunction  consists  of  one,  or 
a  linear  combination,  V'  =  ^  terms  of  the  form : 


i»i 

where  each  is  a  spatial  orbital,  Ji.  is  the  antisymmetrizing  operator,  0,  is  the  spin-projection 
operator  for  spin-quantum  number  S,  and  is  a  product  of  a  and  one-electron  spin 
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functions  of  magnetic  quantum  number  M^.  No  requirement  is  imposed  as  to  the  double  occu¬ 
pancy  of  the  spatial  orbital,  so  linear  combinations  of  the  form  given  by  Equation  (6)  can  describe 
a  completely  general  wavefunction.  The  spatial  orbitals  may  be  whatever  basis  orbitals  have 

been  introduced,  arbitrary  linear  combinations  thereof,  or  specific  linear  combinations  deter¬ 
mined  pursuant  to  the  particular  calculational  method  in  use. 

The  spatial  orbitals  4>it\  >  spin  functions  d^ni>  ^nd  the  coefficients  of  different  W ,  if  a 
linear  combination  of  %  is  used,  may  be  explicitly  determined  by  invoking  the  variational 

principle.  Various  specific  methods  are  described  below  for  determining  wavefunctions.  However, 
we  should  first  observe  that  the  adequacy  of  an  initio  calculation,  or  for  that  matter  any  energy 
calculation,  will  depend  crucially  upon  the  extent  to  which  the  wavefimction  can  be  qualitatively 
appropriate.  Some  of  the  considerations  surrounding  the  choice  of  a  wavefunction  are  the 
following: 

(1)  Necessity  that  the  wavefunction  possess  sufficient  flexibility  to  be  able  to  describe 
dissociation  to  the  correct  atomic  and  molecular  fragments  as  various  intemuclear  separations  are 
increased . 

(2)  Maintenance  of  equivalent  quality  of  calculation  for  nuclear  geometries  differing  in  the 
nature  or  number  of  chemical  bonds. 

(3)  Ability  to  describe  degenerate  or  near^degenerate  electronic  states  when  they  are 
pertinent . 

(4)  Ability  to  describe  different  electronic  states  to  equivalent  accuracy  when  their  interrela¬ 
tion  (e.g.,  crossing)  is  relevant,  in  particular,  ability  to  describe  ionic-valence  state  mbdng  • 

(5)  Ability  to  represent  changes  in  the  coupling  of  electron  spins  as  bonds  are  broken  or 
reformed. 

The  foregoing  considerations  indicate  that  it  will  often  be  necessary  to  consider  wavefunc¬ 
tions  with  more  than  a  minimum  number  of  singly-occupied  spatial  orbitals,  and  that  there  will 
be  many  potential  curves  or  surfaces  for  which  a  wavefunction  consisting  of  a  single  W  cannot 

suffice.  It  will  then  be  necessary  to  allow  mixing  of  with  different  degrees  of  orbital  spatial 

occupancy  so  as  to  obtain  smooth  transitions  from  the  occupancies  characteristic  of  separated 
atoms  or  molecules  (or  ions)  to  those  characteristic  of  a  comjjound  system  or  a  different  fragmen¬ 
tation. 

Another  implication  of  the  considerations  surrounding  the  choice  of  a  wavefunction  is 
related  to  the  treatment  of  electron  spin.  Not  only  is  it  necessary  to  require  that  the  wavefunction 
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be  an  eigenfunction  of  and  Sz,but  it  is  also  necessary  to  take  account  of  the  fact  that  under 
many  conditions,  there  will  be  more  than  one  spin  eigenfunction  of  given  S  and  ms-  The  different 
spin  eigenfunctions  correspond  to  different  couplings  among  the  individual  spins.  Since  reactive 
processes  involve  the  breaking  and  forming  of  electron-pair  bonds,  they  must  necessarily  be 
accompanied  by  reorganizations  of  the  spin  coupling.  A  failure  to  take  account  of  this  will  lead  to 
qualitatively  inappropriate  wavefunctions. 

In  Hartree-Fock  calculations,  y(R)  is  restricted  to  a  single  ^  which  is  assumed  to  consist 
as  nearly  as  possible  of  doubly-occupied  orbitals.  The  orbitals  <f>^i  are  then  selected  to  be  the 
linear  combinations  of  basis  orbitals  best  satisfying  the  variational  principle.  Writing 
=  2  aviZv.  the  avi  are  determined  by  solving  the  matrix  Hartree-Fock  equation  : 


X  =  Ci  X  Sivavi  (each  A)  (7) 

V  V 


where  (j  is  the  orbital  energy  of  . 

The  Fock  operator,  F^v.has  been  thoroughly  discussed  in  the  literature  (Reference  16)  and  de¬ 
pends  upon  one-  and  two-electron  molecular  integrals  and  upon  the  avi  ■  This  makes  Equa¬ 
tion  (7)  nonlinear  and  it  is,therefore,solved  iteratively.  UTRC  has  developed  programs  for  solving 
Equation  (7)  for  both  closed  and  open-shell  systems,  using  basis  sets  consisting  of  either  Slater- 
type  or  Gaussian-type  atomic  orbitals.  Examples  of  their  use  are  in  the  literature  (Reference  17). 

In  conHguration  interaction  calculations,  the  overall  wavefunction  has  more  than  one  term; 
W  and  the  c^  are  determined  by  invoking  the  variational  principle  to  obtain  the  secular 

equation: 


^(H^v-ws^v)cv  =  0  (each  ft) 

V 

where: 


(8) 


(9) 
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Equation  (8)  is  solved  by  matrix  diagonalization  using  either  a  modified  Givens  method  (Refer¬ 
ence  18)  or  a  method  due  to  Shavitt  (Reference  19)  or  Raffenetti  (Reference  20). 

The  matrix  elements  and  may  be  reduced  by  appropriate  operator  algebra  to  the 
forms: 


H^v  =  Xfp(dM|OsP|0M)  (  fl  %i(ri,RN)|K(RN)P|  fl  'Pvi(ri,RN))  (10) 

P  i=l  i=l 

(^|OsP|0M)(n%i(ri.RN)lPl  n'I'yi(r,RN))  (11) 

P  i  =  l  i-1 

where  P  is  a  permutation  and  €p  its  parity.  The  sum  is  over  all  permutations.  (dn)|0,P|d„)  is  a 

“Sanibel  coefficient”  and  the  remaining  factors  are  spatial  integrals  which  can  be  factored  into 
one-  and  two-electron  integrals.  If  the  <pid  are  orthonormal.  Equations  (10)  and  (U)  become  more 

tractable  and  the  H^y  and  may  be  evaluated  by  explicit  methods  given  in  the  literature  (Ref¬ 
erence  15).  Computer  programs  have  been  developed  for  carrying  out  this  procedure,  and  they 
have  been  used  for  problems  containing  up  to  106  total  electrons,  10  unpaired  electrons,  and  sev¬ 
eral  thousand  configurations. 

The  Cl  studies  described  above  can  be  carried  out  for  any  orthonormal  set  of  for  which 
the  molecular  integrals  can  be  calculated.  Programs  developed  by  UTRC  make  specific  provision 
for  the  choice  of  the  as  Slater-type  atomic  orbitals,  as  Gaussian-type  orbitals,  as  symmetry 

molecular  orbitals,  as  Hartree-Fock  orbitals,  or  as  more  arbitrary  combinations  of  atomic 
orbitals. 

The  one-  and  two-electron  integrals  needed  for  the  above  described  method  of  calculation 
are  evaluated  for  STO  s  by  methods  developed  by  this  Center  (Reference  21).  For  Gaussian 
orbitals,  either  the  Camegie-Mellon  integral  package  (Reference  22)  or  the  integral  routines 
incorporated  in  the  HONIX)  or  GAMESS  programs  (Reference  23)  can  be  employed. 

2.1.3  Configuration  Selection 

Using  a  double-zeta  plus  polarization  basis  set  of  one-electron  functions,  a  typical  system 
can  easily  have  of  the  order  of  10^  configurations  in  full  Cl  (that  resulting  from  all  possible  orbital 
occupancies).  It  is, therefore, essential  to  identify  and  use  the  configurations  describing  the 
significant  part  of  the  wavefiinction.  There  are  several  ways  to  accomplish  this  objective.  First, 
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atomic-orbital  occupancies  may  be  screened  to  eliminate  configurations  with  excessive  formal 
charge.  Alternatively,  in  a  molecular-orbital  framework,  configurations  may  be  eliminated  with 
excessive  numbers  of  antibonding  orbitals.  A  third  possibility  is  to  carry  out  an  initial  screening  of 
configurations,  rejecting  those  whose  diagonal  energies  and  interaction  matrix  elements  do  not 
satisfy  energy  significance  criteria. 

Another  common  method  of  classifying  configurations  is  to  examine  the  total  number  of 
orbitals  in  the  wavefunction  that  differ  from  the  SCF  reference  wavefunction.  We  write  the  Cl 
wavefunction  as : 


V'a  =  Co^o  +  (12) 

'  j 

where  a  single  excitation  function,  <p  f  differs  from  the  SCF  reference  V'o  by  one  orbital  and  the 
double  excitation  function  <p^  by  two,  and  so  on.  The  Cl  coefficients,  Cj,  are  then  determined 

variationally  to  yield  the  lowest  possible  total  energy.  In  the  limit  of  a  complete  basis  set  (N-*(»), 
and  where  all  possible  substitutions  are  included  in  V’a .  die  variational  energy  approaches  the 
correct  nonrelativistic  Bom-Oppenheimer  result.  Errors  arise  from  a  truncation  of  the  functions 
used  to  determine  the  SCF  reference  wavefunction  and  from  the  truncation  of  the  excitation  func¬ 
tions  series  in  V'ci  •  The  reference  wavefunction  V’o  will  typically  be  the  same  for  a  Cl  or  Many 

Body  Perturbation  Theory  (MBPT)  calculation;  however,  even  for  small  triatomic  molecules, 

V;ci  becomes  a  function  of  a  very  large  number  of  basis  functions.  Additionally,  the  reference 

state  itself  may  be  multidimensional  in  order  to  describe  systems  such  as  diradicals  or  systems 
with  orbital  degeneracy.  Because  of  this,  one  must  truncate  the  expansion  and  eliminate  unimpor¬ 
tant  configurations.  In  general  one  can  show  that : 


(tpM>'^)  =  =  0 


(13) 


In  words,  the  matrix  elements  between  the  reference  wavefunction  and  triple  and  quadruple  exci¬ 
tations  is  zero.  Thus  to  first  order,  only  single  and  double  excitations  contribute.  Because  of  this, 
many  Cl  calculations  attempt  to  include  all  single  and  double  excitations  in  the  expressions  for 
.  To  go  beyond  this,  generally  more  than  one  reference  wavefunction  is  used.  Programs  to  han¬ 
dle  configurations  on  all  the  above  criteria  are  available  at  UTRC. 
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Other,  potentially  more  elegant  methods  of  configuration  choice  involve  formal  approaches 
based  on  natural  orbital  (Reference  24)  or  multiconfiguration  SCF  (Reference  25)  concepts.  To 
implement  the  natural-orbital  approach,  an  initial  limited  Cl  wavefunction  is  transformed  to 
natural  orbital  form,  and  the  resulting  natural  orbitals  are  used  to  form  a  new  Cl.  The  desired 
result  is  a  concentration  of  the  bulk  of  the  Cl  wavefunction  into  a  smaller  number  of  significant 
terms.  The  multiconfiguration  SCF  approach  is  more  cumbersome,  but  in  principle  more 
elective.  It  yields  the  optimum  orbital  choice  for  a  preselected  set  of  configurations.  This 
approach  works  well  when  a  small  number  of  dominant  configurations  can  be  readily  identified. 
The  method  is  described  briefly  below. 

2.1.4  Multiconfiguration  -  Self-Consistent  Field  Method  (MC-SCF) 

Tlie  Hartree-Fock  self-consistent  field  method  has  been  proven  to  be  a  powerful  tool  for  the 
calculation  and  understanding  of  many  ground  state  properties  of  molecules  in  the  vicinity  of 
their  equilibrium  structure.  However,  in  most  cases  the  one  configuration  Hartree-Fock  approxi¬ 
mation  is  not  adequate  to  properly  describe  the  dissociation  of  molecular  bonds.  Also,  many 
excited  states  caimot  be  represented  by  a  single  configuration  wave  function.  In  order  to  calculate 
properties  for  such  states,  or  to  investigate  the  formation  of  molecular  bonds,  one  often  needs 
multiconfiguration  wave  functions  for  which  both  the  linear  coefficients,  C,  of  the  configuration 
expansion: 


(14) 


as  well  as  the  set  of  orthonormal  molecular  orbitals  {^j}  ,  from  which  the  configurations  are 
constructed,  are  optimized  according  to  the  variational  principle.  As  is  well  known,  this  “MC- 
SCF”  problem  presents  many  more  difficulties  than  the  simple  one  configuration  Hartree-Fock 
case  and  much  work  has  been  devoted  to  obtaining  convergent  solutions  during  the  last  decade. 

The  difficulties  mainly  arise  from  the  fact  that  for  general  MC-SCF  wave  functions  the 
energy  is  not  invariant  with  respect  to  rotations  between  occupied  orbitals.  Hence,  instead  of  a 
relatively  simple  pseudo-eigenvalue  equation  in  the  one  determinant  case,  the  set  of  coupled  Fock 
equations: 
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with  the  hermiticity  conditions: 


fi)  = 


(16) 


has  to  be  solved,  ejj  are  Lagrange  multipliers  which  account  for  the  orthonormality  constraints 
imposed  on  the  orbitals.  The  Fock  operators  Fij  depend  on  the  orbitals  [^,j  and  the  set  of  Cl  co¬ 
efficients  (Cij .  In  analogy  to  the  one  determinant  case,many  attempts  have  been  made  to  solve 

these  equations  iteratively  by  keeping  the  Fock  operators  fixed  in  each  iteration  step.  Then  the 
Lagrange  multipliers  can  be  expressed  by  coupling  operators  constructed  so  that  the  Fock 
equations  are  transformed  into  pseudo-eigenvalue  equations  yielding  the  improved  orbitals.  These 
are  used  in  a  second  step  to  determine  new  Cl  coefficients  by  diagonalizing  the  Cl  matrix 

■  The  convergence  of  these  algorithms,  however,  has  often  been  found  to  be  poor. 


A  second  group  of  MC-SCF  methods  is  based  on  the  generalized  Brillouin  theorem.  In  these 
methods,  the  orbital  changes  are  derived  from  the  coefficients  of  a  Cl  expansion  consisting  of  the 
MC-SCF  wavefunction  and  all  one-electron  singly  excited  configurations.  A  computer  program 
(ALIS)  implementing  this  method  has  been  developed  by  Ruedenberg,  et  al.  (Reference  26).  A 
somewhat  more  elegant  program  (GAMESS)  which  also  incorporates  analytical  gradients  is  also 
available  (Reference  23).  More  recently,  various  methods  have  been  proposed  which  are  based  on 
direct  minimization  of  the  energy,  avoiding  the  Fock  operators  altogether.  Such  methods  are  now 
being  studied  in  our  laboratory  and  will  be  incorporated  into  our  existing  computer  programs  if 
they  prove  to  be  highly  efficient. 

2.1.5  Many-Body  Perturbation  Theory  (MBPT) 

In  MBPT  we  again  begin  with  the  SCF  wavefunction  (%)  as  our  reference  and  attempt  to 

account  for  Ec,  the  correlation  energy.  The  concept  of  excitation  functions  described  in  the  above 
section  on  Cl  calculations  carries  over  to  MBPT  calculations.  In  MBPT  one  can  write  the 
wavefunction,  ^p,as: 


V'p  =  e'T^o) 


(17) 
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where  T  is  an  excitation  operator  defined  as : 


T  =  Ts  +  Td  +  Tx  +  Tq...  (18) 

where  S,  D,  T,  and  Q  refer  to  single,  double,  triple,  and  quadruple  substitutions.respectively.  One 
can  write  Tn,  in  general,  where  n  refers  to  the  number  of  excitations  as : 

"•  ijk... 
abc... 

where  a,b,c,...  are  excited  orbitals  and  i,j,k,...  are  orbitals  occupied  in  •  The  total  energy  is  now 
given  by : 


Embpt  =  (^olKe^|0o)  (20) 

Tb  evaluate  Embpt,  the  ‘  from  above  must  be  determined.  An  equivalent  expression  that 
makes  the  perturbation  expansion  clearer  is: 


Embpt  =  X  (^ol36{(Eo-Ko)”^K]^|^o)  (21) 

k»0 

where  the  sum  is  over  only  so-called  linked  diagrams,  Dto  has  eigenfunctions  0©  >  and  the  expan¬ 
sion  is  of  orders  in  the  perturbation,  V  =  (36-36o)  •  The  k  =  0  term  gives  the  reference  energy 
and  for  k  >  0  correlation  corrections  are  included.  In  practice  the  MBPT  total  energy  is  calcu¬ 
lated  by  truncating  the  T  operator  expansion  and  projecting  Ke^*|0o)  onto  the  appropriate  n- 

space.  This  leads  to  a  set  of  nonlinear  coupled  equations  for  the  t*  ’  coefficients  which 

correspond  to  the  Cl  expansion  coefficients.  The  equations  are  solved  iteratively  and  Embpt 
evaluated.  In  practice, T4  is  an  upper  limit  that  corresponds  to  quadruple  substitutions.  The  series 
is  an  oscillatory  convergent  sum,  which  in  practice  has  proven  to  be  at  least  as  accurate  as  Eq 
with  single  and  double  excitations  included. 

The  best  possible  method  to  use  would  be  a  hill  Cl  (all  possible  configurations)  with  a 
complete  basis  set.  However,  since  the  number  of  configurations  is  proportional  to  (nf ,  where  n 
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is  the  number  of  basis  functions  and  I  is  the  level  of  excitation,  it  would  be  prohibitive  to  even 
include  all  single  and  double  excitations.  This  truncation  causes  the  loss  of  size  consistency, 
which  implies  that  the  energy  calculated  for  A  and  B  as  a  molecular  system,  but  dimensionally  far 
apart,  is  the  sum  of  the  energy  calculated  for  A  and  B  separately.  In  a  size  extensive  calculation  the 
energy  is  proportional  to  the  size  of  the  system.  These  properties  are  very  important  if  one  wishes 
to  compute  correct  relative  energies  on  a  potential  energy  surface,  a  necessary  criteria  for  defining 
the  reactive  pathways  of  interest  in  this  research  program.  MBPT,  on  the  other  hand,  is  guaranteed 
to  have  the  correct  size-dependence  because  the  expansions  contain  only  the  so-called  linked 
diagrams.  In  addition,  because  of  the  computational  efficiency  of  MBPT,  calculations  can  be 
performed  up  to  fourth  order  in  the  perturbation  expansion  and  include  single,  double,  triple,  and 
quadruple  excitations  in  the  calculation  of  the  correlation  energy.  Such  calculations  are  usually 
performed  with  at  least  a  split  valence  plus  polarization  basis  set.  For  the  light  element  compounds 
which  exhibit  ionic  bonding,  the  inclusion  of  diffuse  basis  functions  and  possibly  higher 
polarization  d  or  f-functions  may  be  required  for  a  more  quantitative  treatment. 

2.2  Ab  Initio  Gaussian  Wavefunction  Electronic  Structure  Codes 

Owing  to  the  complexity  of  evaluating  multicenter  electron  repulsion  integrals  over  Slater- 
type  exponential  orbitals,  various  groups  have  adopted  a  computational  approach  to  electronic 
structure  calculations  based  on  gaussian  orbitals.  A  highly  developed  series  of  programs,  named 
GAUSSIAN  92,  is  available  from  Camegie-Mellon  University  (Reference  22).  A  second  series  of 
programs  has  evolved  from  the  original  version  of  Dupuis  and  King's  (Reference  27)  HONDO 
code.  This  code  has  been  further  developed  by  separate  groups  as  GAMESS  (Reference  23), 
CADPAC  (Reference  28)  and  HONDO  7  (Reference  29).  Finally,  a  new  code  named 
COLUMBUS  (Reference  30  and  31)  has  been  developed  by  the  Ohio  State/  Argonne/Batelle 
group.  UTRC  has  been  acting  as  a  beta  test  site  for  this  new  development.  A  brief  description  of 
the  features  of  these  codes  follows. 

2.2.1  GAUSSIAN  92 

GAUSSIAN  92  is  a  connected  system  of  programs  for  performing  ab  initio  molecular 
orbital  (MO)  calculations.  It  represents  further  development  of  the  GAUSSIAN 
70/76/80/82/86/88/90  systems  already  published.  The  contributors  to  this  program  include;  M.  J. 
Frisch,  G.  W.  Trucks,  M.  Head-Gordon,  P.  M.  W.  Gill,  M.  W.  Wong,  J.  B.  Foresman,  B.  G. 
Johnson,  H.  B.  Schlegel,  M.  A.  Robb,  E.  S.  Replogle,  R.  Gomperts,  J.  L.  Andres,  K. 
Raghavachari,  J.  S.  Binkley,  C.  Gonzalez,  R.  L.  Martin,  D.  J.  Fox,  D.  J.  Defrees,  J.  Baker,  J.  J. 
P.  Stewart,  and  J.  A.  Pople.  GAUSSIAN  92  was  originally  implemented  on  the  chemistry 
department  DEC  VAX  1 1/780  computer  at  Camegie-Mellon  University.  Since  then,  this  program 
has  been  installed  on  a  number  of  different  computers. 

GAUSSIAN  92  was  designed  with  a  transparent  input  data  stream,  making  this  program 
very  user  friendly.  All  of  the  standard  input  is  free-format  and  mnemonic.  Reasonable  defaults  for 
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input  data  were  provided,  and  the  output  is  intended  to  be  self-explanatory.  Mechanisms  are 
available  for  the  sophisticated  user  to  override  defaults  or  interface  their  own  code  to  the 
GAUSSIAN  system.  In  this  respect,  we  used  GAUSSIAN  92  as  a  fundamental  framework  for 
.several  applications.  Options  were  incorporated  into  this  code  to  provide  capabilities  beyond 
Hartree-Fock  and  various  perturbation  theory  options. 

The  capabilities  of  the  GAUSSIAN  92  system  include: 

a)  Calculation  of  one-  and  two-electron  integrals  over  s,  p,  d,  and  f  contracted  gaussian  functions. 
The  basis  functions  can  either  be  cartesian  gaussians  or  pure  angular  momentum  functions  and 
a  variety  of  basis  sets  are  stored  in  the  program  and  can  be  requested  by  name. 

b)  Self-consistent  field  calculations  for  restricted  closed-shell  (RHF),  unrestricted  open-shell 
(UHF),  and  open-shell  restricted  (ROHF)  Hartree-Fock  wavefunctions  as  well  as  those  types 
of  multiconfigurational  wavefunctions  that  fall  within  the  Generalized  Valence  Bond-Perfect 
Pairing  (GVB-PP)  formalism. 

c)  Evaluation  of  various  one-electron  properties  of  the  Hartree-Fock  wavefunction,  including 
Mulliken  population  analysis,  multipole  moments,  and  electrostatic  fields. 

d)  Automated  geometry  optimization  to  either  minima  or  saddle  points,  and  analytical  or  numerical 
differentiation  to  produce  force  constants,  polarizabilities,  and  dipole  derivatives.  This  feature 
can  be  used  to  develop  minimum  energy  reaction  paths  along  a  complicated  many-dimensional 
potential  energy  surface. 

e)  Correlation  energy  calculations  using  Mdller-Plesset  perturbation  theory  carried  to  second, 
third,  or  fourth  order. 

f)  Correlation  energy  calculations  using  configuration  interaction  (Cl)  with  either  all  double 
excitations  (CID)  or  all  single  and  double  excitations  (CISD). 

g)  Correlation  energy  calculations  using  coupled  cluster  theory  with  double  substitutions  (CCD). 

h)  Correlation  energy  calculations  using  quadratic  convergence  SCF  (QCSCF).  This  is  a  new 
highly-efficient  size-consistent  method  recently  developed  by  Pople. 

i)  Analytic  computation  of  the  nuclear  coordinate  gradient  of  the  RHF,  UHF,  ROHF,  GVB-PP, 
MP2,  CID,  and  RCISD  energies. 

j)  Computation  of  force  constants  (nuclear  coordinant  second  derivatives),  polarizabilities, 
hyperpolarizabilities,  dipole  derviatives,  and  polarizability  derivatives  either  analytically  or 
numerically. 
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k)  Harmonic  vibrational  analysis. 

l)  Determination  of  intensities  for  vibrational  transitions  at  the  HF,  MP2,  Cl  levels. 

m)  Testing  the  SCF  wavefunctions  for  stability  under  release  of  constraints. 

n)  Correlated  electron  densities  and  properties. 

o)  Minimum-energy  pathway  following  from  products  to  reactants  through  a  transition 
state- intrinsic  reaction  coordinate  finder. 

2.2.2  GAMESS 

A  wide  range  of  quantum  chemical  computations  is  possible  using  GAMESS  (Reference  23),  a 
refinement  by  M.  Schmidt  and  S.  Elbert  of  the  original  HONDO  code  (Reference  27). 

The  capabilities  of  this  code  include: 

a)  Calculations  of  RHF/UHF/ROHF/GVB-SCF  molecular  wavefunctions. 

b)  Calculations  of  multiconfiguration  SCF  (MCSCF)  wavefunctions. 

c)  Calculations  of  Cl  wavefunctions  using  the  unitary  group  method. 

d)  Optimization  of  molecular  geometries  using  an  energy  gradient  in  terms  of  flartesian  or 
internal  coordinates. 

e)  Searches  for  potential  energy  surface  saddle  points. 

f)  Tracing  the  intrinsic  reaction  path  from  a  saddle  point  to  reactants  or  products. 

g)  Computation  of  normal  modes  and  vibrational  frequencies. 

h)  Calculation  of  the  following  properties: 

1 .  Dipole,  quadrupole,  and  octupole  moments. 

2.  Electrostatic  potentials 

3 .  Electric  field  and  electric  field  gradients 

4.  Electron  density  and  spin  density 

5 .  Mulliken  and  Ldwdin  population  analysis 

6.  Localized  orbitals  by  the  Boys  method 

7 .  Virial  theorem  and  energy  components. 
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GAMESS  is  a  synthesis  with  many  major  modifications  of  several  programs.  A  hirge  part  of  the 
program  is  from  HONDO.  For  pure  sp  basis  sets,  the  HONDO  symmetry  and  supermatrix 
procedure  were  adapted  from  GAUSSIAN  76  and  GRADSCF  integrals,  both  for  the  SCF  and 
gradient  parts.  The  GVB  section  is  a  heavily  modified  version  of  GVBONE.  A  Boys  localization 
algorithm  is  implemented  from  a  heavily  modified  version  of  Streitweiser's  QCPE  program. 

The  Cl  module  is  based  on  Brooks  and  Schaefer's  unitary  group  program  which  was 
modified  to  run  within  GAMESS,  using  a  Davidson  eigenvector  method  written  by  S.  T.  Elben. 
The  MCSCF  module  is  a  Newton-Raphson  procedure,  developed  at  NRCC,  based  on  the  unitary 
group  Cl  package.  The  intrinsic  reaction  coordinate  pathfinder  was  written  at  North  Dakota  State 
University. 

2.2.3  CADPAC 

The  Cambridge  Analytic  Derivatives  Package  (CADPAC)  (Reference  28)  is  a  group  of 
programs  which  has  been  under  development  at  Cambridge  University,  UK.  It  originated  as  a 
version  of  Dupuis  and  King's  HONDO  program.  From  its  initial  state  as  an  SCF  gradient  package 
for  closed-shell  and  UHF  wavefunctions,  this  program  was  extensively  modified  with  many  of  the 
old  features  being  enhanced  and  many  new  features  being  added.  The  input  data  is  now  in  free 
format  with  a  "keyword"  system,  to  make  the  program  easier  to  use.  The  integral  routines  use 
essentially  the  same  methods  as  those  in  HONDO,  being  based  on  the  Rys  polynomial  method,  but 
were  extended  to  cover  f-functions.  These  routines  were  vectorized  to  take  advantage  of  the  much 
greater  availability  of  supercomputers  such  as  the  Cray  in  recent  years.  The  SCF  programs  are  by 
now  a  blend  of  techniques,  but  still  contain  a  few  features  from  the  initial  HONDO  program, 
particularly  in  the  way  symmetry  is  handled.  There  were  many  modifications  to  improve  efficiency 
and  the  addition  of  level-shifting  and  damping  techniques  and  the  implementation  of  the  DIIS 
method  to  aid  convergence.  A  high-spin  open-shell  SCF  program,  and  a  completely  general  open- 
shell  SCF  program  were  included. 

The  gradient  routines  were  considerably  altered  from  those  in  the  original  HONDO,  the 
method  was  changed  to  one  which  is  more  efficient  and  easier  to  extend  to  higher  order 
derivatives.  These  routines  also  work  for  f-functions.  The  original  program  capabilities  for  the 
calculation  of  force-constants  by  numerical  differentiation  of  gradients  and  for  the  optimization  of 
geometries  remain  essentially  intact.  However,  there  is  now  a  choice  of  two  optimization 
algorithms  with  the  inclusion  of  Schlegel's  method. 

In  addition  to  the  above  modifications,  a  range  of  extra  facilities  were  developed.  These 
included  a  more  powerful  method  of  calculating  one-electron  properties,  analyzing  molecular 
charge  distributions,  and  various  "post-SCF"  stages,  beginning  with  a  4-index  transformation. 
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These  new  sections  include  M<J)ller-PIesset  perturbation  theory  for  total  energies  and  molecular 
properties,  and  coupled  Hartree-Fock  calculations  of  molecular  polarizabilities.  The  polarizability 
routines  can  calculate  dynamic  properties  at  real  or  imaginary  frequencies,  and  obtain  dispersion 
coefficients.  It  is  also  possible  to  use  CHF  theory  to  obtain  the  perturbations  due  to  nuclear 
displacements.  These  can  be  used  to  obtain,  analytically,  all  the  dipole  and  quadrupole  moment 
derivatives  of  a  molecule.  They  also  form  part  of  the  most  important  addition  to  the  package,  the 
section  which  calculates  analytic  second  derivatives  of  the  energy.  This  is  a  powerful  technique 
whose  speed  and  accuracy  represents  a  considerable  improvement  over  numerical  differentiation. 

The  capabilities  of  the  latest  CADPAC  code  (Version  4)  include 

a)  The  evaluation  of  one-  and  two-electron  integrals  over  contracted  cartesian  gaussian 

basis  functions  of  type  s,  p,  d.  or  f. 

b)  SCF  calculations  for  closed  shell,  open-shell,  UHF  and  generalized  open-shell  tech¬ 
niques. 

c)  Calculation  of  one-electron  properties  for  these  types  of  wavefunction,  including  a  dis¬ 
tributed  multipole  analysis. 

d)  Calculation  of  the  gradients  of  the  SCF  energy. 

e)  Use  of  the  gradients  for  automatic  geometry  optimization,  and  for  the  calculation  of 
force  constants  by  numerical  differentiation.  There  is  a  choice  of  two  optimization  al¬ 
gorithms. 

0  Transformation  of  the  integrals  from  the  atomic  orbital  to  the  molecular  orbital  basis. 

g)  M<})ller-Plesset  perturbation  theory  calculations  to  third  order  in  the  energy  and  sec¬ 
ond  order  in  the  one-electron  properties. 

h)  Coupled  Hartree-Fock  calculations  of  polarizabilities,  including  frequency  dependence, 
and  magnetizabilities. 

i)  Coupled  Hartree-Fock  calculations  of  the  perturbation  due  to  nuclear  displacements. 

j)  Calculation  of  the  dipole  and  quadrup>ole  moment  derivatives. 

k)  Calculation  of  the  second  derivatives  (force  constants)  of  the  energy  by  analytic  meth¬ 
ods. 

l)  Analytic  calculations  of  polarizability  derivatives. 

m)  Calculation  of  infrared  and  Raman  intensities,  and  the  study  of  vibrational  circular 
dichroism. 
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n)  Calculation  of  MP2  gradients,  dipole  moment  derivatives,  polarizabilitites  and  force 

constants  using  analytic  algorithms. 

o)  Spin-projected  UHF  MP2  energies 

2.2.4  HONDO 

The  HONDO  program,  (Reference  27)  originally  developed  by  Dupuis  and  King  at  NRCC 
has  recently  been  refined  and  updated  by  M.  Dupuis  at  IBM-Kingston  (Reference  29). 

The  following  features  are  available  in  the  present  version  (HONDO  7)  of  the  program; 

a)  Single  configuration  self-consistent-field  wavefunctions  (closed  shell  RHF,  spin  unre¬ 
stricted  UHF,  restricted  open  shell  ROHF),  generalized  valence  bond  GVB  and  gen¬ 
eral  multiconfiguration  self-consistent-field  MCSCF  wavefunctions,  and  configuration 
interaction  Cl  wavefunctions  can  be  calculated. 

b)  The  electron  correlation  correction  to  the  energy  of  closed  shell  RHF  wavefunctions 
can  be  calculated  by  means  of  M4>ller-Plesset  (MP)  perturbation  theory  applied  to 
second-,  third-,  and  fourth-order  (with  or  without  the  effects  of  triple  excitations). 

c)  The  effective  core  potential  approximation  can  be  used. 

d)  Optimization  of  molecular  geometries  using  the  gradient  of  the  energy  with  respect  to 
nuclear  coordinates  is  possible  with  all  but  the  Cl  and  MP  wavefunctions.  Optimiza¬ 
tion  can  be  carried  out  in  the  cartesian  space  or  in  the  internal  coordinate  space  with 
the  possibility  of  freezing  some  cartesian  or  internal  coordinates. 

e)  The  force  constant  matrix  in  the  cartesian  space,  and  the  vibrational  spectrum,  includ¬ 
ing  infrared  and  Raman  intensities,  can  be  calculated  with  all  but  the  Cl  and  MP 
wavefunctions. 

f)  Calculation  of  the  dipole  moment  and  polarizability  derivatives  with  respect  to  the  nu¬ 
clear  coordinates  is  possible,  for  use  with  a  previously  calculated  force  constant  ma¬ 
trix.  The  force  constant  matrix  can  be  transformed  to  the  internal  coordinate  basis. 

g)  Itansition  state  structures  can  be  determined  with  all  but  the  Cl  and  MP  wavefunc¬ 
tions  by  talcing  advantage  of  the  energy  gradients. 

h)  The  “Intrinsic  Reaction  Coordinate’'(IRC)  pathway  can  be  determined  with  all  but  the 
Cl  and  MP  wavefunctions  by  taking  advantage  of  the  energy  gradients. 

i)  Molecular  energies  for  several  points  on  a  potential  energy  surface  can  be  calculated  in 
a  single  run. 
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j)  Nongradient  optimization  of  basis  function  exponents  is  possible.  The  source  code 
can  be  modified  to  cany  out  optimization  of  other  nonlinear  parameters,  for  example, 
contraction  coefficients  and  geometrical  parametei^. 

k)  The  following  electronic  properties  can  be  extracted  from  the  waveftinction: 

1.  Dipole  moment 

2.  Quadrupole  moment 

3.  Mulliken  population,  bond  order  and  valency  analyses 

4.  Spin  density  maps 

5.  Electron  density  maps 

6.  Electrostatic  potential  maps 

7.  Localized  orbitals  via  Boys’  method 

8.  Static  dipole  polarizability 

9.  Static  first  and  second  hyperpolarizabilities. 

l)  The  potential  due  to  finite  point  charges  for  a  classical  representation  of  an  environ¬ 
ment,  or  a  uniform  electric  field  can  be  incorporated  into  the  one-electron  Hamil¬ 
tonian. 

2.2.5  COLUMBUS 

The  COLUMBUS  code  (References  30  and  31)  is  a  continuing  development  based  on  a  joint 
project  at  Ohio  State  University  (I.  Shavitt),  University  of  Karlsruhe  (R.  Ahlrichs),  and  Argonne 
National  Laboratory  (R.  Shepard).  The  unique  features  of  this  code  are  the  incorporation  of 
relativistic  or  nonrelativistic  core  potentials  to  permit  analysis  of  heavy  atom  molecular  systems 
and  the  inclusion  of  a  sophisticated  Cl  package  based  on  the  unitary  group  approach  of  Shavitt 
(Reference  32).  The  basic  programs  included  in  the  COLUMBUS  code  are  as  follows; 

1.  AO  Integrals 

This  is  R.  Pitzer’s  version  of  the  integral  package  from  HONDO  with  partial  vectorization 
of  auxiliary  integral  routines.  This  package  can  handle  up  to  g-functions. 

2.  Integral  Tfansformation 

The  transformation  algorithm  is  written  over  all  the  orbitals  which  constitute  a  shell,  using 
contraction  coefficients  defined  over  the  primitive  basis  functions.  Limited  vectorization  is 
possible  over  the  innermost  loops  for  dyadic  operations. 
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3.  SCF  Gradients 

These  routines  are  based  on  the  HONDO  version  but  include  checks  on  integral  symmetry 
to  avoid  operations  over  zero  or  near  zero  elements.  At  this  time,  only  a  first  derivative 
analysis  has  been  included. 

4.  SCF  Energy 

The  SCF  routines  are  typical  of  those  found  in  the  molecular  structure  codes.  They  incor¬ 
porate  level  shifting,  damping,and  the  incorporation  of  Pulay’s  DUS  convergence  accelera¬ 
tion  procedure. 

5.  MCSCF  Analysis 

The  MCSCF  package  is  due  to  R.  Shepard  and  incorporates  extensive  vectorization.  The 
output  vectors  can  be  taken  directly  ft-om  this  package  and  transferred  to  the  Cl  program. 

6.  Multireference  Direct  Cl 

This  set  of  routines  is  based  on  the  graphical  representation  of  the  unitary  group  ap¬ 
proach  (GUGA)  for  constructing  a  Cl  wavefunction.  This  code  differs  from  conventional 
Cl  in  that  a  matrix  representation  of  the  hamiltonian  is  not  explicitly  computed.  The  uni¬ 
tary  group  approach  can  be  much  more  efficient  in  most  applications  since  ~80  percent  of  the 
hamiltonian  matrix  elements  are  usually  zero  but  must  still  be  included.  This  leads  to  the 
storage  and  diagonalization  problems  associated  with  large,  sparse  matrices.  However, 
little  work  was  done  on  analysis  of  excited  states  with  the  same  molecular  symmetry  as 
the  ground  state.  An  approach  that  looks  promising  is  to  define  an  approximate  vector 
corresponding  to  the  excited  state  wavefunction  and  to  iteratively  improve  upon  this  solu¬ 
tion  using  the  direct  Cl  contraction.  The  uncertainties  in  the  procedure  are  applications  to 
situations  where  degeneracies  or  near-degeneracies  arise.  Further  studies  of  this  case  are 
in  progress. 

2.3  Spin-Projected  Unrestricted  Hartree-Fock  Method 

The  unrestricted  Hartree-Fock  (UHF)  method  developed  by  Pople  and  Nesbet  (Refer¬ 
ences  33  and  34)  yields  the  best  single-determinant  approximation  to  the  exact  wavefunction  for 
an  atomic  or  molecular  system.  Such  a  wavefunction  incorporates  correlation  by  allowing  orbitals 
of  different  spin  to  adjust  to  spatially  different  forms,  thus  breaking  the  symmetry  restrictions  of 
the  conventional  (RHF)  method  (References  35  and  36).  It  is  necessary,  however,  to  project  from 
such  a  wavefunction,  a  properly  antisymmetrized  spin  and  angular  momentum  state  in  order  to 
define  eigenstates  and  eigenenergies  corresponding  to  observable  spectroscopic  states. 

Let  |a),  refer  to  the  doubly-occupied  molecular  orbitals  (MOs),  <Pa,  <t>p .  Let  |y),  |d),  \x) 
refer  to  the  singly  occupied  MOs.  We  assume  that  all  MOs  have  been  subjected  to  a 


transformation  to  orthogonal  form.  Let  Cyj  be  the  coefficient  associated  with  the  permutation 
y  - — ►  <5  (this  permutation  is  (-1)  times  the  overlap  of  the  permuted  spin  eigenfunction  with  the 
original  spin  eigenfunction),  and  adopt  the  convention  Cyy  =  -1.  Then  the  expectation  value  of 
the  nonrelativistic  Hamiltonian : 


X^(0+SV(i.j) 

‘  ><j 


(22) 


is  given  by : 


(K)  =  25;(a|U|a)  +  £{)iU|y)  +  2£(a^) 
a  Y  ofi 


afi  ay  ay 

yi  yi 


where  i,  j  refer  to  electron  numbers,  a,  fi..  y,  6  etc.,  refer  to  MOs  and  the  sums  run  over  all 
a,  /5.... 

The  equations  for  determining  the  optimum  MOs  can  be  derived  from  the  following 
variational  form: 


<5  (3fe)  -  2  21  (a  1^)  -  (yW))  '  Z  ® 

afi  afi  yfi  ~ 


(24) 
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These  equations  are  : 


(U  +  2Jc-Kc  +  Jo-^Ko)|a)  = 


(25) 


(U  +  2Jc-Kc  +  Jo)|y)  +  =  X%l^)  +  )  (26) 

6  fi  6 


where: 


l(  aiviv).!)  =  I  dr«;(f)v(r.f>MiV«(f) 
(w|v|W')  =  |drw*(i)v(f,P)w’(r) 


(27) 


p 


h  - 

6 


Kc|W  =  Xl(^|V|W)/5)  KolW  =  ) 

P  6 


Left-multiplying  Equation  (25)  by  |  y]  and  Equation  (26)  by  ( a)  ,  these  equations  can  be  com¬ 
bined  to  yield : 


f)-  =  X(2C,^  +  l)(y«5|V|3a) 


(28) 


Noting  that: 


2;CyX<5r|V|r>-)  =  XC».(<)r|Vlry>, 


(29) 
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Equations  (2S)  and  (26)  can  be  recast  as  follows.  From  Equation  (25) : 


with: 


(30) 


M  =  2(2Cfc+l)|«  )(t(«|V|t)| 

dr 


From  Equation  (26) : 


Z%l/*)  =■  2(»V+  l)l/>)(/?r|V|4)'}  -  OcM’ly} 

/»  /h 


(31) 


with; 


fi  6 


(32) 


Using  Equations  (31)  and  (32)  in  (25)  and  (26),  we  End; 


(U  +  2Jc-Kc  +  Jo-|Ko-|M)la)  =  J^€pa\fi) 


(33) 


(U  +  2Jc-Kc-*-Jo-^Ko+iM^-OcM^))y)=  ) 


(34) 
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Finally,  Equations  (33)  and  (34)  can  be  rearranged  to  yield  the  same  Hermitian  operator  on  the 
left-hand  side  with  the  resulting  definition  of  the  one-electron  eigenenergies  for  the  closed  and 
open-shell  eigenstates: 


6a'  ={alU  +  ZIc-Kc  +  Jo-^Kola  ) 


(35) 


6/  =(ylU  +  2Jc-Ko  +  Jo-|Ko  +  ) 


(36) 


Equations  (35)  and  (36)  represent  the  correct  one-electron  energy  expressions  for  the  spin- 
projected  eigenstate. 

This  formalism  was  developed  to  extend  the  conventional  HF  method  to  include 
split-shell  correlation  and  proper  spin  and  symmetry  projections.  It  is  being  incorporated  into  a 
computer  program  using  gaussian-type  orbitals  (GTOs  )  as  the  elementary  basis  functions.  A 
crucial  feature  of  this  method  is  that  dissociation  always  follows  the  lowest  energy  pathway 
thereby  permitting  a  proper  description  of  bond  formation  and  breakage.  In  contrast,  RHF  or 
MC-SCF  methods  often  exhibit  improper  dissociation  character  or  exhibit  size  inconsistencies 
owing  to  conelation  energy  changes  in  going  from  molecular  geometries  to  separated  atom- 
molecule  or  atom-atom  dimensions. 

2.4  Dissociative  Recombination  Calculations 

The  theory  of  the  capture  of  an  electron  by  a  positive  molecular  ion  was  developed  for 
both  direct  attachment  processes  (Reference  37),  and  for  several  possible  indirect  processes  which 
involve  the  formation  of  an  intermediate  excited  Rydberg  state  belonging  to  the  molecular  ion 
core  configuration  (References  38-40).  Let  V'f  (L*  represent  the  continuum  wavefunction 
describing  the  free  electron  with  energy  6  plus  the  n-1  electron  bound  system  of  the  molecular 
ion.  Let  V'd  (L'  represent  the  wavefunction  for  an  eigenstate  of  the  neutral  molecule  which 
can  be  written,  in  the  Bom-Oppenheimer  approximation,  as  the  product  of  an  electronic  and 
nuclear  wavefunction  in  the  form: 


^d(r,R)  =  v^'(r,R)ed(R) 


(37) 
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The  cross  section  for  electron  capture  in  dissociative-recombination  is  determined  by  the 
asymptotic  form  of  the  nuclear  wavefunction  ^d(R).  Let  the  incident  flux  of  electrons  per  unit 
area  be ; 


2jiH 


6€ 


(38) 


where  <5^  is  a  measure  of  the  energy  spread  of  the  incoming  beam  and  k  is  the  wave  number 
associated  with  the  incoming  electron.  The  outgoing  flux  of  atoms  is  given  by : 


/i 


l^d(R)l' 


(39) 


where  —  is  the  relative  velocity  of  the  separating  neutral  atoms.  Equations  (37)  and  (39)  yield 
for  the  cross  section : 


o(€,Q)  =  lim 


_ ,2 


R— 


ied(R)|- 


(40) 


Equation  (40)  must  Anally  be  averaged  over  all  rotational  orientations  as : 


|o(€.! 


o(€)  =  ;^  (  o(€.Q)dQ 


(41) 


Assuming  that  the  total  cross  section  can  be  written  as  the  product  of  a  resonant  capture  cross 
section  and  a  survival  factor,  in  the  form ; 


0(f)  =  Ocap(f)S(€) 


(42) 
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we  have: 


(43) 


where  the  factor  0}o/(o+  is  the  ratio  of  the  electronic  degeneracies  of  the  neutral  and  ionic  states, 
respectively.  The  capture  width,  Tc  (R,  6) ,  is  determined  by  integrating  over  the  coordinates  of 
the  electronic  wavefunctions: 


re(R,6)  =  Tx  (v^i{r,R)  f  (e.r.R)l36ei-€|V^.R))"  (44) 

Thus,  the  capture  cross  section  can  be  written  as : 

or  equivalently : 

0(e)  =  (^+(R)v^i(R)V(KF.R,€)|V(R,e)|V'd(r,R)?o(R)  (46) 

zm€  (0+  '  ' 

The  nuclear  wavefunction,  |o>  is  energy  normalized  so  that  asymptotically: 


(|rx^)’^^sin(kr  +  ?;) 


(47) 
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Equation  (46)  can  be  cast  into  the  computational  fonn : 


o(6,V)craVsec  -  ^  (r»(R)l 


(48) 


where  the  electron  energy  €  and  the  capture  width  Fc  are  in  atomic  units. 

Assuming  a  Maxwellian  temperature  distribution  for  the  electrons,  the  rate  coefficient  can 
be  written  as : 


or  equivalently : 


2J2  f  * 


(49) 


(50) 


The  capture  width  Fc  can  be  calculated  by  examining  the  high  members  of  the  Rydberg  series  of 
neutral  states  which  have  the  structure  of  a  ground  state  molecular  ion  coupled  to  an  electron  in 
diffuse  hydrogenic  orbitals  of  large  effective  principal  quantum  numbers.  A  direct  calculation  of 
the  capture  width  involves  the  knowledge  of  the  continuum  wavefunction  for  the  electron  as  a 
function  of  the  interparticle  coordinates.  This  approach  is  computationally  very  difficult  and  at 
the  present  stage  of  development  would  probably  lead  to  errors  of  at  least  the  same  magnitude  as 
an  extrapolation  procedure  of  the  conesponding  neutral  Rydberg  states. 

2.4.1  Indirect  Models 

An  indirect  model  for  calculating  electron-molecular  ion  recombination  coefficients  can  also 
be  formulated.  Here  we  postulate  the  formation  of  a  collision  complex  (AB*  for  purposes  of 
illustrating  this  for  a  diatomic  system)  during  the  collision  and  its  subsequent  decay  to  form  A*, 
B**,  or  AB"*"  +  e  (References  38-40).  The  approaches  differ  by  the  nature  of  the  (AB*)  complex. 
The  high-energy  approximation  assumes  that  all  of  the  energy  of  the  incident  electron  goes  into 
exciting  one  of  the  electrons  in  the  AB"*"  core,  resulting  in  the  capture  of  the  incident  electron  into 
a  doubly  excited  state  of  AB.  This  mechanism  is  unlikely  to  be  applicable  to  AB"*"  recombination 
at  thermal  energies  because  of  the  energy  required  to  promote  an  electron  out  of  the  AB'*'  cor.,.  A 
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second  approach  is  more  applicable  to  thermal  energy  collisions,  since  it  requires  that  the 
electron-molecular  ion  collision  results  in  a  vibrationally  excited  molecular  ion  which  captures  an 
electron  in  the  Coulomb  field  of  AB"*".  This  collision  complex  is  equivalent  to  a  vibrationally  or 
rotationaliy  excited  Rydberg  state  of  AB.  Experimental  evidence  for  the  existence  of  such  states 
and  their  autoionization  was  independently  obtained  from  photoionization  spectra  of  AB, 
within  a  few  electron  volts  of  the  threshold  (Reference  41).  Wth  this  model,  the  prediction  of  the 
atomic  products  of  the  recombination  reaction  reduces  to  calculating  the  Rydberg  states  of  AB 
that  are  most  likely  to  be  populated  during  the  initial  recombination,  and  also  those  atomic  states 
to  which  they  subsequently  decay. 

As  contrasted  with  the  direct  recombination  process  which  is  governed  by  the  configuration 
interaction  strength,  the  nonadiabatic  coupling  of  the  electron  and  nuclear  motion  is  the  dominant 
mechanism  in  both  stages  of  the  indirect  recombination  process.  In  the  first  stage,  it  is  the 
vibrational  (or  rotational)  excitation  of  the  AB"*"  core  which  results  in  the  initial  capture  of  the 
incident  electron.  The  subsequent  decay  of  the  collision  complex  is  due  to  a  transition  to  a 
repulsive  molecular  state  which  crosses  the  Rydberg  state.  It  is  then  the  nonadiabatic  coupling  of 
the  electronic  and  nuclear  motion  which  results  in  the  transition  between  these  two  molecular 
states.  Actually,  it  is  the  same  term  in  the  Hamiltonian  for  the  entire  molecular  system  that  causes 
the  key  transition  in  both  stages  of  the  recombination  process.  This  becomes  obvious  in  the 
quantum  mechanical  formulation  of  the  problem.  A  possible  competing  process  for  the  decay  of 
the  AB*  complex  is  the  radiative  transition  to  lower-lying  states  of  the  molecule.  The  importance 
of  this  effect  can  be  estimated  from  calculations  of  the  band  transition  probabilities. 

Implicit  in  this  model  of  dissociative-recombination  processes  is  the  existence  of  three 
excited  electronic  states  of  AB.  Since  this  recombination  process  is  envisioned  as  occurring  in  two 
stages,  we  are  concerned  with  only  two  of  these  states  at  a  time.  The  state  common  to  both  stages 
of  the  recombination  process  is  the  highly  excited  Rydberg  state  (AB*),  which  is  bound  with 
respect  to  dissociation  to  AB'*’  and  a  free  electron.  For  the  calculation  of  the  initial  rate  of 
recombination,  a  wavefunction  for  a  continuum  Rydberg  state  is  needed.  This  wavefunction  has 
the  asymptotic  behavior  of  a  free  electron  moving  in  the  field  of  AB'*'. 

To  determine  the  transition  rates  in  the  two  stages  of  the  recombination  process,  it  is 
necessary  to  identify  the  terms  in  the  Hamiltonian  that  are  responsible  for  the  transition.  The 
Hamiltonian  for  the  complete  molecular  system  with  respect  to  the  center  of  mass  is : 


K  =  3fco-  — V|-  — XZVV 

1  '  } 


2M 


(51) 
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where  R  is  the  vector  joining  nucleus  B.  3€o  is  the  Hamiltonian  for  the  system  when  the  nuclei  are 
held  fixed,  fi  is  the  reduced  mass  of  the  two  nuclei,  M  is  the  sum  of  the  nuclear  masses,  and  the 
summations  extend  over  all  of  the  electrons.  Since  all  of  the  electronic  eigenfunctions  under 
discussion  are  implicitly  dependent  on  the  intemuclear  separation  |  R  |  ,  it  is  the  term 

=  -  {tp-/24i)V\  which  plays  a  dominant  role  in  coupling  the  electronic  and  nuclear  motion 

giving  rise  to  the  appropriate  transition  (References  42-45).  This  becomes  more  obvious  when  we 
write  the  total  wavefunction  for  the  collision  complex,  AB*,  as : 

V'c*=;Cc(?e.|Rl)Fc(R)  (52) 

where  Xc  (Fe,  1 R  | )  is  the  electronic  wavefunction  for  the  appropriate  Rydberg  state  of  AB  at 
intemuclear  separation  |R|  and  Fc(R)  is  the  wavefunction  associated  with  the  vibration  and 

rotation  of  the  nuclei.  In  a  similar  manner,  the  total  wavefunction  associated  with  the  repulsive 
state  is : 

V'r  =  zX?e.|R|)F,(R)  (53) 

where  F,(R)  is  the  wavefunction  for  the  neutral  atomic  products  as  they  separate.  With  the  use 

of  Fermi’s  golden  rule,  the  transition  rate  for  the  decay  of  the  collision  can  be  obtained  by 
evaluating  the  matrix  elements; 

=  (54) 

The  rate  of  recombination  into  the  highly  excited  Rydberg  states  can  be  calculated  in  a  similar 
manner. 

2.5  Dissociative  Attachment  and  Vibrational  Excitation  Calculations 

A  model  that  has  been  quite  successful  in  explaining  the  dissociative  attachment  (DA)  and 
vibrational  excitation  (VB)  of  diatomic  molecules  is  the  resonance  model  in  which  the  projectile 
electron  is  temporarily  trapped  by  the  target  molecule.  The  molecular  anion  (or  the  resonant 
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state)  thus  fonned  has  a  hnite  lifetime  and  it  can  either  autodetach  leading  to  VE  of  the  molecule, 
or  if  the  lifetime  is  sufficiently  large,  it  can  lead  to  DA  forming  a  neutral  atom  and  an  atomic 
anion. 


A  few  points  of  this  model  are  the  following;  first,  the  dissociative  electron  attachment  and 
the  vibrational  excitation  of  the  molecule  are  two  possible  decay  channels,  apart  from  electronic 
excitation  etc.,  resulting  hrom  a  particular  resonance  state.  Thus  a  calculation  of  the  cross  section 
for  the  dissociative  electron  attachment  to  a  molecule  will  provide  resonant  contributions,  of  that 
particular  resonant  state,  to  the  cross  sections  for  vibrational  excitation  of  the  molecule  and  vice 
versa.  Second,  in  explaining  vibrational  excitation  by  this  resonance  model,  it  is  implicitly 
assumed  that  the  transition  from  the  resonant  state  V~(R)  to  the  neutral  state  Vo(R)  after 
autodetachment  is  a  Franck-Condon  transition,  that  is,  an  instantaneous  transition  with  no 
change  in  nuclear  velocities  or  positions.  This  is  a  so-called  local  complex  potential  model.  This  is 
true  only  if  the  energy  of  the  projectile  electron  is  much  greater  than  the  vibrational  spacing.  At 
very  low  impact  energies  or  if  the  vibrational  spacing  of  the  molecule  is  relatively  large,  a 
description  of  the  DA  and  VE  processes  using  a  nonlocal  complex  potential  for  the  resonant 
state  is  essential  since  the  neutral  molecule,  for  its  vibrational  excitation,  can  accept  only  quanta 
of  vibrational  energy.  Third,  in  the  cases  of  some  molecules,it  might  be  possible  for  certain 
electron  energies  to  form  more  than  one  intermediate  resonant  state.  Alternatively,  the  resonant 
anion  state  may  decay  into  more  than  one  electronic  state  of  the  neutral  molecule.  In  such  cases, 
the  total  width  F  (R)  is  the  sum  of  various  partial  widths,  each  partial  width  corresponding  to  a 
certain  transition  between  the  resonant  anion  state  (or  states)  and  the  neutral  molecular  state  (or 
states). 


2.5.1  Dissociative  Attachment 

The  theory  of  the  capture  of  an  electron  through  the  process  of  dissociative  attachment  (DA) 
to  a  molecule  parallels  that  develo{)ed  above  for  the  dissociative  recombination  process.  The 
process  of  dissociative  attachment  to  a  neutral  molecule,  XY,  can  be  written  as : 


e  +  XY- 


X 

X- 


+  Y- 
+  Y 


(55) 


where  X,  Y  are  the  dissociation  fragments.  The  cross  section  for  dissociative  attachment  is 
obtained  by  comparing  the  flux  of  the  incoming  electrons  per  unit  area  to  the  outgoing 
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ion-neutral  pair  flux  defined  by  the  asymptotic  continuum  wavefunction.  We  find,  for  energy 
normalized  asymptotic  functions: 


oda 


2jih^K  . 

— hm  ^j(R) 


(56) 


where  kj  is  the  wave  number  associated  with  the  incoming  electron  and  —  is  the  relative  velocity 
of  the  asymptotic  ion-neutral  pair.  The  partial  wave,  ^j(R)  satisfies  the  radial  equation: 


=  -ZVvje(R)?vj(R) 


dw-inj  dR'?vj{R')V*(R')Svi(R') 


(57) 


Here: 


Vvje<R)  =  (V^^IE(R)|KcI-E|^R))  (58) 


VV/e(R)  is  the  wavefunction  of  the  neutral  molecule  in  vibrational  state  V  in  the  presence  of  the 
scattering  electron.  0d(R)  is  the  wavefunction  of  the  anion  at  short  intemuclear  separations  as 

described,  for  example,  by  a  stabilization  calculation.  If  the  electron  energy  is  large  or  the 
vibrational  spacing  is  small,  one  can  replace  E  -  Ey  by  the  local  classical  electron  energy  (in  au); 


V-(R)-Vo(R) 


k^R) 

2 


(59) 
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Then  using  the  closure  property; 


X5vj(R)Cvj(R')  =  5(R-R') 

V4 


(60) 


Equation  (57)  reduces  to  a  local  width  model: 


=  5v^)ir^R)/2;tko(R)]>/" 


(61) 


where  Cvi  J  (R)  is  the  wavefiinction  of  the  initial  rovibrational  state  of  the  target  molecule.  The 
resonant  wavefunction  ^j(R)  can  be  obtained  by  numerical  integration  of  Equation  (61)  using,  for 
example,  the  Numerov  technique. 

2.5.2  Vibrational  Excitation 

Vibrational  excitation  of  a  molecule  can  be  written  schematically  as : 


e(k,)  +  XY(V,)-e(kO  +  XY(V,) 


(62) 


When  the  electron  is  far  from  the  molecular  target,  we  have  V[e  +  XY(Vi)]  =  0.  The  initial 
wavefunction  can  then  be  written  as: 


yfi  =  A(ki)e-*'^  Zv, 


(63) 


el 

where  A(ki)  is  the  initial  plane  wave  amplitude,  is  just  the  electronic  wavefunction  of  the 

unperturbed  target.  In  a  similar  fashion,  we  can  write  for  the  asymptotic  wavefunction  after 
vibrational  excitation : 
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ipi  =  Xy, 


(64) 


Conservation  of  energy  requires : 


k; 

E(Vi)  +  ^  = 


E(V,)  +  ^ 


(65) 


Finally  the  cross  section  for  vibrational  excitation  can  be  written  as : 


<iV.-v,  -  II  dki|T„|^ 


(66) 


where  the  transition  matrix  element  from  the  initial  to  the  final  excitation  state  is  obtained  from: 


Tif  = 


4;r^ 

[kikf]*/2 


I  diVf’Ve+XY  V'i 


(67) 


Equation  (67)  can  be  further  decomposed  into  resonant  and  nonresonant  contributions  and  a 
partial  wave  expansion  corresponding  to  the  outgoing  plane  wave. 

A  computer  program  for  calculating  dissociative  attachment  cross  sections  has  been 
developed  based  on  the  local  width  model.  This  model  is  generally  valid  when  the  electron 
collision  energy  is  large  relative  to  the  vibrational  and  rotational  spacings  of  the  neutral  target. 

For  the  problem  of  interest  in  this  research,  the  dissociative  attachment  of  electrons  to  silane,  the 
local  width  approximation  should  be  valid  for  collision  energies  greater  than  - 1.0  eV  This  is 
approximately  five  times  the  vibrational  quantum  of  the  fundamental  mode  in  silane.  An  upper 
limit  is  determined  by  the  convergence  properties  of  the  outgoing  plane  wave  representation  of  the 
dissociation  species,  X'  +  Y  or  X  +  Y~.  In  practice,  this  upper  limit  is  10  to  20  eV,  depending  on 
the  extent  of  the  long-range  potential.  For  dissociative  attachment  at  lower  collision  energies,  the 
nonlocal  width  model  described  by  Bardsley  (Reference  46)  must  be  employed.  In  either  method, 
the  calculation  of  the  complex  part  of  the  interaction  potential  (capture  width)  is  the  most  difficult 
aspect  of  the  theoretical  analysis. 
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2.6  Discussion  of  Charge  Transfer  Calculations 

Low-to-intermediate  energy  ( <  lOOeV)  ion-atom  and  ion-molecule  charge  transfer 
reactions  play  an  important  role  in  many  physical  processes.  Until  recently,  comparatively  little 
effort  was  devoted  to  theoretical  studies  of  the  appropriate  cross  sections.  In  the  past,  both 
theoreticians  and  experimental  physicists  found  it  easier  to  study  high-energy  1  keV) 
collisions.  At  these  energies,  the  two  colliding  particles  preserve  their  identities,  and  it  is  possible 
to  treat  the  interaction  between  them  as  a  perturbation.  There  is  no  guarantee  that  this  procedure, 
known  as  the  Bom  approximation,  will  always  converge  to  the  correct  result  (Reference  47).  As 
the  energy  of  the  colliding  particles  decreases,  it  is  necessary  to  take  account  of  the  distortion 
these  particles  undergo  during  the  collision.  The  method  of  perturbed  stationary  states  was 
developed  for  calculating  charge  transfer  and  electronic  excitation  cross  sections  in  relatively  slow 
collisions  between  heavy  particles  (Reference  48).  The  method  was  presented  in  both  wave 
and  impact  parameter  formalisms.  In  the  first  of  these,  the  entire  system  is  treated  quantum-me- 
chanically,  while  in  the  latter  the  nuclei  are  assumed  to  behave  as  classical  particles,  traveling 
along  straight  line  trajectories,  and  the  time-dependent  Schrddinger  equation  is  solved  to 
calculate  the  probability  of  various  types  of  electronic  transitions  (References  49-51).  Forcing  the 
particles  to  travel  along  straight  lines  limits  the  validity  of  the  impact  parameter  method  to 
collisions  of  several  hundred  electron-volts  or  greater  (Reference  52).  The  wave  formulation  of  the 
method  of  perturbed  stationary  states  appears  to  be  one  practical  method  of  calculating  thermal 
energy  charge-exchange  cross  sections.  A  semiclassical  close-coupling  method  (Reference  53), 
based  on  an  averaging  scaling  procedure,  also  offers  utility  for  low  to  intermediate  collision 
energies. 

Approximate  quanta!  treatments  of  reactive  charge  transfer  at  intermediate  energies  fall 
generally  into  three  basic  categories: 

(a)  Bom  approximation  ignoring  distortion 

(b)  Distorted  wave  Bom  approximations  (adiabatic  perturba¬ 
tions) 

(c)  Quantal  impulse  approximation 

(a)  In  the  Bom  approximation,  distortion  of  the  target  by  the  projectile  ion  is  ig¬ 
nored  and  the  reactive  and  product  wavefunctions  are  replaced  by  their 
asymptotic  forms.  Since  no  explicit  mechanism  for  either  charge  transfer  or 
dissociation  is  possible  within  this  end-state  apprcntimation,  the  method  is  not 
particularly  useful  for  ion-molecule  collisions  at  low  collision  energies.  The 
qualitative  behavior  of  the  Bom  approximation  does,  however,  indicate  that 
the  cross  sections  should  decrease  very  rapidly  with  increasing  collision 
energy. 
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(b)  The  distorted  wave  Bom  approximation  is  based  on  the  exact  solution  for  the 
reactive  cross  section  in  terms  of  the  T-matrix  elements: 

Tpa  =  {Xfi\^p\Xa)  (68) 

where  is  a  solution  of  the  Schrddinger  equation  using  incoming  boundary  conditions  and  Xa 
the  solution  of  the  complete  Hamiltonian  for  the  outgoing  reaction  channel,  including  all  possible 
internal  product  states.  Wp  is  the  reaction  perturbation  potential  which  includes  distortion, 
polarization,  etc.  In  the  distorted  wave  Bom  approximation,  xp  Xa  replaced  by  simpler 

approximations,  usually  corresponding  to  a  single  internal  state  only.  Choi  and  Ikng  (References 
S4  and  SS)  discussed  the  numerical  problems  associated  with  evaluating  matrix  elements 
arising  from  this  approximation  and  several  applications  to  reactive  collisions  have  been  reported 
(References  56  and  57). 

A  still  simpler  approximation  is  available  by  reducing  the  matrix  elements  given  by  Equation 
(68)  to  a  product  of  an  electronic  matrix  element  coupling  the  reactant  and  product  states  and  a 
Franck-Condon  overlap  of  the  conesponding  nuclear  wavefiinctions; 

Tpa=  M,i{xp\Xa)  (69) 

This  so-called  overlap  or  Franck-Condon  model  was  introduced  by  Halavee  and  Shapiro 
(Reference  58)  and  was  described  in  detail  by  Schatz  and  Ross  (Reference  59). 

(c)  In  the  impulse  approximation,  the  operator  of  the  T-matrix  element  given  by 
Equation  (68)  is  replaced  by  an  approximation  which  describes  the  elastic 
scattering  corresponding  to  a  collision  between  the  projectile  ion  and  the  final 
free  particle  of  the  target  molecule.  This  approximation  is  equivalent  to  the 
assumption  that  the  collisions  between  the  particles  are  sudden  so  that  the  col¬ 
lision  time  is  much  shorter  than  the  vibrational  period  of  the  target  molecule. 
Henglein,  et  al.  (Reference  60)  have  argued  that  for  ion-molecule  collisions  of 
several  eV,  this  model,  also  called  the  spectator  stripping  mechanism,  should 
be  valid.  Bates,  et  al.  (Reference  61)  have  described  a  classical  treatment  of 
this  impulse  approximation. 

In  the  following  section,  a  very  powerful  method  for  solving  the  reactive  coupled-channel 
collision  problem  is  described.  This  approach  is  often  referred  to  as  the  R-matrix  or  reactance 
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matrix  propagation  method.  It  generally  yields  a  fast,  efficient  solution  and  automatically 
produces  unitarity  in  the  S-matrices.  By  freezing  open  channels  or  including  the  effects  of  closed 
channels,  several  levels  of  approximation  are  possible  within  the  same  computational  framework. 

2.7  R-Matrix  Propagator  Solution  for  Charge  Exchange  Reactions 

The  theory  of  heavy  particle  charge  exchange  is  complicated  by  the  fact  that,  in  general,  no 
simple  natural  coordinate  system  exists  for  which  a  uniform  description  of  reactant  and  product 
states  can  be  written.  Considering  the  reaction: 

H*+SiH4  -  SiHl  *A  +  H2 

D  ■’I  g 

(70) 

-SiH^  +2H 

the  internal  coordinates  describing  SiH4  cannot  be  uiriformly  mass-scaled  to  represent  SiH3'^ 
formed  as  a  product.  A  natural  coordinate  system  for  such  a  collision  has  been  described  by 
Marcus  (Reference  62),  but  this  yields  very  complicated  forms  for  the  kinetic  energy  operator  and 
presents  many  computational  difficulties. 

Recently,  Stechel,  et  al.  (Reference  63)  have  described  an  approach  in  which  two  separate 
nonorthogonal  frameworks  are  employed  to  describe  reactant  and  product  states.  The  total 
scattering  wavefunction  is  written  as; 


y>  =  V^R.ea)  +  V'’'(R.ey) 


(71) 


where  a  and  y  represent  reactant  and  product  states,  respectively.  Expanding  each  as  a  linear 
combination  of  wavefunctions  of  the  type: 


=  2;^(R)0](Pi) 


(72) 


they  show  that  a  unique  connection  between  reactant  and  product  states  can  be  achieved, 
provided  the  basis  functions,  {0'(pi)}  are  taken  to  be  nonorthogonal,  i.e.,  S"’'  =  <  >  ?«  0. 
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We  note  that  Equation  (72),  in  molecular  language,  is  a  product  wavefunction  of  the 
valance-bond  type.  We  also  note  that  an  SCF  or  MCSCF  framework  for  Equation  (72)  may  not  be 
appropriate  for  treating  the  collisional  aspects  of  the  reactant  channel  due  to  the  orthogonality 
constraints  that  must  be  imposed  on  such  wavefunctions.  The  use  of  nonorthogonal  basis 
functions  is  simplified  through  the  application  of  projection  algebra  techniques  for  the  construc¬ 
tion  of  wavefunctions  of  overall  symmetry.  Efficient  computer  codes  for  implementing  such 
techniques  were  developed  this  Rf'search  Center  over  the  past  years  and  can  be  used  in 
constructing  the  proper  asymptotic  reactant  and  product  state  channels. 

The  solution  of  the  reactive  charge  exchange  problem,  including  the  product  dissociation 
channels,  can  be  accomplished  using  a  direct  Numerov  integration  of  the  Schrodinger  equations 
arising  from  the  basis  representation  of  Equation  (72).  Alternatively,  the  R-matrix  propagator 
method,  originally  formulated  by  Diestler  (Reference  64),  can  be  employed.  This  approach  relates 
the  known  asymptotic  solution  in  the  reactant  channel  to  that  in  the  product  chaimel  through  a 
series  of  recursive  steps.  Defining  the  normal  derivative  to  the  surface  R  =  constant  by 
^nip(R,  q),  we  seek  a  propagator  operator  of  the  form: 


V<R,p)  *  RVaV<R,e) 


(73) 


The  explicit  matrix  representation  of  this  operator  can  be  written  as: 


V’ 


=  -t\Vnyf 


+  r  ^  VnV' 


(74) 


V' 


=  -r'jVnV' 


+  r  i  VoV' 


(75) 


Stechel,  et  al.  (Reference  63),  have  given  explicit  e^qjressions  for  constructing  the  components  of 
this  operator,  both  in  the  asymptotic  and  overlap  regions  of  the  interaction.  The  initial  boundary 

condition  in  the  reactant  channel  is  R'  =  1/^i ,  where  A,  is  the  eigenvalue  of  the  reactant  target 

(in  this  case,  SiH4).  The  accuracy  in  developing  the  final  R-matrix  depends  on  the  chosen  step 
size,  hi,  which  is  taken  as  an  adjustable  parameter  that  can  be  optimized  during  the  course  of  the 
propagator  solution. 
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We  have  recently  implemented  this  R-propagator  analysis  at  UTRC.  using  a  modified 
version  of  the  code  developed  by  Stechel,  et  al.  (Reference  63).  One  application  was  the 
low-energy  collisional  de-excitation  reaction; 


N*  +  O  -N  '‘s  +  O 


(76) 


which  occurs  along  a  curve-crossing  potential  surface.  The  results  of  this  study  have 

been  reported  (Reference  65).  A  similar  study  of  the  kinetics  of  the  charge  transfer  reaction, 

O"*"  +  N2  NO"*"  +  N,  has  been  carried  out  (Reference  66). 

In  summary,  a  description  of  the  energetics  of  charge  transfer  reactions  in  a  valence-bond 
framework,  coupled  with  an  analysis  of  the  kinetics  of  these  collisions  using  a  nonorthogonal  basis 
for  reactant  and  product  channels,  within  the  R-propagator  method  of  solution,  represents  an 
optimum  analysis  for  the  low  to  intermediate-energy  kinetics. 

2.8  Collisionally  induced  Transitions 

Consider  the  charge  transfer  process: 


+  B  — C+  +  D 


(77) 


where  A"*"  is  a  charged  atomic  or  molecular  ion  reactively  scattered  by  a  target  gas  B  to  yield 
products,  C'*'  D.  If  the  reactant  and  product  channels  can  be  separately  decoupled,  it  is  often 
possible  to  treat  Equation  (77)  as  a  pseudo-curve  crossing  problem.  For  near-adiabatic  collisions, 
this  type  of  system  can  be  handled  using  a  semiclassical  approach,  directly  solving  the  second 
order  differential  equations  which  couple  the  adiabatic  levels,  possibly  of  different  angular 
momentum. 

The  adiabatic  states  R)  and  R)  are  eigenstates  of  the  nonrelativistic 
Hamiltonian,  3€ ,  for  a  fixed  intemuclear  separation,  R.  We  have; 
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3t0i(r,R)  =  fi^r.R);  i  =  1,2 


(78) 


In  the  two-state  approximation  we  have,  using  impact  parameter  formulation : 


=  Ci(t)0i(r,R)  e 


(79) 


+  C2(t)02(r.R)e 


where  R  -  R(t)  and  is  defined  by  the  collision  trajectoiy. 

We  now  require  the  inner  product  of  ^  and  ^  to  vanish  over  the  electronic  coordinates. 


( (p2\X-ih—\t/))  =  0 


(80a) 


(80b) 
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Combining  Equations  (78),  (79),  and  (80),  we  have; 


‘dCi  ^  1  #1 

-^  +  C,  <#,|^> 


at 


.  302  -i  f 
+  C2  <  01 1  —  >  e 
3t 


(81a) 


+  Cl  <  02l^  >  e 


(81b) 


Converting  from  differentiation  with  respect  to  time,  to  velocity  and  angular  momentum,  we  have; 


^_dR^  dV» 

at  ~  at  aR  A  at  dt 


LT0i 


(82) 


where  Vo  is  the  angular  velocity  of  the  intemuciear  axis  and  Lj  is  the  angular  momentum 
coupling  operator.  The  first  term  on  the  right  hand  side  of  Equation  (82)  leads  to  the  well  known 
Landau-Zener  (Reference  67)  solution  for  states  of  identical  molecular  symmetry.  For  such  cases 
the  second  term  in  Equation  (82)  vanishes. 

For  the  systems  under  study  here,  where  the  second  term  on  the  right  hand  side  of  Equation 
(82)  contributes  in  Equation  (81),  we  are  led  to  the  following  coupled  equations: 
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dCi  i  dVfl  ^ 

—  + --^C2  <  ^i|Lrl^  >  e  Jo  =0 


(83a) 


dC2  ,  i  SVe  , .  ,  “if  (<i-<2)dT 

—  +  --—Cl  <  ^|Lt|0i  >  e  J» 
dt  ft  dt 


=  0 


(83b) 


Assuming  a  linear  dependence  of  f  on  R  near  the  crossing  point  of  the  collision,  we  have: 


€2-€i  =  bA(R-Rx)  =  at 


(84) 


1^  (85) 

fi  dR 

Assiuning  also  that  <0i|Lr|^>-<Lr>  is  essentially  constant  over  the  dominant  region  of 
the  collision,  we  have : 

dci  iw  <  Lr  >  "if  „  (^) 

— -  + - - — C2e  Jo  =0 

at  ft 

^  +  .i^_<LT>  ^  0  (86b) 

dt  ft 

Equations  (86a)  and  (86b)  were  solved  numerically  by  Russek  (Reference  68)  for  the  case  of 
curve-crossing  states  with  AA  =  ±  1.  For  small  values  of  angular  velocity,  Russek  shows  that 
Equation  (86)  reduces  to  a  standard  Landau-Zener  form.  For  large  velocities,  the  general  solution 
of  Equation  (86)  must  be  employed. 

When  more  than  one  product  channel  is  possible,  a  matrix  of  coupled  equations  of  the  type 
given  by  Equation  (83)  arises  and  a  multisurface  mechanism  is  required  for  a  complete  descrip¬ 
tion  of  such  systems.  Often  a  single  reaction  channel  dominates  and  electron  or  H~  transfer 
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occurs  with  nearly  unit  probability  along  a  single  reaction  surface.  For  such  systems  where  a  curve 
crossing  is  evident,  the  harpoon  (Reference  69)  or  spectator  stripping  model  (Reference  70)  is 

often  useful.  This  leads  to  the  simplest  reaction  cross  section,  a  =  jiR^,  where  R  is  the  effective 
separation  for  a  crossing  on  the  reaction  path.  In  a  next  higher  order  approximation,  the  reaction 
cross  section  is  given  by : 

a  =  2jr  J  Pjj(b)bdb  (87) 

where  Pij  is  the  single  crossing  transition  probability  ( ^  1.0)  which  can  often  be  approximated  by 

the  Landau  -  Zener  formula.  An  extension  of  this  simplified  treatment  to  include  internal 
vibrational  states  has  been  given  by  Bauer,  et  al.  (Reference  71). 

2.9  Classical  Ion-Molecule  Collisions:  Ute  Langevin-Gioumousis-Stevenson 
Reaction  Rate 

In  the  simplest  development  of  ion-molecule  kinetics,  it  is  often  assumed  that  the  reaction 
will  occur  at  the  first  collision  between  the  reactant  pair,  provided  that  there  is  at  least  one 
exothermic  reaction  chaimel  available.  Simple  ion-molecule  reaction  rates  are,  therefore, 
equated  with  ion-molecule  collision  rates.  Experimentally,  many  ion-molecule  reactions  do  occur 
at  every  collision  between  the  ion-molecule  pair,  but  there  are  many  exceptions  to  this  simple 
case.  It  is  now  established  that  in  many  cases,  ion-molecule  collisions  may  be  unreactive  even  if 
an  exothermic  channel  is  available  (Reference  72).  Competing  reaction  channels  or  the  presence  of 
barriers  in  the  reaction  path  (ex:  0'*’  +  N2)  may  also  cause  a  drastic  reduction  in  the  collisional 
reaction  rate  (References  73  and  74). 

An  elementary  treatment  of  ion-molecule  reaction  kinetics  can  be  developed  from  a  classical 
description  of  the  scattering  cross  section.  For  scattering  by  a  spherical  potential,  the  reaction 
cross  section  can  be  written  as ; 

oiE)  =  2ji(  P(b,  E)bdb  (88) 

Jo 

where  b  is  the  classical  impact  parameter  and  P(b,  E)  is  the  probability  that  the  reaction  will 
occur  for  a  given  energy,  E,  and  impact  parameter,  b.  For  ion-molecule  reactions,  a  model  can  be 
constructed  based  on  the  assumption  that  for  values  of  the  impact  parameter  less  than  a  critical 
orbiting  parameter,  bo,  the  reaction  probability  can  be  taken  as  unity.  This  critical  orbiting  impact 
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parameter,  is  defined  so  that  ion  orbits  for  which  b  <  bg.  spiral  in  close  to  the  target  and  orbits 
for  which  b  >  bg  come  no  closer  than  bg/  Jl .  For  the  simplest  possible  case  of  a  spherical 

molecule  with  isotopic  polarizability  and  with  no  permanent  moment,  interacting  with  a  point 
charge,  the  interaction  potential  takes  the  form: 


V(R)  = 


In  this  case  the  orbiting  parameter,  bg,  is  found  to  be: 


(89) 


4e^a 

- 


(90) 


where  p.  is  the  reduced  mass  of  the  colliding  pair  and  v  is  their  initial  relative  velocity.  Using  bg  as 
an  upper  integration  limit  in  Equation  (88),  we  find  for  the  reaction  cross  section: 


,  ,  ,1  2n 

(j(y)  -  ;rbs  =  —  — 

V  L  u  J 


1/2 


(91) 


Equation  (91)  was  restated  by  Gioumousis  and  Stevenson  (Reference  75)  in  terms  of  the 
energy,  E,  of  the  primary  ion  of  mass,  m,  as ; 


[2e2am]  1/2 

' ''  L-tt-J  <’ 

which  shows  an  inverse  square  root  behavior  of  the  reaction  cross  section  with  collision  energy. 
Assuming  a  Maxwellian  distribution  function,  f(v),  for  both  the  ions  and  neutral  molecules,  the 
rate  coefficient  conespmnding  to  the  cross  section  given  by  Equation  (91)  becomes: 


k  =  I  f(v)c7vdv  =  2ji  ^  cm^/sec  (93) 

Equation  (93)  is  the  well-known  Langevin  formula  (Reference  76)  for  the  reaction  rate,  valid  for 
close  encounters  and  typically  of  the  order  of  10"’  cm^/sec.  The  validity  of  Equation  (93)  is 
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strengthened  for  systems  where  the  polarizability  of  the  molecule  is  large  and  the  reduced  mass  is 
small,  so  that  few  collisions  occur  outside  of  the  orbiting  radius.  The  Langevin  rate,  which  is 
predicted  to  be  independent  of  the  collision  energy,  is  thus  most  valid  for  collisions  with  small 
target  gases  such  as  He"*"  +  H2.  For  large  target  molecules  and  ion-molecule  systems  with  many 
internal  degrees  of  freedom,  including  target  dissociation,  the  Langevin  rate  will  predictably  be  too 
large.  This  simple  model  also  is  incapable  of  predicting  the  product  distribution  of  an  ion-molecule 
reaction.  It  should  be  noted  that,  in  this  regard,  the  Langevin  formula  is  formally  invalid  since  the 
requirement  of  unitary  in  the  scattering  matrix  (principle  of  detailed  balancing)  is  violated  for  this 
classical  treatment. 

2.10  Ab  Initio  Theoretical  Thermochemistry  Models 

The  development  of  a  predictive  ab  initio  thermochentistry  model  that  is  reliable  to  1-2 
kcal/mol  has  been  the  goal  of  computational  chemistry  for  some  time.  Early  models  by  Pople  and 
coworkers  were  based  on  calculated  energy  differences  at  the  MP4  level  of  theory  for  prescribed 
isogyric  processes  (References  77-79).  A  related  bond  additivity  correction  (BAG)  method 
proposed  by  Binkley  and  Melius  (Reference  80)  uses  calculated  and  experimental  atomization 
energies  for  a  set  of  reference  compounds  to  compute  empirical  correction  factors  for  each  bond 
type.  Ignacio  and  Schlegel  (Reference  81)  report  heats  of  formation  for  SiHn,Fn  molecules  based 
on  the  calculated  energetics  of  isodesmic  reactions  at  the  MP4  level  of  theory  using  extended  basis 
sets. 


More  recently,  a  systematic  procedure  for  prediction  of  molecular  energies,  referred  to  as 
Gaussian-1  (Gl)  theory,  has  been  proposed  by  Pople  et  al.  (Reference  82).  In  this  method, 
systematic  corrections  for  both  basis  set  deficiencies  and  the  level  of  correlation  treatment  are 
applied  to  the  energies  calculated  at  the  QCISD(T)/6-311G(d,p)//MP2(FU)/6-31G(d)  level  of 
theory.  The  Gl  energy  can  be  written  as: 


Eo(Gl)  =  EQcisD(D/6-311G(d,p)  +  tEMP4(FC)/6-311+G(d,p)  ”  EMP4(FC)/6-31IG(d,p)l 

[^MP4(FC)/6-311G(2df,p)  “  EmP4(FC)/6-3  1  lG(d,p)]  ~  0.00614  npgu-  (94) 
-  0.00019  nos +  E,pe. 


A  suggested  improvement  to  the  method,  referred  to  as  Gaussian-2  (G2)  theory  (Reference  83), 
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further  corrects  for  nonadditivity  effects  in  the  treatment  of  the  basis  set: 


Eo(G2)  =  Eo(Gl)  +  [EMP2/6.311+G(3df^p)  "  EMP2/6-311G(2df.p)] 

~  t^MP2/6-311+G(d,p)  ~  ^MP2/6-311G(d,p)l  0-001 14  Ppajp  (95) 

In  the  present  study  we  systematically  examined  the  thermochemistry  of  SiH„’  and  SiF„"  species 
using  both  the  G1  and  G2  approaches.  We  then  compared  our  results  with  previous  calculations 
on  the  neutral  molecules  and  the  limited  experimental  data  on  electron  affinities  of  these  molecules. 
Our  conclusion  is  that  the  G1  method  is  accurate  for  Si-containing  anions  to  ±  3  kcal/mol.  The  G2 
method  appears  to  be  less  reliable  for  the  treatment  of  molecules  with  multiple  similar  bonds  such 
as  in  SiF4. 

The  predicted  G1  and  G2  energies  were  used  to  predict  heats  of  formation  using  two 
different  approaches.  In  the  first  approach  all  data  are  referenced  to  the  calculated  total  atomization 
energies.  For  SiX^/SiXo"  species,  we  have: 


EDo  =  E(SiX„)  -  E(Si)-rE(X) 


(96) 


The  experimental  heats  of  formation  of  Si(g),  H(g)  and  F(g)  are  used  as  reference  states.  In  the 
second  approach  an  isodesmic  reaction  such  as: 


n/4  SiX4  (g)  +  (4-n)/4  Si  (g)  SiXn  (g)  (97) 


is  used  to  derive  heats  of  formation  based  on  the  experimental  heats  of  formation  of  the  parent 
molecule.  Both  approaches  were  employed  in  these  studies. 
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SECTION  3 


DISCUSSION  OF  RESULTS 


The  main  fcKus  of  our  current  studies  was  on  understanding  the  role  that  negative  ions  play 
in  the  discharge  chemistry  pertinent  to  deposition  of  amorphous  silicon  thin  films.  To  this  end,  a 
comprehensive  analysis  of  the  thermochemistry  of  negative  silane  ions  was  first  undertaken.  We 
find  that  significant  errors  exist  in  the  thermochemical  data  base  for  SiF^  and  SiFn"  species. 
Negative  ion  clustering  reactions  were  then  examined  to  uncover  possible  pathways  for  the  build¬ 
up  of  successively  larger  negative  ion  clusters  through  association  with  SiH4  and  elimination  of 
H2.  We  found  such  pathways  to  be  energetically  unfavorable  and  conclude  that  addition  reactions 
to  form  negative  ion  clusters  are  probably  not  of  major  importance.  An  analysis  of  excited 
molecular  state  surfaces  for  SiH4"  indicated  several  pathways  for  dissociation  attachment  of  e  + 
SiH4  to  occur  in  the  7  - 10  eV  region.  No  low  energy  pathways  could  be  located  although  several 
autoionizing  resonance  states  were  found.  Finally,  the  energetics  and  possible  reaction  pathways 
for  SiH5’  and  SiF5"  formation  were  examined.  SiH5"  and  SiF5‘  are  found  to  be  stable  gas  phase 
anions  and  their  role  in  the  kinetics  of  silane  plasmas  is  examined.  A  summary  of  our  results  of 
these  investigations  follows. 

3.1  Thermochemistry  of  Sil^/SiHn"  and  SiFn/SiFn"  Species 

The  first  part  of  our  study  of  negative  ion  silane  reactions  was  concerned  with  an  examination 
of  the  electronic  structure,  predicted  equilibrium  geometries  and  thermochemistry  of  both  neutral 
and  ionic  silane  species.  The  thermochemistry  of  negatively  charged  species  of  silicon,  hydrogen 
and  fluorine  is  at  present  uncertain,  owing  to  a  sparsity  of  experimental  data  and  the  lack  of  a 
systematic  theoretical  study.  Reliable  heats  of  formation  of  SiHn"  and  SiF^"  species  are  of 
considerable  interest  since  theory  may  help  to  determine  the  role  that  such  anions  play  in  the 
deposition  of  amorphous  silicon  films  and  in  the  etching  of  silicon  surfaces.  In  contrast  to  positive 
ions,  which  are  swept  to  the  walls  by  the  ambipolar  fields  in  a  typical  plasma  assisted  chemical 
vapor  deposition  (CVD)  device,  negative  ions  tend  to  become  trapped  with  resultant  long  residence 
times.  Garscadden  (Reference  9)  has  suggested  that  negative  ions  may  form  clusters  through  a 
series  of  ion-molecule  reactions,  similar  to  that  discussed  experimentally  by  Mandich,  Reents  and 
Jarrold  (Reference  11),  and  theoretically  by  Raghavachari  (Reference  12),  for  positive  ion  cluster 
growth.  It  was  not  clear  whether  sequential  clustering  reactions  of  negative  silicon  hydride  and 
silicon  fluorine  ions  could  occur,  as  has  been  demonstrated  with  positive  ions,  and  the  relative 
importance  of  such  anions  in  silicon  plasmas  is  uncertain. 

In  this  research  a  systematic  evaluation  of  the  geometry  and  energetics  of  SiHn"  SiPp" 
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anions,  and  their  corresponding  neutral  parents,  was  carried  out  The  thennochemistry  of  neutral 
SiHn  and  SiFn  species,  and  their  corresponding  cations,  has  been  the  subject  of  several 
experimental  (References  84-90)  and  theoretical  studies  (References  91-100).  In  our  work,  we 
extended  these  studies  to  include  the  thermochemistry  of  the  anions  in  order  to  compile  a  consistent 
set  of  heat  of  formation  data  for  these  compounds. 

3.1.1  Theory 

The  development  of  a  predictive  ab  initio  thermochemisoy  model  that  is  reliable  to  1-2 
kcal/mol  was  the  goal  of  computational  chemistry  for  some  time.  Early  models  by  Pople  and 
coworkers  were  based  on  calculated  energy  differences  at  the  MP4  level  of  theory  for  prescribed 
isogyric  processes  (References  77-79).  A  related  bond  additivity  correction  (BAC)  method 
proposed  by  Binkley  and  Melius  (Reference  80)  uses  calculated  and  experimental  atomization 
energies  for  a  set  of  reference  compounds  to  compute  empirical  correction  factors  for  each  bond 
type.  Ignacio  and  Schlegel  (Reference  81)  report  heats  of  formation  for  SiH^Fn  molecules  based 
on  the  calculated  energetics  of  isodesmic  reactions  at  the  MP4  level  of  theory  using  extended  basis 
sets. 


More  recently,  a  systematic  procedure  for  prediction  of  molecular  energies,  referred  to  as 
Gaussian- 1  (Gl)  theory,  was  proposed  by  Pople  et  al.  (Reference  82).  In  this  method,  systematic 
corrections  for  both  basis  set  deficiencies  and  the  level  of  correlation  treatment  arc  applied  to  the 
energies  calculated  at  the  QCISD(T)/6-311G(d,p)//MP2(FU)/6-31G(d)  level  of  theory.  The  Gl 
energy  can  be  written  as: 


Eo(Gl)  =  EQciSD(T)/6-311G(d.p)  IEMP4(PC)/6-31I+G(d,p)  “  EMP4(FQ/6-311G(d,p)] 

tEMp4(FC)/6-3UG(2df.p)  “  EMP4(FC)/6-311G(d.p)I  “  0.00614  npgir  (98) 

—  0.0(X)19  Pqj  +  E^pg. 

A  suggested  improvement  to  the  method,  referred  to  as  Gaussian-2  (G2)  theory  (Reference  83), 
further  corrects  for  nonadditivity  effects  in  the  treatment  of  the  basis  set: 

Eo(G2)  =  Eo(Gl)  f  [EMP2/6-311+G(3df.2p)  “  EMP2/6-311G(2df4))3 

-  fE^|p2/(i.31,+G((j,p)  -  EMP2/6-3nG(d,p)]  +  0.001 14  npajp  (99) 

In  the  present  study  we  systematically  examine  the  thermochemistry  of  SiHn*  and  SiFn"  species 
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using  both  the  G1  and  G2  approaches.  We  then  compare  our  results  with  previous  calculations  on 
the  neutral  molecules  and  the  limited  experimental  data  on  electron  affinities  of  these  molecules. 

3.1.2  Calculated  Results 

Optimized  molecular  structures  for  all  stable  SiHn/SiHn"  and  SiFn/SiFn"  species  were  found 
at  both  the  HF/6-31G(d)  and  MP2/6-31G(d)  levels  of  theory.  The  HF  results  were  used  as  a  fu^t 
trial  for  the  MP2  optimization  runs  and  as  estimates  of  the  zero-point  energy  corrections.  The 
calculated  optimized  geometries  are  given  in  Table  1.  The  predicted  geometries  at  the  MP2/6- 
31G(d)  level  agree  with  experimental  bond  lengths  to  ±  0.  002  A  and  with  experimental  bond 
angles  to  ±  0.4°,  for  SiH^  molecules.  For  the  SiF^  systems,  the  errors  are  somewhat  larger;  ± 

0.03  A  for  bond  lengths  and  ±  0.7°  for  bond  angles.  Overall  the  agreement  with  experiment  is 
very  good  and  follows  the  observed  pattern  that  MP2  generally  overestimates  bond  lengths  while 
the  HF  level  of  theory  typically  predicts  bond  lengths  shorter  than  experiment.  The  optimized 
MP2  geometries  given  in  Table  1  would  appear  to  be  reliable  estimates  for  carrying  out  the  energy 
correlation  calculations  required  for  the  G1  and  G2  methods. 

The  calculated  energies  for  the  G1  method  are  given  in  Table  2.  We  include  several  SiHn  and 
SiHn"  compounds  previously  reported  by  Curtiss,  et  al.  (References  82  and  83)  for  completeness; 
agreement  is  found  to  ±  0.01  mh  in  the  calculated  energies.  The  MP4/6-31  lG(2df)  contribution  to 
the  G1  energy  for  SiF4  represents  one  of  the  largest  calculations  in  snidies  by  this  method.  The 
nonadditivity  corrections  for  the  G2  method  are  given  in  Table  3.  A  comparison  of  the  G 1  and  G2 
methods  for  over  100  compounds  indicates  that  the  G2  method  offers  an  average  improvement 
over  the  G1  method  of  0.4  kcal/mol  in  calculated  total  atomization  energies  (Reference  83). 

The  predicted  G1  and  G2  energies  were  used  to  predict  heats  of  formation  using  two 
different  approaches.  To  facilitate  comparison  with  the  calculated  thermochemistry  reported  by 
Ignacio  and  Schlegel  (Reference  8 1),  in  the  first  approach  all  data  were  referenced  to  the  calculated 
total  atomization  energies: 


IDo  =  E(SiX„)  -  E(Si)-nE(X) 


(100) 


The  experimental  heats  of  formation  of  Si(g),  H(g)  and  F(g)  were  used  as  reference  states.  In  the 
second  approach  the  isodesmic  reaction: 


n/4  SiX4  (g)  +  (4-n)/4  Si  (g)  SiXn  (g) 


(101) 
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Table  1 .  Optimized  Molecular  Geometries  for  SIHp/SIHp"  and  SIFn/SiF, 


t 


o 

CB 

*i* 

(0 


in 

CJ 

IT) 

CO  CO 

CM 

CO 

05 

to 

t'- 

05 

lO 

If) 

CO  to 

O 

CO 

CM 

CO 

CO 

CM 

CO 

05 

CO 

00 

CO 

csi 

cd 

yml  • 

•  '  cvi  T-' 

I  I  d 

in 

od  '  ' 

I  I  r.; 

cd 

05 

T“ 

▼- 

05  CM 

o 

t— 

o 

05 

CM 

CM 

^  U) 

05 

to 

05 

o 

to 

CM 

CO  CM 

00 

00 

r- 

(O 

CO 

CO 

O 

O  CO 

If) 

05 

y— 

o 

LO 

cd 

od 

cd  ' 

•  •  cd  ^ 

■  i  d 

cd 

Y*  1  1 

•  '  1^ 

cd 

c:5 

■»— 

o 

05  CM 

05 

■t— 

T- 

05 

CM 

T— 

T“ 

o 

CO 

CM 

T— 

CM 

CO 

'M' 

CM 

r- 

00 

y^ 

O 

(35 

CO 

m 

CM 

r^ 

in 

in 

CO 

in 

in 

in 

CO 

d 

T— 

y~ 

T— 

r— 

in 

C7) 

CO 

O 

<Ji 

00 

CO 

CO 

to 

to 

o 

CM 

o 

o 

CO 

r*. 

in 

00 

r- 

CO 

CM 

o 

05 

O 

CM 

o 

h- 

CO 

h- 

<75 

T— 

05 

o 

00 

CO 

CJ 

T“ 

00 

00 

00 

CO 

in 

in 

CM 

CM 

T— 

<35 

00 

CM 

o 

<35 

r>- 

in 

in 

Tf 

in 

in 

in 

in 

CO 

to 

CO 

in 

It) 

CO 

r>. 

cq 

cq 

d 

y— 

T“ 

T“ 

o 

CM 

CO 

05 

05 

CO 

CO 

m 

to 

(75 

h>. 

to 

o 

CO 

o 

00 

y^ 

CM 

<75 

o 

o 

in 

00 

in 

in 

CM 

CM 

CO 

<J5 

CM 

CO 

in 

t>. 

00 

y~ 

CM 

CO 

CO 

•r~ 

o 

t'- 

h- 

in 

in 

CO 

'«*• 

o 

<35 

05 

r*. 

in 

o 

00 

00 

in 

in 

in 

in 

in 

in 

in 

CO 

CO 

in 

in 

1(3 

in 

CO 

cq 

cq 

cq 

cq 

d 

T- 

T- 

•r~ 

y- 

■r^ 

ir- 

T- 

n— 

y- 

T- 

T~ 

1—" 

y^ 

T- 

00 

O 

00 

CO 

in 

CO 

CO 

O 

CM 

CM 

in 

<75 

in 

CO 

CO 

CO 

rv. 

CO 

00 

r- 

CM 

y— 

T— 

00 

05 

CO 

00 

00 

CM 

o 

O 

05 

CO 

CM 

CM 

<35 

(O 

CM 

CO 

CO 

in 

in 

CM 

00 

in 

CO 

CO 

o 

CM 

in 

05 

00 

in 

00 

T“ 

05 

<75 

o 

CO 

00 

CM 

O 

00 

CO 

00 

CO 

05 

CO 

f-- 

q 

o 

CO 

q 

'M- 

O 

O 

q 

q 

CM 

<75 

q 

q 

CM 

q 

C75 

d 

d 

d 

T— 

<75 

05 

d 

d 

<35 

od 

od 

od 

od 

cd 

od 

od 

od 

od 

1 

00 

<75 

<75 

<75 

<75 

00 

00 

<35 

<75 

<75 

00 

00 

CD 

00 

00 

CO 

00 

CO 

00 

00 

1 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

y- 

CO 

■M- 

in 

to 

CO 

CO 

in 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

CO 

CO 

00 

05 

CM 

CO 

CO 

<75 

O 

CD 

CM 

CO 

r— 

o 

00 

in 

rv 

00 

00 

in 

y— 

00 

CM 

CO 

CO 

CM 

00 

CM 

<35 

00 

CM 

CO 

<75 

<75 

CM 

ro 

in 

00 

o 

T— 

CO 

CM 

•M- 

(75 

m 

CO 

in 

T— 

CO 

CM 

o 

<75 

05 

o 

CM 

CM 

r^ 

o 

r^ 

00 

<75 

00 

00 

CO 

T“ 

05 

05 

CO 

q 

q 

O 

<o 

(£5 

CO 

00 

N- 

CO 

(75 

CO 

q 

q 

00 

T— 

05 

05 

<35 

d 

05 

05 

d 

d 

od 

od 

h*' 

(d 

in 

od 

od 

od 

i< 

1 

00 

00 

00 

C75 

<75 

00 

00 

<75 

<75 

05 

00 

00 

00 

00 

00 

CO 

00 

00 

00 

00 

1 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CO 

in 

CO 

CO 

CO 

in 

1 

1 

.  1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

+ 

D>  ^  ^  T“ 

W  C  <  CQ  <  < 


»-  CM  T-  O  <M  •*- 


I  + 

H  N  OQ  < 

CO  CM  »- 


+  ^  T-  ^  T-  T- 

w°t:  <"  oT  ^  < 


CV  T-  CO  ^ 


I  +  ^  T- 

H  M  CD  < 

CO  »■  CM  T“ 


M 

0> 

•mm 

O 

a> 

a 

(0 


CM  CM  (O  'M' 

X  I  X  I 
X  CO  CO  CO  <0  c/3 


(  • 

II  CM  to 


XXIX 

CO  CO  CO  W 


CM  CM  CO 

^U-  ULU-lLlL 
U.  CO  M  CO  CO  M 


u. 


O 


I  I 

II  CM  CO 


U.  U.  U.  LI. 
CO  M  ^  CO 


50 


Energy  in  hartrees,  bond  lengths  in  A,  angles  in  degrees.  Experimental  geometries  (  )  are  from  Refs.98  and  101. 
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Table  3.  Calculated  Energies  and  Corrections  for  G2  Method 


CM 

O 

UJ 


o 

o 

o 

o 

ir> 

d 

I 


UJ 

o 

o 

CD 

CM 

UJ 

CO 

CO 

CM 

o 

|V 

UJ 

O 

CO 

CO 

CD 

rv 

O 

CO 

CD 

CO 

00 

CO 

CO 

CO 

fv 

T“ 

|V 

CD 

UJ 

|v 

rv" 

OC) 

CD 

d 

d 

d 

1 

CO 

00 

00 

CD 

CD 

CD 

CD 

1 

CM 

CM 

CM 

CM 

CM 

CM 

1 

1 

1 

1 

1 

1 

CM 

CO 

|V 

UJ 

UJ 

UJ 

o 

00 

o 

CO 

00 

CO 

CD 

CO 

UJ 

r- 

00 

CO 

CO 

o 

CM 

UJ 

CD 

UJ 

UJ 

CM 

00 

d 

00 

CD 

ai 

d 

d 

1 

00 

00 

00 

CD 

CD 

1 

CM 

CM 

CM 

CM 

CM 

1 

1 

1 

1 

1 

CM 

CD 

o 

00 

CO 

CD 

y— 

00 

CD 

o 

UJ 

o 

UJ 

CD 

CO 

CM 

CO 

IV 

CO 

CD 

CO 

CO 

CO 

CM 

00 

CO 

rr 

UJ 

UJ 

CD 

CO 

CO 

|V 

CO 

UJ 

't 

CO 

p 

CD 

CD 

00 

OO' 

00 

00 

00 

|V 

CD 

CD 

00 

00 

OO 

00 

00 

CO 

1 

CO 

■M- 

UJ 

CO 

I  I  I  I  I  I  I 


o 

CO 

UJ 

CM 

CM 

CO 

•r— 

CM 

CM 

CD 

o 

rv 

CO 

CM 

CO 

CD 

rv 

UJ 

rv 

00 

rv 

CO 

UJ 

CD* 

(xj 

ocj 

ocj 

ocj 

CD 

00 

00 

00 

00 

1 

CO 

CO 

UJ 

I  I  I  I 


O 

0. 

> 

UJ 

< 


00 

CM 

00 

CM 

CO 

00 

CM 

CM 

CO 

CM 

00 

o 

CO 

CM 

00 

CO 

o 

CO 

CM 

CM 

CM 

UJ 

CM 

Tt- 

UJ 

CD 

rv 

CM 

T- 

CO 

CM 

CM 

p 

rv 

p 

CM 

CO 

o  ^ 

CM 

CO 

CM 

CO 

T* 

CM 

CO 

CO 

CO 

r^ 

uj 

d 

CD 

CO 

00 

ocj 

uj 

CD 

c> 

T” 

rr 

gJCMOOOCO’M-tOOJ 

'M-oocoocMcor^r^ 
oidd^cocMCO'^f 
I  I  I  I  I  '  I  I 


T-  00  00  0>  CD 
O)  CM  <0  CO  00 
T-'  CO  CO  Tf  tfi 
I  I  I  I  I 


g>Ttcooo'^ooeoo> 
codcMuocor^cM'r- 

I  I  I  I  I  I 


I 


CM 

I 


a 

Ql 

0) 

u 

c 

© 

re 

> 

c. 

X 


CO  CO  '<* 
CM  0>  CM 

in  cvi 

I  I  I 


(O  o 

IV  00 
CO  00 

I  I 


O 

'o 

UJ 


o 

o 

o 

o 

lO 

d 

I 


O  CO  00  o 
O  <0  tv  CO 
m  CO  <o 
CD  00  CO 
rv  cj)  UJ 

t-‘  IV  fld  <j> 
I  CO  00  00 
,  CM  CM 

I  I 


O  CO  ■«-  CM 
O)  T-  CM  00 
IV  O  CO  00 
CO  CO  IV  ^ 
T-  -r-  IV  Tf 

d  d  d 

0>  O)  O)  OJ 
CM  CM  CM  CM 

I  I  I  I 


CO  CO  |V 

CD  CM 
CD  CM  00 
00  00 
UJ  OT  UJ 

d 


I 


0>  00 
00  O  CO 
CO  CO 
CO  O  CM 
UJ  CM  OO 

cjj  d  d 
00  OJ  o> 

CM  CM  CM 


00  O) 

00  00 
CM  CM 

111(1 


UJ  O  CM 
|V  CM  CM 
CM  UJ  O 
CO  CM  CD 
CD  CO  |v 

0>  0>  00 
(D  CD  00 
T-  CO 


I 


CO  CD  CO  O  CO 

1-  -M-  00  UJ  CO 

CD  00  UJ  CD  CM 

CO  ^  CO  UJ  CO 

CO  UJ  '<0;  CO  o 

ocj  00  00  00  rv 

ra  «  ^  00  CO 


I  I  I  I  I 


CO 

I 


CO  O  UJ 
OJ  CD  00 
CD  O  00 
UJ  CM  CD 
rv  cp  rv 

oil  00  00 

CD  00  00 
CO  CO 


I 


I 


CM  rr  e 
|v  CD  S 
CD  00  ^ 

iv  in  *3 

CD  UJ  5 
00  00  = 
00  00  c 
UJ  3 


I  I  I  I 


o 

(0 

(A 


(/)>rcLa.G<cD<< 

CMi-  COCOCMt-CJCMt- 


CA  CA  w  W  dr^ 

y-  n  T-  eg  T- 


CL  W*C  <*  CO’  ^  <*  < 

C\Ii-CNI»-COCMt-  r- 


CA  w  w  Cd”  <” 

^  CO  T-  CM  »- 


ID 

£ 

u 

o 

a 

(0 


CM  CM  CO  ^ 

CM  ._????? 

IIOCACACACACACA 


II  CM  CO 
I  '  ^  ^ 

I  W  CA  CA  CA  W 


CM  CO  ^  ^ 
U.  tL  U.  LL  li 

LU  uT*^  ^  CO  M  M  O 


CM 

.U.  U. 


I  I  M 

II  CM  CO  h 

U.  U.  U.  3 

L  W  <A  W  CA  W 


52 


was  used  to  derive  heats  of  formation  based  on  the  experimental  heats  of  formation  of 
SiH4(g)[AHi"(0  K)  =  +10.5  kcal/moll  and  SiF4(g)[AHfO(0  K)  =  -384.8  kcal/moll.  Ignacio  and 
Schlegel  [81]  have  applied  Equation  (101)  to  SiH„  and  SiF„  species  using  MP4/6-31++G(3df,  3p) 
energy  differences.  The  results  reported  here  are  derived  from  G1  and  G2  energy  differences, 
which  include  a  higher  order  coupled-cluster  energy  correction.  Equation  (101)  tends  to  balance 
the  errors  in  SiX4/SiXn  calculations,  whereas  Equation  (100)  reflects  possible  accumulation  of 
error  in  the  successive  atomization  of  SiXn. 

Corrections  to  the  standard  state  temperature  of  298.15  K  were  made  using  the  relationship; 


AHfO  (298  K)  =  AHfO  (0  K)  +  [H°  (298  K)  -  H°  (0  K)]siXn 

-  [H°  (298  K)  -  (0  K)]si  -  n/2  [H°  (298  K)  -  H“  (0  K)]x2  (102) 


JANAF  values  (Reference  98)  were  used  for  the  enthalpy  corrections  for  Si,  H2  and  F2.  The 
enthalpies  of  SiH^  and  SiFn  species  were  corrected  using  an  ideal  gas  treatment.  Zero-point 
energies  were  calculated  using  scaled  HF/6-31G(d)  frequencies.  Our  calculated  heats  of  formation 
for  SiHj^SiHn"  and  SiF,^SiFn”  species  are  given  in  Tables  4  and  5,  respectively. 

3.1.3  Analysis  of  Calculated  Heats  of  Formation 

A  comparison  of  our  predictions  for  the  SiHn  species  given  in  Table  4  indicates  agreement 
with  the  best  experimental  values  to  ±  2  kcal/mol  with  the  single  exception  of  SiH2  where  the  G 1 
and  G2  predictions  appear  to  be  in  error  by  ~  3  kcal/mol  for  the  [lAj]  ground  state.  In  contrast, 
calculations  for  the  excited  I^Bj]  state  of  SiH2  agree  with  the  best  experimental  estimate  to  within 
1  kcal/mol.  Previous  studies  (References  99  and  109)  have  suggested  that  a  source  of  error  for  the 
[lAi]  ground  state  may  be  the  inadequacy  of  the  single  HF  reference  state.  The  [lAj]  state  of  SiH2 
was  shown  to  be  best  represented  using  a  two-configuration  valence  reference  state: 


(2b22  lbi2  5ai2)  +  X(2b22  lbi2  2b, 2) 


(103) 


Use  of  this  two-configuration  reference  state  for  SiH2  showed  (Reference  109)  to  yield  improved 
lA  j  -  excitation  energies  for  a  variety  of  correlation  treatments  that  were  examined.  However, 

Grev  and  Schaefer  (Reference  100)  have  recently  shown  that  the  [lA,]  state  of  SiH2  can  be 
adequately  treated  using  a  single  configuration  reference  state,  provided  that  a  sufficiently  flexible 
basis  is  chosen  and  a  high  degree  of  correlation  [CCSD(T)]  is  employed.  Our  G1  calculations 
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Table  4.  Heats  of  Formation  for  Silicon  /  Hydrogen  Species 
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\SO{G2):  Based  on  iscxtesmic  reactions,  n/4  SiH4  (4-n)/4  Si  SiH^,  calculated  at  the  G2  level  of  theory. 

ISO<IS):  Based  on  isodesmic  reactions,  nr4  SiH4  +  (4-n)/4  Si  -» SiHn,  calcuiated  at  the  MP4/6-31++G(3df,  3p). 

//HF/6-3lG(d)  ievei  (Ref.  81]. 

BAG:  Based  on  atomization  energies  calculated  at  the  MP4/6-31  G{d,  p)//HF/6-31  G(d)  level  with  empirical  corrections  tor  Si-H  bonds.  [Ref.  99] 

GS:  Based  on  atomization  energies  calculated  at  the  CCSD(T)  level  with  a  [65321/4321]  ANO  basis  set  [Ref.  100] 
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indicate  a  singlet-triplet  splitting  for  SiH2  of  23.7  kcal/mol,  which  is  larger  than  the  21.0  kcal/mol 
splitting  reported  by  Berkowitz,  et  al.  (Reference  88).  This  conesponds  to  the  observation  that  the 
G1  and  G2  analysis  overestimates  the  atomization  energy  for  the  ground  state  of  SiH2.  A  full 
QCISDCD/b-Sl  1+G(2df,p)  calculation  yields  nearly  identical  results  to  the  G1  analysis,  indicating 
that  basis  set  nonadditivity  effects  in  the  G1  theory  are  negligible  for  SiH2.  Using  an  improved 
silicon  basis  [6-31 1+G(3df,2p)],  with  a  full  QCISD(T)  correlation  treatment,  lowers  the  singlet- 
triplet  splitting  to  22.8  kcal/mol.  We  conclude  that  the  error  in  the  Aj]  state  of  SiH2  arises 
mainly  from  basis  set  deficiencies. 

Grev  and  Schaefer  (Reference  100)  have  carried  out  an  extensive  analysis  of  the 
thennochemistry  of  SiU^  species.  They  conclude  that  the  experimental  heats  of  formation  of  SiH 
and  SiH4  may  both  be  too  high  by  1-2  kcal/mol.  The  G1  and  G2  results  reinforce  this 
observation,  but  our  calculated  results  are  within  the  average  error  limits  (±  2  kcal/mol)  for  these 
methods.  In  general,  the  calculated  heats  of  formation  based  on  Equation  (101)  are  in  somewhat 
better  agreement  with  experiment  than  those  based  on  total  atomization  energies.  The  calculated 
electron  affinities,  based  on  G1  and  G2  energies,  are  also  shown  in  Table  4  and  compared  with  the 
best  experimental  estimates.  Our  reported  heats  of  formation  for  the  SiH„"  anions  are  based  on  the 
ion  convention  that  excludes  the  enthalpy  of  the  electron.  The  agreement  with  experiment  is  seen 
to  be  very  good  (~  0.1  eV). 

In  Table  5  we  report  the  results  of  a  parallel  study  of  the  SiFn  species.  The  agreement 
between  the  G1  and  G2  predictions  is  ~  3  kcal/mol  for  the  calculated  heats  of  formation  based  on 
total  atomization  energies.  Somewhat  better  agreement  (~1  kcal/mol)  is  found  for  the  G1  and  G2 
estimates  based  on  the  isodesmic  reactions  given  in  Equation  (101). 

The  experimental  data  for  SiF^  species  are  more  uncertain  than  those  reported  for  the  SiHn 
system.  Johnson  (Reference  107)  has  recently  reevaluated  the  heat  of  formation  of  SiF4  as  -386.2 
±0.1  kcal/mol.  Accepting  this  value,  our  G1  and  G2  predictions  arc,  respectively,  3.8  and  7.1 
kcal/mol  too  low.  This  corresponds  to  an  underestimation  of  the  total  atomization  energy 
compared  with  SiH4,  To  shed  some  light  on  the  source  of  the  error  in  the  SiF4  calculations, 
especially  the  G2  predictions,  a  parallel  calculation  of  the  atomization  energy  for  CF4  was  carried 
out.  The  calculated  energies  and  corrections  for  CF4,  applying  the  G1  and  G2  methods,  are 
shown  in  Tables  2  and  3,  respectively.  For  this  case,  the  G1  and  G2  methods  overestimate  the 
total  atomization  energy  of  CF4  (experiment  =  465.47  kcal/mol)  by  3.4  and  5.6  kcal/mol, 
respectively.  We  conclude  that  the  magnitude  of  the  G1  and  G2  correction  terms,  which  are 
comparable  for  CF4  and  SiF4,  are  not  responsible  for  the  observed  errors  in  the  SiF4  case  and  that 
basis  set  deficiencies  are  the  most  likely  problem.  For  SiF3,  we  find  good  agreement  (1  kcal/mol) 
between  the  value  reported  by  Walsh  (Reference  86)  and  our  predictions  ba.sed  on  isodesmic 
reactions.  The  experimental  value  of  -259  kcal/mol  reported  by  Farber  and  Srivastava  (Reference 
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106)  appears  to  be  20  kcal/mol  too  low.  For  SiF2  the  calculated  heat  of  formation  is  -152  ±  2 
kcal/mol  averaged  over  G1  and  G2  predictions  using  both  atomization  and  isodesmic  reaction 
methods.  This  value  is  in  good  agreement  with  that  reported  by  Ignacio  and  Schlegel  (Reference 
81)  (-153  kcal/mol)  and  the  BAG  value  (-150  kcal/mol)  reported  by  Ho  and  Melius  (Reference 
91).  The  experimental  heat  of  formation  for  SiF2  seemed  to  be  established  as  -140.5  kcal/mol 
(References  87  and  98),  a  value  10-12  kcal/mol  more  endothermic  than  that  predicted  from  theory. 
A  new  experimental  study  of  SiF2  appears  to  be  warranted  since  the  several  different  calculations 
carried  out  for  SiF2  yield  a  value  of  -152  ±  2  kcal/mol.  The  experimental  heat  of  formation  of  SiF 
(-5  kcal/mol)  also  does  not  agree  well  with  theoretical  predictions.  For  this  system  a  single 
reference  configuration  should  be  adequate  for  correlation  energy  treatments,  which  all  appear  to 
predict  a  heat  of  foimation  of  -14  ±  1  kcal/mol.  These  calculations  suggest  that  the  experimental 
situation  should  also  be  reexamined  for  SiF  and  that  a  more  complete  theoretical  calculation  of  the 
SiF  bond  dissociation  energy  should  be  carried  out 

Finally,  we  report  in  Table  5  our  calculated  electron  affinities,  based  on  G1  and  G2  energies. 
Only  SiF3"  is  known  experimentally  with  an  indicated  electron  affinity  <  2.95  eV  (Reference  108). 
We  find  that  both  SiF  (^fl]  and  SiF2  [^A^]  should  exhibit  positive  electron  affinities,  but  of  smaller 
magnitude  than  the  corresponding  SiH  and  SiH2  molecules.  SiF"  (E^  =  0.83  eV),  however,  should 
easily  be  accessible  experimentally. 

3.1.4  Thermochemistry  of  SiF  and  SiF+ 

As  described  above,  we  find  unsatisfactory  agreement  between  our  calculated  heats  of 
formation  for  SiF  and  SiF2  and  experimental  values.  In  particular,  the  experimental  heat  of 
formation  of  SiF  (-5  kcal/mol)  is  in  poor  agreement  with  all  theoretical  calculations  performed  to 
date.  In  an  effort  to  resolve  this  discrepancy,  we  have  undertaken  large  scale  correlated  energy 
calculations  of  SiF  and  SiF"^  at  the  QCISD(T)  level  of  theory  using  increasing  accurate  basis  sets. 
The  calculations  were  carried  out  at  the  experimental  bond  lengths:  SiF  (rg  =  1.601 1  A)  (Reference 
98);  SiF+  (re  =  1.5265  A)  (References  1 10  and  1 1 1). 

The  smallest  basis  set  examined  was  6-31 1+G(3d0-  The  G1  and  G2  methods  are  based  on 
this  basis  although  additivity  approximations  are  employed  to  achieve  a  practical  computational 
level.  For  this  simple  diatomic  system  we  can  afford  to  choose  much  larger  basis  sets.  We 
examined  in  increasing  size  extent:  2)  aug-cc-pVTZ,  the  triple-zeta-nliffuse  correlation  consistent 
set  proposed  by  Dunning  (Reference  1 12);  3)  the  contracted  9s6p  +  diffuse  +  polarization  basis  of 
Handy  (Reference  28);  4)  basis  set  3  augmented  by  additional  diffuse  and  polarization  functions; 
and  5)  aug-cc-pVTZ,  the  quatro-zeta+diffuse  correlation  consistent  basis  of  Dunning  (Reference 
1 12).  We  were  not  able  to  include  the  g-functions  specified  in  basis  5  and  the  effects  of  their 
neglect  are  uncertain  but  probably  small. 
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As  a  test  of  the  accuracy  of  our  basis  sets  and  procedure  we  first  examined  the  ionizadon 
potential  of  SI  and  the  electronic  affinity  of  F.  The  results  are  shown  in  Table  6  where  we  see  that 
the  Dunning  sets  are  remarkably  consistent  showing  systematic  improvement  going  from  triple  to 
quatro-zeta  form.  Basis  set  4,  which  is  constructed  from  Handy's  basis  set  3  by  geometrically 
splitting  both  the  diffuse  set  and  adding  an  additional  df  polarization  group,  is  our  best 
representation  for  F.  We  conclude  that  basis  sets  4  and  5  should  both  represent  maximum  errors 
of±0.05eV. 

Our  calculated  total  energies  were  corrected  for  zero-point  contributions  using  experimental 
spectroscopic  data.  We  find  Dq®  for  SiF  =  5.95  ±  0.05  eV,  for  SiF"*"  =  6.72  ±  0.05  eV  and  Ip 
for  SiF  =  7.35  ±  0.05  eV.  These  data  translate  into  a  heat  of  formation  for  SiF  of  (0  K)  = 

-  12.2  ±  1.0  kcal/tnol,  a  value  in  good  agreement  with  the  G2  and  previous  theoretical  results  and 
in  disagreement  with  the  values  recommended  by  Walsh  (Reference  85).  Further,  we  find  good 
agreement  with  the  most  recent  experimental  ionization  potential  for  SiF  (7.5  ±  0.4  eV)  reported  by 
Hayes,  et  al.  (Reference  113).  Our  calculated  heat  of  formation  for  SiF+  is  AHf^  (0  K)  =  -h158.0 
kcal/mol,  a  value  in  good  agreement  with  previous  theoretical  studies  (Reference  1 14).  The  heat  of 
formation  of  SiF+  reported  by  Weber  and  Armentrout  (Reference  1 15)  is  ~10  kcal/mol  too  large 
suggesting  a  small  barrier  in  the  Si"*"  +  BF3  — »  SiF"^  +  BF2  reaction  or  possible  errors  in  the 
thermochemistry  of  BF2  or  BF3. 

Parallel  studies  at  the  QCISD(T)  level  of  theory  are  currently  in  progress  using  the  basis  sets 
outlined  above.  To  date,  they  support  a  heat  of  formation  for  SiF2  of  AHf®  (298  K)  =  -152 
kcal/mol,  a  value  some  1 1  kcal/mol  lower  than  the  accepted  experimental  value  (Reference  98). 
Both  SiF  and  SiF2  should  be  reexamined  experimentally  to  verify  our  predicted  thermochemistry. 

3.2  Ion-Molecule  Clustering  Reactions 

In  the  first  part  of  this  research  effort,  an  analysis  of  the  Si"*"  +  SiHj  clustering  reaction  was 
initiated.  Mandich,  et  al.  (Reference  1 1)  have  experimentally  studied  the  silicon  cluster  growth 
reactions  starting  with  Si+  +  SiH4.  They  find  that  the  addition  sequence  stops  at  the  Si5Hio'*' 
structure.  In  a  parallel  theoretical  study,  Raghavachari  (Reference  12)  reponed  ab  initio 
calculations  on  these  silicon  growth  reactions.  The  structures  were  optimized  at  the  HF/6-31G* 
level  of  theory.  The  initial  ion-molecule  complex  that  is  formed  in  the  Si"*"  SiH4  reaction  is 
reported  to  be  a  near  linear  attack  of  Si"*"  along  an  Si-H  bond  as  shown  in  Figure  1 .  The  bridge 
structure,  also  shown  in  Figure  1,  is  reported  to  be  -  1.5  kcal/mol  higher  in  energy.  The  barrier 
between  these  two  structures  has  not  been  calculated  but  is  estimated  to  be  low  owing  to  the 
extended  Si+-SiH4  bond.  Either  of  these  structures  is  stabilized  by  the  polarization  forces  of  Si"^ 
on  SiHt.  These  result  in  a  long-range  [-a/R4],  attractive  force  balanced  by  exchange  repulsion. 
The  relative  stability  of  the  two  structures  shown  in  Figure  1  is,  therefore,  critically  dependent  on 
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Table  6.  SiF  and  SiF-*^  Total  Electronic  Energies  and  Calculated  Thermochemistry 
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ZPE(eV):  SiF  (0.0530);  SiF+ (0.0650);  [  ]  experimental;  (  )  including  counterpoise  correction 


the  value  of  the  polarizability  of  SiH4.  Structure  2  should  be  favored  for  a  larger  value  of  the 
polarizability.  It,  therefore,  becomes  important  to  carefully  examine  the  polarizability  of  SiHj. 

Haaland  (Reference  10)  has  pointed  out  that  the  calculated  polarizability  of  SiH4  is  too  small 
(3.60  A})  even  at  the  MP2/6-31++G**  level  of  theory.  The  best  experimental  estimate  for  a  is 
4.62  A3,  according  to  Haaland.  We  have  computed  the  polarizability  using  several  basis  sets  at  the 
MP2  level  of  theory  available  in  the  CADPAC  program.  The  results  are  summarized  in  Table  7 
which  illustrates  that  the  use  of  diffuse  basis  functions  on  both  Si  and  H  is  required  for  a 
quantitative  estimate  of  the  polarizability.  Both  the  MP2/6-31G**(sp,d,Sn)  and  MP2/6- 
31  lG**(sp,d,Sn)  calculations  yield  a  calculated  value  for  the  polarizability  within  90  percent  of 
Haaland's  estimate.  The  calculated  value  at  the  MP2/6-31G*  level  of  theory  (a  =  3.21  A3)  is 
significantly  smaller  and,  therefore,  underestimated  the  long  range  attractive  force  in  the  Si+  + 

SiH4  reaction.  This  adversely  affects  the  relative  energetics  of  structures  1  and  2  shown  in 
Figure  1. 

To  assess  the  role  of  the  polarization  potential  in  determining  the  relative  stability  of  these 
loose  charge  transfer  complexes,  a  series  of  calculations  were  undertaken.  Using  the  6-31G**(sp, 
d,Sn)  basis,  optimized  UHF  calculations  were  carried  out  for  both  the  Cg  Si-H  bond  structure  and 
the  C2v  bridge  structure.  The  results  are  given  in  Table  8  which  again  indicate  that  the  calculated 
energies  are  too  close  at  the  UHF  level  to  choose  between  these  two  clustering  mechanisms. 
However,  all  calculations  at  a  correlated  level  indicate  that  the  C2v  bridge  structure  is  lower  in 
energy.  These  new  data  shed  doubt  on  the  suggested  Si2H4'''  equilibrium  structure  (Reference  12). 
Insertion  of  an  Si  atom  into  the  long  Si-H  bond  in  the  Cs  ion-molecule  structure  yields  the  (Cg) 
equilibrium  structure  shown  at  the  top  of  Figure  2.  However,  since  the  C2v  bridge  structure 
shown  at  the  bottom  of  Figure  2  is  lower  in  energy  at  long  Si-H2  bond  lengths,  this  may  actually 
represent  the  more  stable  equilibrium  conformation. 

The  study  of  reaction  pathways  for  the  corresponding  negative  ion  clustering: 

StnH2n-2' Sin=iH2n‘+H2  (n=l,4)  (104) 

was  begun  by  considering  the  initial  Si'  +  SiH4  reaction  surface.  The  interaction  of  Si"  with  SiH4 
can  proceed  along  two  pathways  leading  initially  to  the  ion-molecule  complexes  shown  in  Figure  3 
Studies  of  the  SI  structure  were  carried  out  at  the  HF/6-31+G**  level  of  theory  with  the  Si-Si 
bond  treated  as  a  reaction  parameter.  Our  results  are  given  in  Table  9.  We  find  a  weak  polarization 
attraction  at  a  Si-Si  bond  length  of  ~  4.5  A.  However,  there  is  no  significant  distortion  of  the  Si-H 
bond  in  silane  along  this  reaction  coordinate.  This  result  is  in  contrast  to  the  Si"*"  +  SiH4  reaction  in 
Cs  symmetry  where  there  is  clearly  a  lengthening  of  the  Si-H  bond  to  -  1.6  A  along  the  reaction 
coordinate.  Further  shortening  of  the  Si-Si  bond  gives  rise  to  a  strongly  repulsive  interaction. 

The  transition  state  for  insertion  of  Si"  into  SiH4  thus  appears  to  be  highly  energetically  in  Cj 
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Table  7.  Calculated  MP2  Polarizabmties  of  SiH4 


Method 

R(Si-H)(A) 

E(hartrees) 

a(A3) 

MP2/6-31G* 

t.4829 

-291.3168497 

3.20 

MP2/6-31G** 

1.4726 

-291.3498595 

3.17 

MP2/6-31++G** 

1.4719 

-291.3531119 

3.59 

MP2/6-31G**(sp,d,Sn) 

1.4730 

-291.3619018 

4.12 

MP2/6-311++G** 

1.4737 

-291.4967731 

3.74 

MP2/6-311G**(sp,d,s„) 

1.4759 

-291.5030978 

4.21 

EXP. 

4.62 

Ojp  =  0,025795;  =  0.15;  =  0.036 
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Table  8.  Energetics  of  the  Si'*'  +  SiH4  Clustering  Reaction 


Method 

Geometry 

Energy 

h 

e 

Cj  Structure 

UHF/6-31G** 

1.460 

1.601 

1.802 

180. 

-  579.80980 

UHF/6-31G**(sp,d,Sn) 

1.461 

1.600 

1.784 

180. 

-579.81459 

MP2/6-31G** 

MP4/6-31G**/AJHF 

1.460 

1.588 

1.753 

180. 

-579.98130 

-  580.02340 

C2v  Structure 

UHF/6-31G** 

1.457 

1.530 

1.971 

-  579.80650 

UHF/6-31G**(sp,d,Sn) 

1.457 

1.534 

1.975 

-579.81249 

MP2/6-31G** 

MP4/6-31G**/AJHF 

1.458 

1.537 

2.148 

-  579.98620 

-  580.02720 

Bond  lengths  in  A,  bond  angles  in  degrees,  energy  in  hartiees. 
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Bond  Structure 


Bridge  Structure 


Figure!.  Si'*'-SiH4  Ion-Molecule  Complexes. 
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Figure  3.  Si’'>SiH4  lon-MoIecule  Complexes. 
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Table  9.  Potential  Energy  for  Si*  *■  SiH4  lon-Mcrfectile  Surface 
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symmetry.  Calculations  were  also  carried  out  in  C2v  symmetry  with  similar  conclusions,  a  weak 
polarization  attraction  followed  by  a  strongly  repulsive  region. 

The  Si~  +  SiHi  reaction  was  then  examined  in  C}  symmetry.  The  lowest  reaction  pathway  is 
illustrated  in  Figure  4  where  we  compare  the  transition  stmcture  with  the  stable  silylsilylene 
structure.  In  contrast  to  the  Si+  +  SiH4  reaction,  where  the  transition  state  for  the  insertion  of  Si^ 
is  exothermic  relative  to  Si"*"  +  SiH4,  we  find  a  positive  reaction  barrier  for  the  corresponding 
negative  ion  reaction.  To  explore  this  reaction  step  further,  we  next  carried  out  calculations  of  the 
stationary  points  of  Si2H2"  and  Si2H4"  at  the  UHF/6-31+G**  level  of  theory.  Six  stationary 
structures  for  Si2H2"  and  four  stationary  structures  for  Si2H4"  were  located.  TTie  calculated 
energetics  and  description  of  these  structures  are  given  in  Table  10.  The  minimum  energy  ground 
state  (all  stable  vibrational  frequencies)  of  Si2H2"  corresponds  to  a  nonplanar  ttans  structure  in  C2 
symmetry.  A  planar  trans  [C2h]  structure,  lying  only  0.14  eV  above  the  ground  state,  is  also 
found  to  be  stable.  However,  the  geometry  of  this  structure  appears  to  be  dependent  cm  the  chosen 
basis  set  and  this  state  may  be  an  artifact  of  the  6-314G**  basis.  Two  other  stable  Si2H2" 
structures  were  Icxrated:  a  C2v  bridge  structure  and  a  C2v  planar  cis  structure.  Two  other  stable 
stationary  points  are  identified  as  saddle  regions  on  the  potential  surface.  The  optimized  structures 
for  Si2H2"  are  illustrated  in  Figure  5.  Five  stationary  points  for  Si2H4"  were  located.  The 
calculated  energetics  and  description  of  these  structures  are  given  in  Table  10.  Both  the  nonplanar 
trans  ethyelene-like  structure  [C^hl  l^sl  silylsilylene  stmcture  exhibited  all  stable 

vibrational  frequencies.  All  other  conformations  corresponded  to  transiticm  or  degenerate  saddle 
regions  on  the  potential  energy  surface.  The  optimized  structures  for  Si2H4"  are  illustrated  in 
Figure  6. 

An  examination  of  the  stmcture  illustrated  in  Figures  4  and  6  leads  to  the  following  reaction 
step  sequence: 


Si-  +  SiH4  ^  SiH3  -  SiH*  (TS) 

AE  ~  +  10  kcal/mol 

(105) 

SiHj-SiH-  (TS)  SiH2^SiH2- 

AE  =  +  57  kcaVmol 

(106) 

SiH3--SiH-  (TS)  ->  SiH3-SiH- 

AE  =  -  29  kcal/mol 

(107) 

SiH3-SiH-  ->  SiH-SiH-  +  H2 

AE  =  +  41  kcal/mol 

(108) 

In  contrast  to  the  positive  ion  reaction  sequence,  the  addition  of  Si“  to  SiH4  to  form  Si2H2*  + 
H2  is  an  overall  endothermic  process.  The  major  anions  to  be  expected  under  most  discharge 
conditions  are  SiH3",  SiH2'  and  SiH",  formed  from  dissociative  attachment  of  e  +  SiH4  and 
possibly  Si2H4'  which  may  exhibit  a  long  lifetime  as  a  disilene  stmcture.  Independent  studies  by 
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transitim  state  of  Si~-fSi2H4 


H 


silysilylene  structure  of  Si2H4 


Figure  4.  Si'  +  SiH4  Reaction  Pathway  in  [C|]  Symmetry. 
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Table  10.  Calculated  Stability  and  Energetics  for  Si2Ho'  Conformations 
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Figure  5.  Optimized  Structures  of  Si2H2'‘. 
70 


Structure 


H 

106.46®  X  114.03° 


^  H 


Sl-i - Si 

2.292 


Symmetry 
Stgte  Group 


^^2  C2v 


Figure  6.  Optimized  Structures  of  Si2l 
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Tada  (Reference  116)  have  shown  that  the  disilane  radical  anion,  Si2H5~,  is  not  thermodynamically 
bound  relative  to  Si2H5  -t-  e  but  may  exist  for  several  ^sec  as  a  metastable  kinetically  bound  anion. 
It  thus  appears  that  negative  ions  have  little  potential  for  addition  reactions  of  the  type  illustrated  in 
Equation  (104). 

3.3  e  +  SiH4  Collisional  Dissociative  Attachment 

Ab  initio  calculations  of  the  dissociative  reaction  surfaces  for  negative  ion  formation  in  silane 
were  undertaken.  In  Table  1 1,  we  show  the  possible  product  branching  for  e  +  SiH4.  Two 
possible  kinetic  routes  must  be  considered  for  conditions  that  exist  in  the  cathode  sheath  region  of  a 
silane  glow  discharge.  The  first  assumes  that  an  unstable  resonance  state  of  SiH4"  is  formed  by 
primary  electron  capture.  Subsequently,  this  unstable  negative  ion  can  decay  into  the 
thermodynamically  accessible  product  channels  which  correspond  to  dissociation  products.  Table 
1 1  illustrates  that  reactions  (2)  -  (S)  fit  this  criteria  but  that  (6)  thru  (8)  do  not  Dissociation  of 
SiH4"  should  exhibit  thresholds  of  1.25, 2.58, 2.93  and  3.23  eV  for  formation  of  SiH2',  SiH3", 
Si~  and  H~,  respectively,  provided  that  there  are  no  barriers  for  dissociation  in  the  exit  channels. 
Although  taken  with  limited  energy  resolution,  the  data  of  Potzinger  and  Lampe  (Reference  117) 
and  Ebinghaus,  Krauss,  Mtiller-Duysing,  and  Neuert  (Reference  118)  appear  to  be  in  agreement 
with  these  predictions  for  formation  of  SiH2’  and  SiH3",  although  a  unique  reaction  pathway 
cannot  be  determined  from  these  experiments.  Srivastava  (Reference  1 19)  has  also  observed 
resonances  at  ~  2.5  eV  which  he  attributes  to  SiH2’  and  SiH3"  formation.  The  process  of  negative 
ion  formation  following  decay  of  a  resonance  state  of  SiHj"  was  first  described  by  Haaiand 
(Reference  120). 

In  order  to  analyze  the  possible  dissociation  products  for  low  energy  collisions,  the  reaction 
surfaces  governing  the  dissociative  attachment  (DA)  of  electrons  to  SiHa  were  analyzed  using  th ' 
intrinsic  reaction  coordinate  procedure  of  Schmidt,  Gordon,  and  Dupuis  (Reference  121).  In 
Figure  7,  we  show  the  energetics  of  the  reaction  coordinate  leading  to  dissociation  of  SiH4  to  SiH3 
+  H.  At  the  HF/6-31++G**  level  of  theory,  we  find  a  bond  dissociation  energy  of  3.27  eV.  At 
the  MP2/6-31++G**  level,  a  bond  energy  of  3.85  eV,  very  close  to  the  experimental  value,  is 
predicted.  Shown  on  Figure  7  are  data  for  the  SiH3'  +  H  and  SiH3  +  H"  reaction  surfaces.  We 
find  an  avoided  crossing  of  these  two  surfaces,  very  close  to  the  intemuclear  separation  where  the 
ion  curves  cross  the  neutral  SiH3  +  H  +  e  surface,  but  probably  in  the  region  where  the  ion 
potentials  are  complex.  Thus,  the  excited  state  ion  channel,  leading  to  SiH3  +  H'  may  depq)ulate 
through  this  avoided  crossing  leaving  SiH3~  +  H  as  the  only  significant  products  of  the  DA 
reaction.  Experimental  data  by  Haaiand  (Reference  120)  and  by  Potzinger  and  Lampe  (Reference 
1 17)  indicate  that  H*  is  probably  not  formed  by  direct  DA  of  e  +  SiHa.  Srivastava  (Reference  1 19) 
reports  a  small  cross  section  (~2xl0''9  cm^)  for  formation  of  H"  by  direct  DA  of  e  +  SiH4  at  a 
collision  energy  of  ~9.0  eV.  The  coupling  of  the  SiH3"  +  H  and  SiH3  +  H"  reaction  surfaces 
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Table  11.  Thermodynamics  of  Negative  Ion  Formation  in  Silane 


e  +  SiH4[^A2] 


Reaction  AH  (eV) 


SiH4-[2T2] 

(+3.98) 

(1) 

-> 

SiH2-[2Bi] 

+  H2[»Ig+] 

+1.35 

(2) 

-> 

SiH3-[lAi] 

+  H[2S] 

+2.58 

(3) 

-> 

Si-[4S] 

+  2H2[ii:g+] 

+2.93 

(4) 

SiH3[2Ai] 

+  H-[*S] 

+3.23 

(5) 

SiH-[3E-] 

+  H2[i2:g+1 

+  H[2S] 

+4.52 

(6) 

SiH[2n] 

+  H2[lXg+] 

+  H-[1S] 

+5.04 

(7) 

— > 

Si[3p] 

+  H2[l2g+] 

+  H[2S] 

+  H-[lS] 

+8.08 

(8) 
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should  lead  mainly  to  SiH3'  +  H  products  since  the  exit  reaction  channels  are  strongly  mixed.  One 
explanation  is  that  Srivastava  is  probing  the  resonant  excited  states  of  SiH4"  resulting  from  s-»p 
promotion  on  the  silicon  atom.  The  statistical  breakup  of  these  states  into  reaction  product 
channels  could  lead  to  H‘  formation,  albeit  by  the  indirect  attachment  mechanism  described  below. 

To  facilitate  analysis  of  the  asymptotic  behavior,  we  have  calculated  the  polarizability  of  SiH3 
at  the  HF/6-31++G**  level  and  find  a  =  3.95  A^.  The  corresponding  value  of  the  polarizability  of 
SiH4  at  this  level  of  theory  is  3.25  A^.  The  leading  term  of  the  long-range  potential  of  SiH3  +  H’ 
and  SiH3'  +  H  can  be  calculated  according  to  the  form:  V(R-^)  =  -  a/2R^.  These  asymptotic 
curves  are  shown  on  Figure  7  which  clearly  indicate  that  an  avoided  crossing  between  these  two 
channels  must  occur. 

Preliminary  studies  of  dissociative  attachment  (DA)  of  electrons  to  silane  were  undertaken 
through  analysis  of  the  possible  product  channels  and  quantum  coupling  with  the  reactants.  The 
lowest  adiabatic  attachment  channel  correlates  to  SiH2'[^Bi]  +  H2[XlZg''']  at  1.35  eV.  The  SiH^' 
resonance  state  corresponding  to  this  channel  appears  to  be  of  the  Feshbach  type  with  a  small 
capture  width  from  reactants.  In  addition,  as  studies  described  below  will  indicate,  there  is  a  large 
barrier  on  the  reaction  surface  leading  to  SiH2"  +  H2  products.  The  second  DA  channel 
corresponds  to  the  formation  •.>  SiH  '!  ’Aj]  +  H[2S]  at  2.58  eV.  This  channel  corresponds  to  a 
shape  resonance  but  also  appear,  to ,  ve  a  small  DA  cross  section  for  vibrationally  cold  SiH4. 
Srivastava  (Reference  119)  finds  a  resonance  in  the  SiH3’  +  H  channel  at  -  2.7  eV  but  this 
experimental  energy  appears  to  be  the  same  for  formation  of  Si‘,  SiH*,  SiH2*  and  SiH3*,  which 
should  exhibit  quite  different  threshold  behavior.  Recently,  Wan,  Moore,  and  Tossell  (Reference 
122)  have  reported  electron  scattering  cross  sections  for  silane  in  the  range  0.2— >12.0  eV.  They 
observe  a  broad  maximum  in  the  total  electron  scattering  cross  section  at  ~  2.9  eV  which  they 
attribute  to  a  virtual  transition  to  the  lowest  unoccupied  o*(t2)  level  in  SiH4.  This  feature  has  a 
large  cross  section  (~  50  A^)  and  may  be  related  to  the  features  seen  by  Srivastava  in  dissociative 
attachment  (DA)  at  ~  2.5  eV.  Based  on  these  known  low-energy  quantum  processes,  our 
conclusion  is  that  low-energy  direct  electron  dissociative  attachment  processes  should  have  low 
probability. 

The  second  ki.ietic  route  for  DA  of  e  +  SiH4  was  examined  by  looking  at  the  excited  state 
surfaces  of  SiHj*.  The  DA  mechanism  for  this  route  may  be  the  initial  capture  of  an  electron  into  a 
high  lying  SiHj*  state  with  subsequent  statistical  branching  into  the  allowed  product  channels.  For 
a  mechanism  of  this  type,  the  product  distribution  is  uncorrelated  with  the  initial  dynamics  of  the 
reaction  channel.  The  data  of  Srivastava  suggest  that  a  stochastic  formation  of  product  states,  all 
formed  from  the  same  state(s)  of  SiH4*  lying  at  8-9  eV  is  occurring. 

In  order  to  identify  the  excitation  channels  for  e  -»■  SiH4,  extensive  configuration  interaction 
(Cl)  studies  of  SiH4"  were  next  carried  out  to  identify  the  excited  state  reaction  surfaces  that  play  a 
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role  in  dissociative  attachment  of  electrons.  In  Tables  12-14,  we  show  the  excitation  spectrum  for 
SiH4'  calculated  at  the  equilibrium  geometry  of  SiH4  (R  =  1.477  A).  These  spectra  consist  of 
low-lying  resonance  stmctures  below  ~  8  eV  and  a  large  manifold  of  excited  states  corresponding 
to  electron  capture  by  a  Rydberg  state  of  SiH4-  These  latter  states  are  most  likely  responsible  for 
the  high  energy  (>7  eV)  dissociative  recombination  of  e+SiH4,  outlined  in  our  previous  work. 
The  spectra  are  given  in  Tables  12-14  for  three  levels  of  basis  sets.  The  improvement  in  basis 
from  Table  12  to  Table  14  should  mainly  affect  the  location  of  the  Rydberg-like  states  since  the 
more  flexible  triple  split-zeta  basis  primarily  improves  description  of  the  Si  electrons.  The 
improvement  in  basis  from  Table  13  to  Table  14  consists  of  added  diffuse  functions  to  the 
hydrogen  atoms.  This  permits  a  better  description  of  large  diffuse  charge  clouds  characteristic  of 
shape  resonances.  We  show  the  comparative  spectra  in  Figure  8  which  illustrates  the  stability  of 
the  states  lying  above  ~  8  eV  to  basis  set  improvement. 

At  the  equilibrium  geometry  of  SiH4,  the  low-lying  excited  state  surfaces  are  triply 
degenerate  since  the  full  point  group  with  this  symmetry  is  Tj.  In  the  C3V  subgroup,  the  low-lying 
states  of  SiH4~  are  all  of  symmetry,  which  is  representative  of  a  loosely  bound  electron, 
symmetrically  distributed  about  SiH4.  These  are  pseudoeigenstates  of  the  system,  strongly 
dependent  on  the  chosen  basis  set,  and  are  not  physically  meaningful.  However,  the  third 
excitation  (at  9.10  eV)  is  characterized  by  a  compact  charge  distribution  resulting  from  s-»p 
promotion  on  the  Si  atom.  We  identify  this  state  as  2A2  in  C3V  symmetry  (^A"  in  Cg)  and  believe 
that  it  represents  the  first  identification  of  the  resonant  excited  states  in  this  system.  Our  location  of 
this  state  at  9.1  eV  is  in  good  agreement  with  previous  experimental  data  placing  the  important 
dissociative  attachment  states  in  the  7-10  eV  range  (Reference  120).  To  give  some  feel  for  the 
magnitude  of  these  calculations,  our  configuration  interaction  (Cl)  expansion  is  composed  of  all 
single,  double,  triple,  and  quartic  excitations  in  an  active  orbital  space  of  the  lowest  16  natural 
orbitals.  This  results  in  a  Cl  comprising  some  28,000  configurations.  The  identification  of  the  9.1 
eV  state  as  ^A2  suggests  that  dissociative  attachment  to  this  state  should  not  be  strong  owing  to 
symmeoy  restrictions  which  weakly  couple  this  stare  to  the  ^Aj  channel  corresponding  to  e  + 

SiH4.  Based  on  this  initial  study,  a  more  thorough  examination  of  excited  state  surfaces  of  SiH4* 
was  undertaken.  Examination  of  the  eigenvalues  \|r4  -4  Yio  indicated  a  high  density  of  states  in 
the  7.1  to  10.5  eV  region.  The  next  higher  eigenstates  appear  beyond  14  eV  and  must  represent 
autoionizing  states  of  SiH4*. 

The  ^Aj  states  lying  in  the  7.1  - 10.5  eV  region  can  correlate  to  the  totally  symmetric 
products  SiH4"  SiH3"  +  H.  The  states  correlate  to  a  mixture  of  SiH3"  +  H  and  SiH2"  +  H2, 
the  latter  formed  as  the  products  of  a  Rydberg  state  with  SiHa"*"  core-like  structure.  The  direct 
fomiation  of  SiH'  rom  DA  of  e  +  SiH4  is  much  less  probable.  An  ab  initio  calculation  of  the 
cross-section  for  DA  in  e  +  SiH4  is  extremely  difficult  since  knowledge  of  both  the  excited  states 
of  SiH4*  and  the  width  for  initial  electron  capture  are  required.  The  cross  section  for  electron 
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Table  12.  Excitation  Spectra  of  SiH4“.  (CIS/6-31+G**  Level  of  Theory) 


Excited  State 

Tc(eV) 

X  (nm) 

f-number  * 

1 

0.4739 

2616.13 

0.1373 

2 

0.4750 

2610.25 

0.1376 

3 

0.4760 

2604.58 

0.1380 

4 

3.1698 

391.14 

0.0178 

5 

3.1874 

388.99 

0.0174 

6 

3.2049 

386.86 

0.0170 

7 

7.0946 

174.76 

0.0000 

8 

9.0813 

136.53 

0.0836 

9 

9.0943 

136.33 

0.0838 

10 

9.1071 

136.14 

0.0837 

11 

9.2068 

134.66 

0.0002 

12 

9.3586 

132.48 

0.0004 

13 

9.3684 

132.34 

0.0004 

14 

9.3793 

132.19 

0.0003 

15 

9.6702 

128.21 

0.0000 

16 

9.6883 

127.97 

0.0000 

17 

9.7034 

127.77 

0.0000 

18 

9.7174 

127.59 

0.0000 

19 

9.7261 

127.47 

0.0002 

20 

9.7372 

127.33 

0.0000 

21 

9.9425 

124.70 

0.0122 

22 

9.9525 

124.58 

0.0121 

23 

9.9608 

124.47 

0.0123 

24 

10.2594 

120.85 

0.0000 

25 

10.2869 

120.52 

0.0000 

26 

10.2988 

120.39 

0.0000 

27 

10.3109 

120.25 

0.0000 

28 

10.3233 

120.10 

0.0000 

29 

10.3988 

119.23 

0.0000 

30 

10.7386 

115.46 

0.0000 

*  The  f-number  corresponding  to  the  transition  from  the  ground  state  SiHt"  structure  to 
the  corresponding  excited  state. 
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Table  13.  Excitation  Spectra  of  SiH4'.  (CIS/6-3114-G**  Level  of  Theory) 


Excited  State 

Te(eV) 

X  (nm) 

f-number  ♦ 

1 

0.2284 

5427.21 

0.0680 

2 

0.2299 

5393.27 

0.0683 

3 

0.2315 

5356.24 

0.0685 

4 

2.3163 

535.26 

0.0096 

5 

2.3281 

532.56 

0.0097 

6 

2.3397 

529.91 

0.0099 

7 

5.0138 

247.29 

0.0001 

8 

8.0567 

153.89 

0.0028 

9 

8.0632 

153.76 

0.0030 

10 

8.0695 

153.64 

0.0031 

11 

9.0059 

137.67 

0.0428 

12 

9.0105 

137.60 

0.0390 

13 

9.0125 

137.57 

0.0353 

14 

9.1453 

135.57 

0.0000 

15 

9.2237 

134.42 

0.0118 

16 

9.2396 

134.19 

0.0149 

17 

9.2575 

133.93 

0.0192 

18 

9.3690 

132.33 

0.0001 

19 

9.3829 

132.14 

0.0002 

20 

9.3866 

132.09 

0.0003 

21 

9.3969 

131.94 

0.0005 

22 

9.4049 

131.83 

0.0002 

23 

9.4073 

131.79 

0.0002 

24 

9.7755 

126.83 

0.0528 

25 

9.7920 

126.62 

0.0512 

26 

9.8101 

126.38 

0.0501 

27 

9.9629 

124.44 

0.0001 

28 

9.9713 

124.34 

0.0001 

29 

9.9801 

124.23 

0.0001 

30 

9.9863 

124.15 

0.0000 

31 

9.9947 

124.05 

0.0001 

32 

10.3587 

119.69 

0.0000 

33 

10.5192 

117.86 

0.0000 

34 

10.5277 

117.77 

0.0000 

35 

10.5361 

117.67 

0.0000 

36 

10.5592 

117.42 

0.0000 

37 

10.5767 

117.22 

0.0000 

38 

10.8243 

114.54 

0.0254 

39 

10.8395 

114.38 

0.0248 

40 

10.8548 

114.22 

0.0242 

*  The  f-nutnber  corresponding  to  the  transition  from  the  ground  state  SiH4"  structure  to 
the  corresponding  excited  state. 
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Table  14.  Excitation  Spectra  of  SiH4~.  (ClS/6-311-»-»-G**  Level  of  Theory) 


Excited  State 

Te(eV) 

X(nm) 

f-number  * 

1 

0.6718 

1845.6 

0.1932 

2 

0.6741 

1839.15 

0.1957 

3 

0.6764 

1832.88 

0.1981 

4 

1.2496 

992.22 

0.1103 

5 

1.2507 

991.31 

0.1109 

6 

1.2519 

990.40 

0.1116 

7 

3.4008 

364.57 

0.0000 

8 

4.0376 

307.07 

0.0055 

9 

4.0515 

306.01 

0.0054 

10 

4.0653 

304.98 

0.0052 

11 

5.5362 

223.95 

0.0000 

12 

8.7745 

141.30 

0.1427 

13 

8.7780 

141.24 

0.1425 

14 

8.7813 

141.19 

0.1423 

15 

8.9660 

138.28 

0.0006 

16 

8.9677 

138.26 

0.0007 

17 

8.9701 

138.22 

0.0009 

18 

9.1449 

135.58 

0.0001 

19 

9.3693 

132.33 

0.0014 

20 

9.4057 

131.82 

0.0282 

21 

9.4115 

131.74 

0.0293 

22 

9.4206 

131.61 

0.0291 

23 

9.4409 

131.33 

0.0003 

24 

9.4507 

131.19 

0.0007 

25 

9.4540 

131.14 

0.0000 

26 

9.4637 

131.01 

0.0001 

27 

9.4745 

130.86 

0.0001 

28 

9.7793 

126.78 

0.1255 

29 

9.7871 

126.68 

0.1259 

30 

9.7957 

126.57 

0.1258 

31 

9.8841 

125.44 

0.0001 

32 

9.8955 

125.29 

0.0004 

33 

9.9060 

125.16 

0.0000 

34 

9.9538 

124.56 

0.0000 

35 

9.9737 

124.31 

0.0000 

36 

10.4508 

118.64 

0.0000 

37 

10.4584 

118.55 

0.0000 

38 

10.4660 

118.46 

0.0000 

39 

10.4826 

118.28 

0.0000 

40 

10.4968 

118.12 

0.0000 

*  The  f-number  corresponding  to  the  transition  from  the  ground  state  SiH^*  structure  to 
the  corresponding  excited  state. 
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Excitation  Level  (eV) 


6-3UG**  6-31  UG‘*  6-31 1++G** 


Figure  8.  Vertical  Excitation  Spectra  of  SiH4'. 
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capture  in  DA  is  determined  by  the  asymptotic  form  of  the  nuclear  motion  wavefunction,  ^(R). 
For  electrons  with  initial  energy  we  have  for  the  cross  section; 


ObA 


4ff^mK 

Mki 


lim  1^(R)|2 
R— >oo 


(109) 


Here  ^(R)  represents  the  nuclear  modon  wavefunction  along  the  minimum  energy  pathway. 
Treating  electron  capture  as  a  direct  process  and  neglecting  interference  with  possible  electron 
vibrational  excitation  processes,  we  can  evaluate  Equation  (109),  provided  some  estimate  of  the 
capture  width  for  electron  attachment  can  be  made  for  the  region  of  SiH4"  to  the  left  of  the  crossing 
points  of  the  product  channels.  Using  Fermi's  golden  mle,  we  can  estimate  the  resonance  width 
as: 


r  =  V(2^ 


(110) 


where  AE  is  the  calculated  vertical  separation  of  the  negative  ion  state  and  that  of  SiH4  and  k  is  an 
average  value  of  the  coupling  matrix  element  for  the  electronic  states  of  the  neutral  and  negative  ion 
species.  We  find  an  average  value  of  F  -  0.01,  a  value  characteristic  of  coupling  to  Rydberg 
states.  We  have  approximately: 


<Tda  =  3.7 X 10“^''  r(Re) <  W  V^f  (111) 

e 

where  <  VblVf  >  is  an  unnormalized  overlap  matrix  element  of  order  unity.  For  a  threshold 
collision  energy  of  9. 1  eV,  we  find: 


Oda  =  0.64  X  10*18  cn:i2 

for  SiHj-  +  H 

(112) 

Oda  =  0.89  X  10’ 18  cm2 

for  S1H2”  +  H2 

(113) 

Although  these  cross  sections  are  of  the  correct  order-of-magnitude,  most  experimental  data 
(Reference  120)  indicate  that  ©^^(SiHj"  +  H)  >  aDA(SiH2"  +  H2).  Further  studies  are  required 
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to  classify  this  discrepancy.  The  problem  may  lie  in  tracing  the  proper  number  of  states  which 
correlate  with  SiH3“  +  H  and  SiH2'  +  H2  products  or  with  our  estimates  of  the  capture  widths. 

The  problem  of  DA  of  e  +  SiH4  at  higher  (>  4  eV)  collisional  energies  can  be  viewed  using  a 
different  mechanism.  In  this  case,  where  the  resonance  state  has  a  long  lifetime,  a  two-step 
mechanism  may  best  describe  the  attachment  process: 

e  +  AB  <=>  [AB**]~  :  <Vab°  Vab**>  (1^^) 

[AB**]-  ->  A-(*,t)  +  B(*,t)  :  <VAB**lHD'VA-(*.t)v|rB(*.t)>  (115) 

Equation  (114)  governs  the  rate  of  capture  into  the  resonance  state  of  the  negative  ion  and  the 
corresponding  rate  of  autoionizadon.  Equation  (115)  governs  the  rate  of  dissociation  of  the 
resonance  state  into  the  possible  product  channels.  This  dissociation  can  occur  along  a  number  of 
exit  channels,  subject  to  the  usual  conservation  laws  of  total  energy  and  angular  momentum.  The 
competition  among  the  various  possible  modes  of  dissociation,  however,  does  not  depend  on  the 
manner  in  which  the  resonance  state  was  initially  formed.  It  follows  that  the  overall  attachment 
probability  can  be  written  as  a  product  of  two  terms: 


W[\)/’j(£r),  VI/'jJj]  —  W^p[\jrj(ER)]  WjJgggy  [V(j] 

(116) 

Wcap[¥i(eR)l  =  271  r(.aplVi(eR)l^ 

(117) 

W(jecayt¥a'l  ~  f  ^tot » 

(118) 

where  Er  is  the  energy  corresponding  to  a  particular  resonance  state  in  the  incident  scattering 
channel,  is  the  energy  width  for  coupling  the  incident  channel  into  the  resonance  state,  and  Fjj 
is  the  partial  width  corresponding  to  the  probability  that  the  resonant  state  will  dissociate  into 
products  a.  Obviously,  the  total  decay  width,  I'tot  “  ^  I'o’  is  the  sum  of  the  partial  widths  for  all 
open  dissociation  channels.  This  mechanism  is  often  described  in  the  literature  of  nuclear  physics, 
where  the  resonance  state  is  referred  to  as  a  "compound  nucleus"  (References  123  and  124).  The 
independent  probabilities  for  the  formation  and  subsequent  decay  of  the  resonance  state  are  thus 
interpreted  as  a  model  whereby  the  resonance  state  has  "forgotten"  its  mode  of  formation  before  it 
decays. 
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This  mcxiel  for  DA  can  also  be  understood  in  terms  of  a  semiclassical  analysis  of  scattering 
theory  (References  125  and  126).  In  such  a  description,  the  cross  section  for  DA  can  be  written  as 
a  product  of  two  factors; 


<yDA(£)  =  <ycap(e)S 


(119) 


where  cj(^(e)  is  the  cross  section  for  the  capture  of  an  electron  to  form  the  resonant  negative  ion 
state.  S  is  the  survival  probability,  which  is  classically  related  to  the  decay  width  of  the  resonant 
state: 


S  =  exp(-|’^*  ^^dR)  (120) 

V 

where  the  region  of  interaction  extends  from  the  classical  capture  radius,  R^,  to  the  effective  radius 
of  stabilization  of  the  ion-atom  or  ion-molecule  products,  for  which  V'^Rj)  =  VO(Rs).  Thus  the 
resonant  state  decays  exponentially  with  n^an  life  l/T. 

If  the  resonance  state  is  long-lived,  a  statistical  breakup  may  occur  with  nearly  equal 
probabilities  (a  priori)  for  all  accessible  product  channels.  This  alternate  kinetic  route  can  be 
derived  from  the  assumption  of  LTE  (or  at  least  quasi-steady  state).  In  this  case,  electron 
attachment  can  occur  not  only  to  SiH4,  but  also  directly  to  the  neutral  silane  fragments  which  result 
from  the  decompositions  of  the  negative  ion  resonance  state.  The  product  reaction  channels  arc 
most  likely  those  correlating  to  repulsive  Rydberg  or  high  lying  valence  states  of  SiH4.  In  Table 
15,  we  illustrate  electron  attachment  to  SiH3,  SiH2  and  SiR  A  wide  spectrum  of  threshold 
energies  is  spanned  by  these  molecules  (0.45  -  2.57  eV)  and  all  possible  negative  ion  species, 
including  H‘,  should  be  energetically  accessible.  The  relative  importance  of  these  several 
processes  is,  of  course,  dependent  on  the  branching  probabilities  for  decomposition  of  the 
resonance  state  and  direct  attachment  process  to  the  various  neutral  silane  fragments.  Knowledge 
of  the  neutral  composition  of  the  plasma  in  the  sheath  region  now  becomes  a  requirement  for 
analyzing  the  relative  importance  of  direct  versus  statistical  breakup  electron  attachment. 

3.4  Negative  Ion-Molecule  Reactions 

3.4.1  Negative  Ion-Molecule  Reaction:  H“  +  SiH4 

Several  studies  have  been  reported  which  examine  the  mechanisms  of  the  H"  +  SiR^ 
reaction.  This  is  believed  to  be  a  classical  second  order  nucleophilic  substitution  reaction  (5^2), 
similar  to  the  H‘  +  CH4  and  F'  +  CH3F  systems  which  involve  substitution  at  the  carbon  center 
with  inversion  of  the  MO  configuration.  Sheldon,  Hayes,  and  Bowie  (Reference  127)  explored 
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Table  15.  Thermodynamics  of  Negative  Ion  Formation  in  SiH/SiH2/SiH3 

Reaction  AH  (eV) 


-  1.441 


e  +  SiH3[2A,] 

SiH3-[’Ail 

+ 

SiH2-[2Bi] 

+ 

-» 

SiH2[»Ai] 

+ 

SiH[2n] 

+ 

e  +  SiH2[lAi] 

SiH2-[2Bi] 

Si-[4S1 

+ 

— > 

SiH-[3r] 

+ 

— ) 

SiH[2n] 

+ 

Sil3p] 

+ 

e  +  SiH[2n] 

-4 

SiH-[3Z-] 

Si-[^S] 

+ 

— > 

Si[3p] 

+ 

H2[»i:g+] 

+  0.53 

H[2S1 

+  1.89 

H-[1S] 

+  2.26 

HI2S]  + 

H-[lS] 

+  5.58 

-  1.12 

H2(J2g+J 

+  0.45 

H[2S] 

+  2.04 

+  2.57 

H[2S]  + 

H-[1S] 

+  5.60 

-  1.28 

HI2S] 

+  1.65 

+  2.28 
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the  fully  relaxed  motion  of  H‘  +  SiH4  leading  to  three  possible  kinds  of  products:  a)  deprotonation 
to  form  H2  +  SiH3",  b)  equatorial  addition  to  form  SiH5",  and  c)  apical  addition  to  form  SiH5". 
They  find  no  barriers  for  direct  attachment  of  H‘  to  silane.  Once  formed,  however,  SiH5"  was 
found  to  exhibit  a  barrier  of  ~  70  kcal.mol  relative  to  SiH3‘  +  H2  decomposition.  These  authors 
fail  to  report  their  calculated  stability  of  SiH5'  relative  to  SiH4  +  H"  reactants.  Earlier  studies  of 
this  reaction  have  been  reported  by  Wilhite  and  Spialter  (Reference  128),  Keil  and  Ahlrichs 
(Reference  129),  and  Baybutt  (Reference  130).  It  was  assumed  that  the  SN2-reaction  of  H"  + 

SiH4  proceeds  in  C3V  symmetry,  finally  forming  SiH5"  as  a  031,  (trigonal  bypyramid)  structure. 
There  was  uncertainty  about  the  stability  of  SiH5"  relative  to  H"  +  Silij  and  some  evidence  of  a 
barrier  along  this  reaction  path.  Keil  and  Ahlrichs  (Reference  129)  and  Baybutt  (Reference  130) 
predict  no  bairier  for  this  reaction  whereas  Wilhite  and  Spialter  (Reference  128)  predict  a  large  (8.6 
kcal/mol)  barrier  along  the  Sn2  reaction  path,  similar  to  that  found  for  the  H'  +  CH4  reaction. 

A  careful  study  of  the  H"  +  SiH4  reaction  shows  that  the  calculated  interaction  energy  is  very 
sensitive  to  the  chosen  basis  set.  In  the  conformational  change  from  H"  +  SiH4  to  SiH5"  to  SiH3" 
+  H2,  the  diffuse  character  moves  from  the  hydrogen  center  onto  the  silicon  center.  Thus  a 
balanced  description  of  this  reaction  requires  diffuse  polarization  functions  located  on  all  centers. 

In  addition,  a  diffuse  polarization  function  is  needed  on  the  hydrogen  centers  to  smoothly  permit 
charge  transfer  to  silicon.  In  Table  16  we  give  a  summary  of  the  energetics  of  this  reaction  with 
several  basis  sets  and  approximations.  We  find  that  SiH5’  is  stable  by  ~  20  kcal/mol  relative  to 
SiHa  +  H*,  in  good  agreement  with  the  most  accurate  studies  reported  in  Table  16.  However,  it  is 
clear  that  the  presence  of  a  large  (+39.4  kcal/mol)  reaction  barrier  for  formation  of  SiH3"  +  H2 
rules  out  the  8^2  pathway.  The  energetics  of  the  H"  +  SiHa  reaction  are  illustrated  in  Figure  9. 

Brandemark  and  Siegbahn  (Reference  131)  and  Sheldon,  et  al.  (Reference  127)  have  made 
the  interesting  suggestion  that  the  direct  proton  abstraction  reaction  may  be  the  preferred  pathway: 

H'  +  SiHa  ^  [H--H+-SiH3-]  H2  +  SiH3-  (121) 

They  find  a  much  lower  barrier  (+17.0  kcal/mol)  for  this  direct  abstraction  reaction,  apparently 
ruling  out  the  formation  of  SiH5’  as  an  intermediate.  An  alternate  possible  route  is  to  abstract  H2 
from  the  Cav  pseudorotation  transition  state  geometry  which  lies  only  -  2  kcal/mol  above  the  SiH5' 
minimum  energy  structure  (Reference  132).  At  the  suggestion  of  Mark  Gordon  (Reference  133), 
we  have  examined  dissociation  pathways  staning  from  the  Cay  transition  state  of  SiH5’,  which 
represents  the  intermediate  state  in  the  D3J,  -+  Cay  -» 03^  pseudorotation  process.  As  mentioned, 
the  energy  barrier  to  form  the  Cay  transition  state  structure  is  only  ~  2  kcal/mol,  indicating  that 
dissociation  from  this  geometry  may  be  energetically  possible.  Calculations  to  date,  however,  do 
not  indicate  a  low  barrier  reaction  pathway  from  the  Cay  transition  state.  Both  a  concerted  (C2v) 
and  relaxed  (Cg)  dissociation  lead  to  high  energy  barriers  for  the  removal  of  H2.  We  conclude  that 
the  experimentally  observed  formation  of  SiH3"  must  arise  from  a  proton  transfer  reaction  along 
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Table  16.  Energetics  of  the  H’  +  SiH4  -♦  SiH5“  SiH3’  +  H2  Reaction 


AE  (kca]/mol) 


a 

b 

£ 

d 

£ 

f 

H‘  +  SiH4 

0. 

0. 

0. 

0. 

0. 

0. 

[SiHs  ]  reactants 

0. 

8.6 

0. 

0. 

0. 

0. 

SiHs' 

-  13.8 

-  18.6 

-20.3 

-20.7 

-  10.2 

[SiH5  ]  projiucts 


+  71.7 


+  39.4 


SiH3"  +  H2  -3.1  —  —  --  -9.5  -4.1 


a)  Ref.  127 

b)  Ref.  128 

c)  Ref.  130 
(1>  Ref.  129 
e)  Ref.  131 


HF/6-21G 

CNDO/2;  Basis:  Si:(12s9p2d/6s4p2d],  H:[6s/3s] 

HF/STO-3G  and  several  laiger  basis  sets  developed  from  Si(lOs^)  and  H(5s)  primitive  sets. 
CEPA/Iarge  sp  basis  plus  2  sets  of  d-functions 
Cl/large  sp  basis  plus  2  sets  of  d-functions 


0  Ref.  10  MP4/6-31++G**//MP2/6-3I++G** 
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the  Si-H  bond  to  form  H2  +  SiH3",  rather  than  abstraction  of  H2  from  SiH5"  along  the  Sn2 
reaction  pathway. 

3.4.2  Negative  Ion>MoIecuIe  Reactions:  H'  -1-  SiF4 

The  SN2-reaction  of  H"  +  SiF4  was  studied  in  C3V  symmetry.  The  corresponding  reaction  of 
H'  +  SiHt  is  described  above.  For  that  system,  we  did  not  find  any  barrier  along  the  reaction  path 
in  the  entrance  channel.  This  result  contrasts  with  the  H"  +  CH4  system  which  experimentally 
exhibits  a  large  (~  8.6  kcal/mol)  barrier  toward  reaction.  In  Table  5  we  give  the  enthalpies  for  the 
silane/fluorine  system,  including  the  negative  ion  results  calculated  at  the  G1  and  G2  levels  of 
theory.  It  can  be  seen  that  all  SiF^  (n=l-3)  species  support  a  stable  negative  ion  and  that  SiF3*  is 
particularly  stable  with  an  electron  affinity  of  nearly  -  2.5  eV.  The  reaction; 

H"  +  SiF4  HF  +  SiF3-  (122) 

is  exothermic  by  19.3  kcal/mol.  At  the  MP2/6-31+G*  level  of  theory,  we  find  an  exothermicity  of 
15.0  kcal/mol,  in  fair  agreement  with  experiment,  considering  the  limitations  of  the  basis  set. 

Calculations  of  the  H'  +  SiF4  reaction  surface  \wre  carried  out  in  C3V  symmetry.  Figure  10 
clearly  shows  that  the  reactants  and  products  are  represented  by  two  different  molecular  orbital 
descriptions,  characteristic  of  a  SN2-reaction  pathway.  The  repulsive  character  along  the  entrance 
channel  is  similar  to  that  found  for  the  H"  +  CH4  reaction.  The  barrier  shown  in  Figure  10  is 
clearly  too  high  since  a  two-configuration  MCSCF  or  GVB  description  of  this  system  is  required 
to  carry  the  description  smoothly  from  the  H"  +  SiF4  reactant  channel  to  HF  +  SiF3'  products. 

We  estimate  that  a  multiconfiguration  treatment  will  lower  this  barrier  by  approximately  a  factor  of 
two,  yielding  a  somewhat  larger  barrier,  than  found  for  the  H"  +  CH4  system.  A  noore  detailed 
investigation  of  this  system  is  indicated,  particularly  in  light  of  the  great  stability  of  the  SiF3"  ion. 
Based  on  the  results  of  this  initial  study,  we  predict  that  the  H~  +  SiF4  reaction  would  not  proceed 
for  collision  energies  below  ~  0.8  eV.  At  higher  collision  energies,  the  rate  should  approach 
Langevin  (k  ~  5x10*^  cmVsec)  with  IIF  +  SiF3"  as  the  dominant  products. 

3.4.3  Negative  Ion-Molecule  Reactions:  F*  SiF4 

Studies  of  SiF4  at  low  ion  source  pressures  have  indicated  the  formation  of  several  negative 
ions  (Reference  134).  Recently,  during  the  1992  Gaseous  Electronics  (inference  held  in  Boston, 
Professor  Lawrence  Overzet  inquired  about  the  stability  of  a  possible  silicon  fluoide  anion  of  mass 
1 23.  Such  a  species  was  observed  in  the  sheath  region  of  low  pressure  discharges  in 
silane/fluorine  mixtures  (Reference  135).  The  concentration  of  this  ion  seems  to  track  with  F2" 
concentration.  In  addition,  mass  123  increases  as  mass  85  decreases  which  suggests  the  following 
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X  (bohrs) 


Figure  10.  Reaction  Path  for  HC  +  SiF 


overall  mechanism: 


SiF4  +  F2’  ->  SiFs'  +  F  (123) 

The  species  SiF4*  is  Jahn-Teller  distorted  to  a  much  lower  symmetry  species,  which  can  be 
represented  as  SiF3‘  •  F2.  Previous  studies  by  Marat  and  Janzen  (Reference  136)  have  examined 
the  intemuclear  fluorine  exchange  in  four-,  five-,  and  six-cowdinate  silicon  conqx)unds,  but  the 
structure  of  SiF5~  has  not  been  examined  in  detail. 

Several  possible  reactions  may  lead  to  the  fortnatitm  of  SiF5~.  It  is  assumed  that  the  primary 
ions  in  the  discharge  are  F",  SiF3"  and  possibly  F2~,  formed  by  the  dissociative  attachment 
processes: 


e  +  SiF4  — >  SiF3”  +  F 

(124) 

— >  SiF3  +  F” 

(125) 

— >  SiF2  +  F2 

(126) 

The  last  reaction  is  probably  not  important  since  our  studies  of  SiH[4  showed  that  decomposition  to 
form  SiH2  +  H2  is  energetically  unfavorable. 

In  order  to  evaluate  the  possible  ion-molecule  reactions  leading  to  SiF5~,  we  carried  out  ab 
initio  studies  of  SiF5‘  at  the  HF  and  MP2  levels  of  thetwy  using  several  basis  sets.  This  anion  is 
isoelectronic  with  PF5  and  should  exhibit  thermodynamic  stability  relative  to  SiF3"  +  F2.  Our 
optimized  structure  at  the  highest  level  of  theory,  MP2/6-31 1+G*,  yields  a  trigonal  bypyramid 
structure  with  an  axial  bond  length  of  1 .68  A  and  an  equatorial  bond  length  of  1.65  A.  These  can 
be  compared  with  the  values  1 .58  A  (ax)  and  1.53  A  (eq)  for  the  more  compact  PF5  molecule. 

A  preliminary  evaluation  of  the  thermochemistry  of  SiF5‘  was  carried  out  Fot  the  reaction: 

Si(s)  +  5/2  F2(g)  +  e  SiFj-  (g)  (127) 

we  obtain  AHfO  (0  K)  =  -  504  kcal/mol.  This  value  is  based  on  energies  calculated  at  the 
QCISD(T)/6-3 1 1  +G*//MP2(FC)/6-3 1 1  +<j*  level  of  theory.  This  can  be  compared  with  the 
JANAF  values  for  the  heat  of  formation  of  PF5  of  -  378.5  kcal/mol,  indicating  great  stability  for 
this  pentacoordinated  anion.  The  possible  ion-molecule  reactions  are  indicated  in  Table  17.  We 
.sec  that  reactions  (1),  (2)  and  (5)  are  all  exothermic  leading  to  SiF5”  formation.  The  experimental 
data  of  Overzet  (Reference  135)  indicate  that  reaction  (2)  dominates  since  SiF5"  formation  tracks 
F2‘  depletion.  The  process  represents  F"  abstraction  which  probably  occurs  with  no  barrier  in  the 
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reactant  channels.  Reaction  (5),  the  most  exothermic  for  formation  of  SiF5",  probably  has  a  small 
cross  section  since  decomposition  studies  of  SiH5“  and  NF5  showed  significant  barriers  for 
formation  of  a  diatomic  product 


Table  17.  SiF5'  Ion*Molecuie  Thermochemistry 


(1) 

F'  + 

SiF4 

-> 

SiFg- 

• 

AHr  =  -74.03 

kcal 

(2) 

Fg-  + 

SiF4 

— > 

SiFg- 

+  F  : 

AHr  =  -48.87 

kcal 

(3) 

SiF3"  + 

SiF4 

SiFg- 

+  SiF2  • 

AHr  =  +  4.30 

kcal 

(4) 

SiF-  + 

SiF4 

SiFg- 

+  SI 

AHr  =  -2.65 

kcal 

(5) 

SIF3"  + 

F2 

SiFs' 

AHr  =  -  209.53 

kcal 
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SECTION  4 


CONCLUSIONS 

Several  qualitative  pictures  of  both  ion  and  neutral  fragment  production  in  silane  discharges 
were  developed.  Secondary  ion-silane  reactions  to  form  disilane  ions  were  proposed  to  explain  the 
observed  ion  currents.  A  separate  model  was  proposed  by  Turban,  et  al.  (Reference  7)  in  which 
primary  silane  ions  are  formed  by  low-energy  electron  impact  dissociation,  followed  by  secondary 
ion-molecule  reactions.  Chapman  and  Gallagher  (Reference  5),  however,  have  criticized  the 
applicability  of  this  model  to  cathode  sheath  conditions  where  it  is  believed  that  high-energy 
electron  processes  (10-1(X)  eV)  govern  the  overall  reaction  kinetics.  Their  interpretation  of  the 
discharge  ion  chemistry  is  that  electron  impact  dissociation  of  silane  occurs  as  a  primary  process 
forming  SiHn"’"  (n=0,3),  with  n=2  or  3  heavily  favored  over  more  highly  stripped  silane.  Neutral 
SiHn  radical  production  as  a  primary  process  must  also  be  considered.  Based  on  the  known  rate  of 
H  atom  abstraction  from  silane  and  the  direct  electron  impact  dissociation  reaction,  SiH3  radicals 
are  believed  to  predominate  in  the  discharge  region.  SiH2'*’  and  more  highly  stripped  silane 
cations,  are  probably  driven  to  SiH3  by  fast  reactions  of  the  type  SiHn+  +  SiH4  SiHn+j  + 
SiH3‘‘'.  The  intent  of  this  technical  program  was  to  address  these  uncertainties  in  the  chemistry  of 
silane  ion-molecule  reactions  through  a  theoretical  analysis  of  the  governing  reactions.  Since 
dissociative  reactions,  with  or  without  accompanying  charge  transfer,  may  occur,  our  approach 
was  based  on  detailed  quantum  mechanical  studies  of  the  reaction  surfaces  to  yield  all  possible 
reaction  product  channels.  Subsequent  to  this  analysis  of  the  reaction  energies,  cross  sections  for 
the  various  positive  ion-molecule  processes  were  calculated  within  a  quantum  frameworic.  These 
results  were  analyzed  and  described  in  the  Interim  Report  WP-TR-91-2022  dated  July  1991. 

The  initial  research  effort  carried  out  under  this  program  had  clearly  identified  several  areas 
where  further  theoretical  studies  should  be  carried  out  These  areas  included  a  more  detailed 
analysis  of  ion-molecule  clustering  reactions  and  an  indepth  study  of  the  role  that  negative  ions 
play  in  silane  plasmas.  The  premise  that  positive  ion-molecule  reactions  in  silane  plasmas  follow  a 
characteristic  growth  pattern,  forming  large  cluster  ions  which  eventually  form  observable  dust, 
was  questioned  by  both  experimental  (Reference  11)  and  theoretical  work  (Reference  12).  A 
further  argument  against  such  a  mechanism  is  that  positive  ions  have  short  residence  times  in  the 
discharge  since  they  are  swept  to  boundaries  of  the  reactor  by  the  ambipolar  sheath  fields 
(Reference  137).  Garscadden  (Reference  9),  Haaland  (Reference  138),  and  others  have  observed 
that  such  fields  act  to  constrain  negative  ions  in  the  homogeneous  reactor  region,  thereby, 
amplifying  the  importance  of  their  gas  phase  chemistry. 

As  an  extension  of  the  research  problems  addressed  under  this  technical  program,  additional 
research  was  carried  out  to  gain  a  better  understanding  of  the  role  that  negative  silane  ions  play  in 
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semiconductor  deposition  discharges.  This  research  included;  1)  an  examination  of  the 
thermodynamics  of  SiHn",  Si2Hn“,  SiF„"  and  Si2F„jHn"  anions,  2)  studies  of  sequential 
clustering  reactions  of  silane  anions,  3)  the  role  of  dissociative  attachment  of  electrons  to  silane 
species  forming  radical  neutral  and  anion  products,  and  4)  studies  of  the  formation  and 
decomposition  pathways  for  negative  ion-molecule  reactions  such  as  H'  +  SiHn  and  H"  +  SiFn- 
An  outline  of  our  results  and  conclusions  is  given  below; 

1 .  The  thermochemistry  of  neutral  SiH^  and  SiF„  species,  and  their  corresponding  cations  was 
the  subject  of  several  experimental  (References  139-144)  and  theoretical  studies  (References  145- 
150).  In  contrast,  little  was  known  about  the  corresponding  anions,  particularly  the 
thermochemistry  of  SiFn'  species.  Under  our  current  research  program,  accurate  thermochemistry 
was  developed  for  all  SiHn/SiHn"  and  SiFn/SiFn"  species.  Significant  errors  (>10  kcal)  were 
found  in  the  experimental  data  for  SiF  and  SiF2.  The  present  results  represent  the  most  reliable 
and  consistent  set  of  data  for  these  systems. 

2 .  A  second  research  task  dealt  with  the  analysis  of  sequential  clustering  reactions  among 
negative  silane  ions.  As  discussed  above,  Mandich  (Reference  11)  and  Raghavachari  (Reference 
12)  have  examined  the  Si"*"  -f-  SiH4  clustering  sequence,  leading  to  successive  formation  of 
Si2H2‘*’,  Si3H4'*',  Si4H6'^  and  Si5H2o'*‘*  They  find  that  Si5Hio'*‘  does  not  undergo  H2  elimination 
and,  thereby,  becomes  a  bottleneck  against  growth  of  larger  size  clusters.  The  basic  chemistry  of 
these  clustering  reactions  involved  addition  of  the  SiH4  group  across  the  Si-Si  bond  which  leads 
to  branched  chain  products  rather  than  linear  species.  Studies  of  silicon  anion  clustering  indicate 
that  clustering  reactions  which  represent  H2  elimination  are  not  energetically  favorable.  Several 
silicon  anions  such  as  Si2H2"  and  Si2H4*  indicate  stability  but  do  not  follow  any  systematic 
growth  pattern. 

3 .  A  third  research  task  was  to  examine  dissociative  attachment  of  e  +  SiHt  in  greater  detail. 
Our  preliminary  studies  of  dissociative  attachment  (DA)  of  electrons  to  silane  were  concerned  with 
an  analysis  of  the  possible  product  channels  and  quantum  coupling  with  the  reactants.  The  lowest 
adiabatic  attachment  channel  coaelates  to  SiH2"[2Bi]  +  H2[X*Zg'*'].  The  SiH^’  resonance  state 
corresponding  to  this  channel  appears  to  be  of  the  Feshbach  type  with  a  small  capture  width  from 
reactants.  In  addition,  there  is  a  large  barrier  on  the  reaction  surface  leading  to  SiH2~  +  H2 
products.  The  second  DA  channel  corresponds  to  the  formation  of  SiH3’[^Ai]  +  H[2si.  This 
channel  corresponds  to  a  shape  resonance  but  also  would  appear  to  have  a  small  DA  cross  section 
for  vibrationally  cold  SiH4.  Based  on  these  known  low-energy  quantum  processes,  our 
conclusion  under  the  present  program  was  that  low-energy  direct  electron  dissociative  attachment 
processes  should  have  low  probability. 

The  problem  of  DA  of  e  +  SiH4  at  higher  (>  7  eV)  collisional  energies  was  examined  by 
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looking  at  excited  state  surfaces  of  SiH4".  In  this  case,  where  the  resonance  state  has  a  long 
lifetime,  a  two-step  mechanism  best  describes  the  attachment  process.  Step  one  involves  the 
capture  of  an  electron  into  a  resonance  state  of  the  negative  ion  with  lifetime  so  that  this  resonance 
state  has  forgotten  its  mode  of  formation  before  it  decays.  Step  two  governs  the  rate  of 
dissociation  of  the  resonance  state  into  the  possible  product  channels.  This  dissociation  can  occur 
along  a  number  of  exit  channels,  subject  to  the  usual  conservation  laws  of  total  energy  and  angular 
momentum.  The  competition  among  the  various  possible  modes  of  dissociation,  however,  does 
not  depend  on  the  manner  in  which  the  resonance  state  was  initially  formed.  An  estimate  of  the  DA 
cross  section  was  made  for  e  +  SiH^  at  energies  corresponding  to  threshold  for  these  high-lying 
excited  states.  We  find  cross  sections  in  good  agreement  with  experimental  data  but  fail  to 
correctly  predict  the  relative  branching  to  products. 

An  alternate  mechanism  for  formation  of  negative  ions  by  electron  impact  on  molecules  is  that 
of  polar  dissociation,  a  mechanism  leading  to  ion-pair  product  channels  plus  a  slow  electron. 
Srivastava  (References  119  and  151)  reports  a  large  and  persistent  cross  section  for  formation  of 
the  ion  pairs:  (SiH3-  +  H-*-),  (SiH2-  +  H2+),  (SiH*  +  H+  +  H2),  and  (H’  +  SiH+  +  2H)  for  e  + 
SiH4,  up  to  electron  impact  energies  of  50  eV.  For  this  system  the  ion-pair  channels  are  open  at 
16-18  eV,  in  agreement  with  Srivastava's  data.  The  persistence  of  the  cross  section  up  to  higher 
electron  collision  energies  is  in  agreement  with  a  model  which  couples  the  ion-pair  channel  to  many 
neutral  dissociation  product  channels  which  lie  asymptotically  below  and  exhibit  numerous  curve 
crossings.  An  unobserved  channel  is  open  at  1 1.5  eV  leading  to  H"  +  SiH3'*'  but  the  cross  section 
is  predicted  to  be  very  small  owing  to  the  unlikelihood  of  internal  charge  separation  in  Silij  to 
form  H".  Similar  studies  by  Srivastava  for  e  SiF4  again  lead  to  large  and  persistent  cross 
sections  for  ion-pair  formation.  For  this  system,  however,  the  experimental  data  exhibit 
appearance  energies  for  the  polar  dissociation  process  which  do  not  correlate  well  with  the  new 
thermodynamic  data  in  this  report. 

4.  Finally,  a  study  of  the  formation  of  negative  silane  ions  through  ion-molecule  reactions, 
such  as  H"  +  SiH4,  H"  +  SiF4  and  F"  +  SiF4'  was  carried  out.  Our  initial  studies  of  the  H‘ 

SiH4  reaction  indicated  that  the  most  probable  reaction  pathway  was  abstraction  of  H+  from  Silij 
to  form  H2  +  SiH3".  The  reaction  sequence  ;  H’  +  Siftj  —>  SiH5"  — »  SiH3"  +  H2  in  contrast  to 
the  H’  +  CH4  system,  does  not  exhibit  a  barrier  to  the  formation  of  SiH5".  However,  the 
abstraction  of  H2  occurs  with  a  57  kcal/mol  barrier,  making  the  overall  reaction  improbable.  In  the 
H"  +  SiF4  reaction,  we  find  a  large  barrier  in  the  entrance  channel  leading  to  SiHF4*,  similar  to  the 
corresponding  hydrocarbon  system.  The  topology  of  this  reaction  surface  is  more  complicated 
than  the  SiH5"  system.  For  the  F‘  +  SiF4  reaction,  we  find  that  SiFs",  a  very  stable  anion,  is 
formed.  Our  estimated  heat  of  formation  for  SiF5"  is  -504  kcal/mol.  This  species  has  recently 
been  identified  by  Overzet  (Reference  1 35)  in  a  plasma  discharge. 

In  summary,  a  quantitative  description  of  many  of  the  major  reactions  that  occur  in  the 
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discharge  region,  including  electron  impact  dissociation  and  ionization,  and  the  chemical  role  of 
both  cations  and  anions  (formed  by  dissociative  attachment)  is  required  to  permit  reliable  modeling 
of  rf  discharge  sheaths.  A  goal  of  this  research  program  was  to  identify  the  ion-molecule 
chemistry  occurring  in  silane  discharges,  and  to  predict  both  kinetic  reaction  rates  and  product 
distributions.  Emphasis  was  placed  on  the  development  of  sound  theoretical  methods  to  predict  the 
collisional  energy  dependence  of  the  reaction  rates  and  on  examination  of  the  energetics  and 
kinetics  of  several  ion-molecule  reactions. 
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Patterson  Air  Force  Base,  July  1991. 

2 .  "Electronic  Structure  and  Thermochemistry  of  Silicon  Hydride  and  Silicon  Ruoride 
Anions."  H.  H.  Michels  and  R.  H.  Hobbs,  to  be  published  in  Chemical  Physics 
Letters,  1993. 

3.  "On  the  Bond  Dissociation  Energy  of  SiF  and  SiF'*'."  H.  H.  Michels  and  J.  A. 
Montgomery,  Jr.,  to  be  submitted  to  the  Journal  of  Chemical  Physics,  1993. 

4.  "Dissociative  Electron  Attachment  of  e  +  SiH4."  H.  H.  Michels,  R.  H.  Hobbs  and  J. 
M.  Wadehra,  to  be  submitted  to  Chemical  Physics  Letters,  1993. 

5 .  On  the  Structure  and  Stability  of  the  SiFs"  Anion."  H.  H.  Michels  and  J.  A. 
Montgomery,  Jr.,  to  be  submitted  to  Chemical  Physics  Letters,  1993. 

B.  Presentations 

1 .  "Potential  Energy  Surfaces  for  Silane  Ion-Molecule  Reactions."  H.  H.  Michels  and  R. 
H.  Hobbs,  presented  at  the  Gordon  Research  Conference  on  Plasma  Chemistry,  Tilton, 
New  Hampshire,  August  15-19, 1988. 

2.  "Dissociative  Charge  Transfer  in  H'*'(H3'*')  +  SiH4  Reactions."  H.  H.  Michels  and  R. 
H.  Hobbs,  presented  at  the  1988  41st  Gaseous  Electronics  Conference,  Minneapolis, 
Minnesota,  October  18-21, 1988. 

3 .  "Electronic  Structure  and  Thermochemistry  of  Silicon  Hydride  and  Silane  Fluoride 
Anions."  H.  H.  Michels  and  R.  H.  Hobbs,  to  be  presented  at  the  1991  Annual 


108 


American  Chemical  Society  meeting,  Atlanta,  Georgia,  April  14-19, 1991. 


4.  "Dissociative  Electron  Attachment  of  e  +  SiH4."  H.  H.  Michels  and  R.  H.  Hobbs, 
17th  ICPEAC,  Brisbane,  Australia,  July  10-16, 1991. 

5 .  "Electron  Attachment  to  Silicon  Hydride  and  Silicon  Fluoride  Species."  H.  R  Michels 
and  R.  H.  Hobbs,  45th  Annual  Gaseous  Electronics  Conference,  Boston, 
Massachusetts,  October  27-30, 1992. 


109 


APPENDIX  B 
ABSTRACTS 

APPENDIX  B-1 

POTENTIAL  ENERGY  SURFACES  FOR  SILANE 
ION-MOLECULE  REACTIONS 

U.  H.  Michels  and  R.  H.  Hobbs 
United  Technologies  Research  Center 
East  Hartford,  CT  06108 

ABSTRACT 

Commercial  silane  discharge  reactors  usually  operate  with  a  mixture  of  silane  and  a  noble 
gas.  A  reaction  set  involving  collisions  between  hydrogen  ions,  noble  gas  ions  and  silane  (or 
disilane)  must,  therefore,  be  considered  for  a  complete  description  of  the  ion-molecule  chemistry  in 
such  devices.  Existing  experimental  data  suggest  that  no  simple  reaction  model  for  collisional 
dissociative  charge  transfer  can  be  inferred  for  this  system.  A  theoretical  quantum  mechanical 
study  has,  therefore,  been  initiated  to  define  the  pertinent  potential  energy  reaction  surfaces.  The 
definition  of  these  multidimensional  reaction  surfaces,  through  electronic  structure  calculations, 
should  yield  information  on  the  energetics  and  branching  probabilities  for  the  various  product 
channels. 

We  have  completed  a  study  of  the  reaction  surfaces  for  the  dissociative  charge  transfer 
reaction: 


H+  +  SiH4  SiHn+  +  mH2  +  /H 
(n=0,3;  n+f+2m=5) 

The  dominant  product  ion  is  found  to  be  SiH3''',  formed  from  a  hydride  abstraction  mechanism. 
Dissociation  channels  connecting  to  SiH2'*'  and  Si"*"  exhibit  surprisingly  significant  reaction  barriers 
relative  to  the  SiH3'*'  product  channel.  A  complete  set  of  the  potential  surfaces  for  this  reaction  and 
a  preliminary  analysis  of  the  kinetic  rates  and  product  distribution  are  described. 
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DISSOCIATIVE  CHARGE  TRANSFER  IN  H+(H3+)  +  SiH4  REACTIONS 


H.  H.  Michels  and  R.  H.  Hobbs 
United  Technologies  Research  Center 
East  Hartford,  CT  06108 


ABSTRACT 

An  ab  initio  study  of  the  potential  energy  reaction  surfaces  for  the  dissociative  charge 
transfer  reaction:  H"^  +  SiHa  — >  SiH^"^  +  mH2  +  ^2»  (n=0,3;  n+/+2m=5)  was  undertaken.  The 
basis  set  employed  was  the  standard  6-3 IG*  representation,  augmented  with  diffuse  functions 
(Si:3d,  H;2p)  for  an  accurate  description  of  the  polarizability  of  SiH4.  Dissociative  channels 
connecting  to  SiH2'*'  and  SiH'’’  exhibit,  surprisingly,  significant  reaction  barriers  relative  to  the 
SiH3'*'  product  channel.  The  dissociative  channel  leading  to  Si"*"  from  either  SiH2'''  or  SiH'*’, 
however,  exhibits  no  reaction  barrier.  The  dominant  product  ion  is  found  to  be  SiH3‘*',  formed 
from  a  hydride  abstraction  mechanism.  A  preliminary  calculation  of  the  rate  for  SiH3'^  formation, 
using  a  modified  R-matrix  propagator  technique,  yields  results  close  to  the  Langevin  prediction, 
Kl  =  5.2  X  10’9  cmVsec.  The  chemistry  of  the  H3'*'  +  SiH4  reaction  is  qualitatively  similar  with 
SiH3'*'  +  2H2  formed  as  primary  products. 
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ELECTRONIC  STRUCTURE  AND  THERMOCHEMISTRY  OF  SILICON 
HYDRIDE  AND  SILICON  FLUORIDE  ANIONS 


H.  H.  Michels  and  R.  H.  Hobbs 
United  Technologies  Research  Center 
East  Hartford,  CT  06108 


ABSTRACT 

An  ab  initio  study  of  the  electronic  structure  of  the  SiH^'  and  SiFn"  anions  was  carried  out 
using  M0ller*Plesset  perturbation  theory.  Optimized  geometries  were  calculated  at  the  SCF  and 
MP2  levels  of  theory  using  several  basis  sets.  Correlation  energy  treatments  included  both  MP4 
and  QCISD  methods.  The  anion  thermochemistry  was  calculated  through  several  isogyric 
processes  for  internal  consistency.  Excellent  agreement  with  experimental  electron  affinities  is 
found  for  the  SiHn'  anions.  We  find  that  all  SiFn*  anions  are  stable  with  the  exception  of  SiF4* 
which  exhibits  two  low-lying  resonance  states.  Corresponding  resonance  states  are  also  found  for 
the  SiH4*  anion.  These  states  adiabatically  correlate  to  SiH3*  H  and  SiH3  H*  products  but 
exhibit  a  potential  energy  curve  crossing  at  ~  2.1  A.  The  role  of  these  anions  in  plasma  chemical 
vapor  deposition  of  silane  and  fluorosilane  is  examined. 
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DISSOCIATIVE  ELECTRON  ATTACHMENT  OF  e  +  SiH4 

H.  H.  Michels*.  R.  H.  Hobbs*  and  J.  M.  Wadch/a^ 

•United  Technologies  Research  Center,  East  Hanford,  CT  06108,  USA 
^Waync  State  University.  Detroit.  Ml  48202,  USA 


Two  possible  kinetic  routes,  which  are  dependent 
on  the  collisional  energy,  must  be  considered  for 
dissociative  electron  attachment  of  e  +  SiH4.  The  first 
assumes  that  a  Jahn-  feller  unstable  C^T2l  resonance  state 
of  SiHa"  is  formed  by  primary  electron  capture. 
Subsequently,  this  unstable  negative  ion  can  decay  into 
the  thermodynamically  accessible  product  channels  which 
correspond  to  dissociation  products.  The  second  route 
assumes  that  a  long-lived  resonance  state  of  the  negative 
ion  is  formed  with  a  statistical  decay  into  the  possible 
modes  of  dissociation. 

The  reaction  surfaces  governing  the  tow-energy 
dissociative  attachment  (DA)  of  electrons  to  SiH4  have 
been  analyzed  using  the  intrinsic  reaction  coordinate 
procedure.  In  Figure  1 ,  we  show  the  energetics  along  the 
reaction  coordinate  leading  to  dissociation  of  SiH4  to 
SiH3  +  H.  At  the  HF/6-31-i-fG**  level  of  theory,  we 


Figure  1.  Intrinsic  Reaction  Pathway  for 
Dissociation  of  SiH4  and  SiH4'. 

find  a  bond  dissociation  energy  of  3.27  eV.  At  the 
MP2/6-3I+4-G**  level,  a  bond  energy  of  3.85  eV,  very 


close  to  the  experimental  value,  is  predicted.  Also  shown 
on  Figure  I  arc  preliminary  data  on  the  SiHy'  +  H  and 
SiH3  +  H‘  reaction  surfaces.  We  fuid  an  avoided  crossing 
of  these  two  surfaces,  very  close  to  the  iniernucTear 
separation  where  the  ion  curves  cross  the  neuval  SiH3  + 
H  +  e  surface,  but  probably  in  the  region  where  ihc  ion 
potentials  are  complex.  Thus,  the  excited  state  ion 
channel,  leading  to  SiHy  +  H'  may  depopulate  throjgh 
this  avoided  crossing  leaving  SiH3'  +  H  as  the  only 
significant  products  of  the  DA  reaction.  Experimenul 
data  by  Haaland’  and  by  Potzinger  and  Lampe^  indicate 
that  H"  is  probably  not  formed  by  direct  DA  of  e  +  S1H4. 
These  dissociative  attachment  studies,  including  the  work 
reported  by  Ebinghaus,  et  al^,  all  indicate  that  low  energy 
(<4  eV)  DA  to  SiHa  must  occur  with  a  small  cross 
section  (o<l0'^®cm^).  Our  calculations  of  the  low- 
energy  structure  of  the  SiH4'[^T2]  resonance  suggest  that 
autoionization  is  strongly  favored  over  DA  since  the 
resonarxie  state  geometry  is  very  close  to  that  of  neutral 
SiHs.  Quantum  calculations  of  the  cross  sections  for  low 
energy  DA  of  e  +  SiHa  are  currently  in  progress.^ 

For  higher  energy  collisions,  if  the  resonance 
state  is  long-lived,  a  statistical  break-up  may  occur  wiih 
nearly  equal  probabilities  (a  priori)  for  all  acccs.siblr 
product  channels.  This  alternate  kinetic  route  can  be 
derived  from  the  assumption  of  LTE  (or  at  least  quasi¬ 
steady  state).  In  this  case,  electron  attachment  can  occur 
10  the  silane  fragments  which  result  from  the  statistical 
decomposition  of  the  negative  ion  resonance  state.  The 
product  reaction  channels  are  most  likely  those  correlating 
to  Rydberg  or  high  lying  valence  stales  of  SiH4.  Studies 
of  these  reaction  surfaces  in  the  7-11  eV  region  are  in 
progress. 
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ELECTRON  ATTACHMENT  TO  SILICON  HYDRIDE  AND 
SILICON  FLUORIDE  SPECIES 

H.  H.  Michels  and  R.  H.  Hobbs 
United  Technologies  Research  Center 
East  Hartford,  CT  06108 


ABSTRACT 


The  thermochemistry  of  negatively  charged  species  of  silicon,  hydrogen  and  fluorine  is  at 
present  uncertain,  owing  to  a  sparsity  of  experimental  data  and  the  lack  of  a  systematic  theoretical 
study.  Reliable  heats  of  formation  of  SiHn"  and  SiFn’  species  are  of  considerable  interest  since 
theory  may  help  to  determine  the  role  that  such  anions  play  in  the  deposition  of  amorphous  silicon 
films  and  in  the  etching  of  silicon  surfaces.  In  this  study  a  systematic  evaluation  of  the  geometry 
and  energetics  of  SiH,,'  and  SiF„~  anions,  and  their  corresponding  neutral  parents,  was  carried  out 
using  Moller-Plesset  perturbation  theory.  Optimized  geometries  were  calculated  at  the  SCF  and 
MP2  levels  of  theory  using  several  basis  sets.  Correlation  energy  treatments  included  both  MP4 
and  QCISD  methods.  The  anion  thermochemistry  was  calculated  through  several  isogyric 
processes  and  compared  with  Pople's  G1  and  G2  method  predictions.  Good  agreement  with 
experimental  affinities  is  found  for  the  SiHn"  anions.  We  find  that  all  SiFn"  anions  (n  =  0-3)  are 
stable,  with  relative  stabilities  similar  to  that  found  for  SiHn". 
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ELECTRON  STRUCTURE  AND  THERMOCHEMISTRY  OF 
SILICON  HYDRIDE  AND  SILICON  FLUORIDE  SPECIES 


H.  H.  Michels  and  R.  H.  Hobbs 
United  Technologies  Research  Center 
East  Hartford,  CT  06108 


ABSTRACT 


An  ab  initio  study  of  the  electronic  structure  of  SiH^"  and  SiFn"  anions,  and  the 
corresponding  neutral  molecules,  was  carried  out  using  perturbation  theory.  Correlation  energy 
treatments  include  both  MP4  and  QCISD  methods.  The  anion  thermochemistry  is  calculated 
through  several  isodesmic  processes  and  compared  with  Pople's  G1  and  G2  method  predictions. 
Good  agreement  with  experimental  election  affinities  is  found  for  the  SiH^"  anions.  Significant 
differences  (~10  kcal/mol)  are  found  between  theoretical  and  experimental  heats  of  formation  for 
SiF  and  SiF2.  We  find  that  all  SiF„"  anions  (n  =  0-3)  are  stable,  with  relative  stabilities  similar  to 
that  found  for  SiHn". 
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ON  THE  BOND  DISSOCIATION  ENERGY  OF  SiF  AND  SiF+ 


H.  H.  Michels  and  J.  A.  Montgomery,  Jr. 
United  Technologies  Research  Center 
East  Hartford,  CT  06108 


ABSTRACT 


Ab  initio  calculations  of  the  bond  dissociation  energy  of  SiF  and  SiF*"  were  carried  out  at 
the  QCISD(T)  correlated  level  of  theory.  Several  basis  sets  of  increasing  accuracy  were  explored 
to  examine  the  convergence  properties.  We  find  Dq®  for  SiF  and  SiF+  as  5.95±0.05  eV.  and 
6.7210.05  eV,  respectively.  These  values  are  significantly  larger  than  previous  experimental  and 
theoretical  estimates.  Our  calculated  adiabatic  ionizadon  potential  for  SiF  is  7.3510.05  eV,  in  good 
agreement  with  recent  electron  impact  studies.  These  data  translate  into  a  heat  of  formation  for  SiF 
of  AHf"  (0  K)  =  -12.211.0  kcal/mol,  a  value  in  good  agreement  with  previous  theoretical  studies 
and  ~7  kcal.mol  more  exothermic  than  current  experimental  estimates. 
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DISSOCIATIVE  ELECTRON  ATTACHMENT  OF  e  +  SiH4 
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ABSTRACT 


Dissociative  electron  attachment  of  e  +  SiH4  was  studied  as  a  mechanism  for  forming  the 
negative  fragment  ions:  SiH3",  SiH2",  SiH'  and  Si’  in  silane  discharge  plasmas.  Detailed 
calculations  of  the  low-lying  dissociative  reaction  surfaces  for  SiHt’  indicate  that  low-energy  direct 
attachment  processes  should  have  low  probability.  At  higher  electron  collision  energies  (7-10 
eV),  a  long-lived  resonance  state  of  SiHa",  which  statistically  decomposes  into  the  accessible 
negative  fragment  ions,  appears  to  be  the  dominant  attachment  process. 
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ON  THE  STRUCTURE  AND  STABILITY  OF  THE  SiFj"  ANION 

H.  H.  Michels  and  J.  A.  Montgomery,  Jr. 
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East  Hartford,  CT  06108 


ABSTRACT 


Ab  initio  studies  of  the  electronic  structure  of  SiF5‘  have  been  carried  out  at  the  HF  and 
MP2  levels  of  theory.  A  vibrationally  stable  D3h  structure  is  found  at  both  levels  of  theory,  similar 
to  that  reported  for  PF5.  The  heat  of  formation  of  SiF5*  is  extended  from  energies  calculated  at  the 
QCIST(T)/6-311+G(d)//MP2(FC)/6-311+G(d)  level  of  theory.  We  obtain  AHfO  (0  K)  =  -504 
kcal/mol,  indicating  remarkable  stability  for  this  anion  and  suggesting  its  importance  in 
silicon/fluorine  plasmas. 
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AB  INITIO  STUDY  OF  THE  HEAT  OF  FORMATION  OF  SiF2 
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ABSTRACT 


The  electronic  structure  and  bond  dissociation  energy  of  SiF2  were  examined  at  the 
QCISD(T)  correlated  level  of  theory.  Several  basis  sets  of  increasing  accuracy  were  explored  to 
examine  convergence  properties.  The  heat  of  formation  of  SiF2  was  calculated  from  several 
possible  isodesmic  reactions.  We  find  an  average  value  of  AHf®  (0  K)  =  -152  kcal/mol,  in  good 
agreement  with  previous  theoretical  studies  and  ~  12  kcal/mol  more  exothermic  than  current 
experimental  estimates. 


Supported  in  part  by  WL  under  Contract  F33615-87-C-2718. 
To  be  submitted  to  Chemical  Physics  Letters,  1993. 

119 


II. S.  Government  Printing  Office  750-113 


